Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,663,474 A→T 7.0% A11A (GCA→GCT)  suhB → inositol monophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,663,4740AT 7.0% 210.4 / 5.9 71A11A (GCA→GCT) suhBinositol monophosphatase
Reads supporting (aligned to +/- strand):  ref base A (36/30);  new base T (3/2);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

TCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACC  >  NC_000913/2663326‑2663619
                                                                                                                                                    |                                                                                                                                                 
tctttcttCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCAATGCATCCGATGCTGAACATCGCCGTgcgcgca                                                                                                                                                   <  1:294622/149‑1 (MQ=255)
                gATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGt                                                                                                                                    >  4:45299/1‑148 (MQ=255)
                       cACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaa                                                                                                                             <  1:28795/148‑1 (MQ=255)
                        aCGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaatt                                                                                                                           >  7:182231/1‑149 (MQ=255)
                        aCGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaat                                                                                                                            <  3:416/148‑1 (MQ=255)
                        aCGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaat                                                                                                                            <  3:154242/148‑1 (MQ=255)
                            aTCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGAGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGcc                                                                                                                        >  2:23641/1‑148 (MQ=255)
                                  tCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAAc                                                                                                                  >  7:131619/1‑148 (MQ=255)
                                    aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGc                                                                                                                                        <  4:141381/124‑1 (MQ=255)
                                    aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGc                                                                                                                                        >  3:141381/1‑124 (MQ=255)
                                    aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGc                                                                                                                                        <  8:94762/124‑1 (MQ=255)
                                    aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGc                                                                                                                                        >  7:94762/1‑124 (MQ=255)
                                    aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                               <  1:138147/149‑1 (MQ=255)
                                     tAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                               >  2:368272/1‑148 (MQ=255)
                                          cGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAcc                                                                                                         <  3:34946/149‑1 (MQ=255)
                                           gAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAccc                                                                                                        >  7:90226/1‑149 (MQ=255)
                                            aaGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAcccc                                                                                                       <  1:377845/149‑1 (MQ=255)
                                                    tGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGGAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTg                                                                                               >  5:54775/1‑149 (MQ=255)
                                                        gATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGa                                                                                           <  8:131619/149‑1 (MQ=255)
                                                        gATCTCGCAGCTATACTCTGCGCCGTTTTCCCGTTCATTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGCTAATTTAATTGCCAAAAACTATGAAACCCAGGACCCTCTAg                                                                                            >  3:50918/1‑148 (MQ=255)
                                                              gctgctATAATCTGCGCCGTGTTCCCGTTCTTTAACATCCAGTGCGAGAGACCGATGCATCAGATGCTGAACATCGCCGGGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCTAAAAACTATGAAACCCCGGACGCTGTAGAAGCGa                                                                                      <  5:192927/148‑1 (MQ=255)
                                                                   tataCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAg                                                                                 >  7:235088/1‑148 (MQ=255)
                                                                     taCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGaaa                                                                              <  1:347292/149‑1 (MQ=255)
                                                                     taCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGaaa                                                                              <  8:69534/149‑1 (MQ=255)
                                                                     taCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGa                                                                                >  3:103736/1‑147 (MQ=255)
                                                                        tctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAgg                                                                            >  4:266760/1‑148 (MQ=255)
                                                                        tctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCAAGAAAgg                                                                            >  6:344506/1‑148 (MQ=255)
                                                                          tGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAc                                                                                                          >  8:337740/1‑116 (MQ=255)
                                                                          tGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAc                                                                                                          <  7:337740/116‑1 (MQ=255)
                                                                          tGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGc                                                                           <  8:182231/147‑1 (MQ=255)
                                                                          tGCGCCGTTTTCACGTTCTTTAACATCCAGTGAGAGATACCGATGCATCCGATGCTGAACATCGCCGTGCGAGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGACAGAAAGCCAg                                                                         >  1:122004/1‑149 (MQ=255)
                                                                            cgcCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCTATGCATCCTATGCTGAACATCGCCGTGCGCGCAGCGAGCAAGGCGGGTAATTTAATTGCCAAATACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCa                                                                          <  1:209202/146‑1 (MQ=255)
                                                                             gcCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTa                                                                                                                >  4:123753/1‑107 (MQ=255)
                                                                                     cccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAACACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCg                                                              >  7:301603/1‑149 (MQ=255)
                                                                                        gTTCTTTAACATCCAGTGAGAGAGACCGATTCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                           <  7:177032/149‑1 (MQ=255)
                                                                                                aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGt                                                    <  1:116680/148‑1 (MQ=255)
                                                                                                aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa                                                   <  1:23641/149‑1 (MQ=255)
                                                                                                aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa                                                   <  7:339060/149‑1 (MQ=255)
                                                                                                aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa                                                   >  1:379206/1‑149 (MQ=255)
                                                                                                  aTCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGa                                                                                                             <  6:92670/89‑1 (MQ=255)
                                                                                                  aTCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGa                                                                                                             >  5:92670/1‑89 (MQ=255)
                                                                                                  aTCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCATAAAGTCAGTAACGATTTCGTGACCAACGTAg                                                  >  5:497/1‑148 (MQ=255)
                                                                                                  aTCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCCGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACg                                                     >  7:57633/1‑145 (MQ=255)
                                                                                                   tCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGa                                                 >  1:296245/1‑148 (MQ=255)
                                                                                                        tGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAAATATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGc                                           <  3:45299/149‑1 (MQ=255)
                                                                                                          agagagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCt                                          <  3:266760/148‑1 (MQ=255)
                                                                                                          agagagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTg                                         >  7:81551/1‑149 (MQ=255)
                                                                                                          agagagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTAGTGACCAACGTAGATAAAGCTg                                         >  8:134701/1‑149 (MQ=255)
                                                                                                           gagagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGc                                        >  7:199085/1‑149 (MQ=255)
                                                                                                            agagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGCCGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGAGTTCGTGACCACCGTAGATAAGGCTGcc                                       >  8:170554/1‑149 (MQ=255)
                                                                                                              agaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCg                                      >  6:174307/1‑148 (MQ=255)
                                                                                                                  cGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGACACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCg                                 >  2:109174/1‑149 (MQ=255)
                                                                                                                  cGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTTGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCg                                 >  4:70776/1‑149 (MQ=255)
                                                                                                                    aTGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGt                               >  8:38576/1‑149 (MQ=255)
                                                                                                                        aTCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTATTTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                           <  4:166889/117‑1 (MQ=255)
                                                                                                                        aTCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                           >  3:166889/1‑117 (MQ=255)
                                                                                                                         tCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTa                          >  7:283006/1‑149 (MQ=255)
                                                                                                                           cGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTAt                         <  7:57289/148‑1 (MQ=255)
                                                                                                                           cGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATc                        <  8:90226/149‑1 (MQ=255)
                                                                                                                                 tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa                                                   >  6:338151/1‑116 (MQ=255)
                                                                                                                                 tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa                                                   <  5:338151/116‑1 (MQ=255)
                                                                                                                                  gAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                           <  7:56203/107‑1 (MQ=255)
                                                                                                                                  gAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                           >  8:56203/1‑107 (MQ=255)
                                                                                                                                   aaCATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTATTTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGATAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGa                 <  4:354699/148‑1 (MQ=255)
                                                                                                                                   aaCATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTa                          <  1:217977/139‑1 (MQ=255)
                                                                                                                                   aaCATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGa                 <  2:296245/148‑1 (MQ=255)
                                                                                                                                   aaCATCGCCGTGCCCGCAGCGCGCAAGGAGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCAGTATAACCCAGCCAGAAAGGCACTAACGAATTCGTGAAAAACGTAGATAAAGCTGACGAAGCGTTGATTATAGatcgcaa                >  2:217977/1‑143 (MQ=255)
                                                                                                                                     cATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTc              <  8:235088/149‑1 (MQ=255)
                                                                                                                                      aTCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAAATATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCg             >  4:285905/1‑149 (MQ=255)
                                                                                                                                       ccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAATGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCg             >  4:289596/2‑148 (MQ=255)
                                                                                                                                         gCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTTGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTa           <  3:70776/148‑1 (MQ=255)
                                                                                                                                            gTGCGCGCAGCGCGCAAGGCGGGTAAGTTAATTGCCATAAACTATGAAACCCCGGACGCGGGAGAACCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGAGTCGTAAATc       >  5:110559/1‑149 (MQ=255)
                                                                                                                                               cgctcagcgcgcaAGGCGGGTAATTTAATTGCCAATAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTTGATAAAGCTGCCGAAGCGGTGATTATCGTCAGGATTCGTAAATCTTa    <  4:50918/149‑1 (MQ=255)
                                                                                                                                               cgcgcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATACAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTa    <  2:379206/149‑1 (MQ=255)
                                                                                                                                               cgcgcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTa    <  1:368272/149‑1 (MQ=255)
                                                                                                                                                gcgcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTAt                         <  5:91549/127‑1 (MQ=255)
                                                                                                                                                gcgcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTAt                         >  6:91549/1‑127 (MQ=255)
                                                                                                                                                gcgcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTAc   <  8:81551/149‑1 (MQ=255)
                                                                                                                                                  gcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTAcc  >  5:327792/1‑148 (MQ=255)
                                                                                                                                                  gcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTAcc  >  2:144828/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATCTTACC  >  NC_000913/2663326‑2663619

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑