Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,965,126 T→A 10.0% K10N (AAA→AAT thyA ← thymidylate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,965,1260TA10.0% 626.8 / 5.9 227K10N (AAA→AATthyAthymidylate synthetase
Reads supporting (aligned to +/- strand):  ref base T (59/131);  new base A (28/2);  total (94/133)
Fisher's exact test for biased strand distribution p-value = 5.07e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 233 total aligned reads displayed.

AGCAGTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATGCTGATGTACTGCACCCAGGC  >  NC_000913/2964977‑2965275
                                                                                                                                                     |                                                                                                                                                     
    gTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACttttt                                                                                                                                                   <  1:3764490/150‑1 (MQ=255)
      tCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGc                                                                                                                                                 <  2:595613/150‑1 (MQ=255)
      tCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGc                                                                                                                                                 <  1:3351183/150‑1 (MQ=255)
      tCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGc                                                                                                                                                 <  2:239587/150‑1 (MQ=255)
      tCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGc                                                                                                                                                 <  2:584837/150‑1 (MQ=255)
      tCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGc                                                                                                                                                 <  1:544820/150‑1 (MQ=255)
       cATGGATGATGGAACGCTGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCa                                                                                                                                                <  1:3396813/150‑1 (MQ=255)
       cATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCa                                                                                                                                                <  1:9659/150‑1 (MQ=255)
       cATGGATGATGGAACGCAGGTGGCAACGATTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATCGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCa                                                                                                                                                <  2:2244715/150‑1 (MQ=255)
        aTGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGGGTGCCGTCGGCGAGCACGTTTTGCAt                                                                                                                                               >  2:3106831/1‑150 (MQ=255)
        aTGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTGTGTGTGTGCCCTCGGCGCGCCCGGTTTGccg                                                                                                                                               >  2:418084/1‑148 (MQ=255)
         tGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTAGGTGTTCCTTCGGCGAGCACTTTTTGCATc                                                                                                                                              >  2:925995/1‑150 (MQ=255)
           gatgatGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGTGTGCCTTCGTCGAGCACAATTTGCATCAg                                                                                                                                            >  1:2922391/1‑150 (MQ=255)
            atgatgGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGt                                                                                                                                           <  1:1294840/150‑1 (MQ=255)
             ggatgGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGtt                                                                                                                                          <  2:4000704/149‑1 (MQ=255)
              gatgGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTc                                                                                                                                         <  1:2160577/150‑1 (MQ=255)
                tgGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTGTGTGTTGCTTGGTCGTGCACATTTTTGATAAGTtata                                                                                                                                       >  1:1749932/1‑150 (MQ=255)
                tgGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTGTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTa                                                                                                                                       >  1:225611/1‑150 (MQ=255)
                tgGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGTGTGCCTTCGTCGAGCACATTTTGAATCAGTTCTa                                                                                                                                       >  1:2318905/1‑150 (MQ=255)
                       cAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTg                                                                                                                                <  2:1202076/150‑1 (MQ=255)
                            ggCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCa                                                                                                                           <  2:2427994/150‑1 (MQ=255)
                            ggCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCa                                                                                                                           <  2:2000691/150‑1 (MQ=255)
                             gCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCAt                                                                                                                          <  1:161116/150‑1 (MQ=255)
                              cAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATg                                                                                                                         <  1:1549225/150‑1 (MQ=255)
                                  gTTTCGTGGTCACCAGCGGGCATCCATCTTGCAGGCTAAAACGCATCTGATGACCAAAACTGGAAAGCGTTCCGGTTCCGGTTCGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTc                                                                                                                     <  2:1772544/150‑1 (MQ=255)
                                    ttAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCt                                                                                                                   >  2:1092561/1‑150 (MQ=255)
                                     tAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTc                                                                                                                  <  1:3363983/150‑1 (MQ=255)
                                      aGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAATCGCTTCTGATTACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCa                                                                                                                 <  2:8844/150‑1 (MQ=255)
                                      aGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCa                                                                                                                 <  2:2631019/150‑1 (MQ=255)
                                        ttGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTGGGTGTGGCTTAGTGGAGCACATTTTGAATAAGTTGTAAATACAGTTTCAAGGTTACTCAgg                                                                                                               >  1:4015050/1‑150 (MQ=255)
                                        ttGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTGTGTGTGCCTTCGGGGAGCACTTTTTTCATCAGTTCTAAATAAATTTTAATGGTTGCTCAgg                                                                                                               >  1:1360286/1‑150 (MQ=255)
                                        ttGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTATGTGTGCCGTCGTCGCGCACATTTTGCATCAGTTGTAAATAATGTTTCATGGGTACTCAgg                                                                                                               >  1:1360229/1‑150 (MQ=255)
                                           tCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGTGTGCCGTCGTCGAGCACATTTTGCATCAGTTCTAAATACTGTTTCATGGTTGCACAGGaaa                                                                                                            >  1:2342693/1‑150 (MQ=255)
                                               cAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGTGTGGCGTAGTCGAGCACATTTTTCATAAGTTCTAAATAAAGTTTAATGGTTCCTCAGGAAACctg                                                                                                        >  1:361161/1‑150 (MQ=255)
                                               cAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTATGTGTGGCGTCGTCGAGCACTTTTTTAATAAGTTATAAATAATTTTTGGTGGTTACACAGGaaaactg                                                                                                        >  1:53862/1‑146 (MQ=255)
                                                aGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTGTTTGTGTGTGGCGTCGTGGGGCACTTTTTGAATCAGTTCTAAAAAAAGTTTAATGGTTAATAAGGaaaactgt                                                                                                       >  1:4094278/1‑145 (MQ=255)
                                                 gCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGtt                                                                                                      <  1:1406855/150‑1 (MQ=255)
                                                 gCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTATGTGTGCCTTCGTCGAGCACTTTTTTCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGtt                                                                                                      >  2:3598804/1‑150 (MQ=255)
                                                    ggAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCt                                                                                                   <  2:527886/150‑1 (MQ=255)
                                                       aTCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCCGTTCTAAATAATGTTTCATGGTTCCTCAGGAAACGGGGTTCtgtg                                                                                                >  2:3949475/1‑150 (MQ=255)
                                                       aTCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCtgtg                                                                                                <  1:3124550/150‑1 (MQ=255)
                                                       aTCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCtgtg                                                                                                <  1:456742/150‑1 (MQ=255)
                                                        tCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTgg                                                                                               <  1:860248/150‑1 (MQ=255)
                                                        tCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTATTTGTGGCTTAGTCGAGCACTTTTTTAATCAGTTCTAAATACATTTTAATTGTTCCTCAGGAAACGTGTTTCTGTgg                                                                                               >  1:2851846/1‑150 (MQ=255)
                                                          cATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGCGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCGGTGGGc                                                                                             >  2:3280407/1‑150 (MQ=255)
                                                          cATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTTCTGTGGGc                                                                                             >  2:2552947/1‑150 (MQ=255)
                                                          cATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTAGTCGAGCACTTTTTGCATCAGTTCTAAATAAATTTTCATTGTTGCTCATGAAACGTGTTGCTGTgggg                                                                                             >  2:185043/1‑149 (MQ=255)
                                                          cATCTTGCAGGTTAAAACGCATCTGATGACCAAAAAAGGAAAGCGTTCAGGTTCCGGTACGGTCGTGTTTCAGGGTGCCTAACTCGGGCACATTTTGCATCACTTGTAAAAACATTTTAATGGTTCAACAGGAAAAGGGTTTGTGTGGGc                                                                                             >  2:1671661/1‑150 (MQ=255)
                                                           aTCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTAGGGGTGCCGGAGGCGGGCACTTTTTTGATCCCGTCTAAATACAGTTTGATTGTTCCTCAGGGAACGTGTTTGGGTGGGcg                                                                                            >  2:840936/1‑149 (MQ=255)
                                                           aTCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTGTTATTTGTGCCATCGTCGCGCACATTTTTGAACAGTTCTAAATAAAGTTTAATTGTTGCTCAGGAAAAGAGTTGCTGTggggt                                                                                            >  1:596142/1‑148 (MQ=255)
                                                               tGCAGGTTAAACCGCATCTGATGACCAAAAATGGACAGCGTTCCGGGTCCGGTACGGTCGTTTGTCTGTGTGCCTTCGTCGAGCACTTTGTGCATCAGTTCTAAATACTGGTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGcga                                                                                        <  1:1197614/150‑1 (MQ=255)
                                                               tGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTGATGGTTCCTCAGGTAACGTGTTGCTGTGGGCGGcga                                                                                        >  2:2840507/1‑150 (MQ=255)
                                                                  aGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGcgacga                                                                                     <  1:274828/150‑1 (MQ=255)
                                                                   ggTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGat                                                                                    <  1:2665187/150‑1 (MQ=255)
                                                                       aaaaCGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATAAAGTTTAATTGTTCCACAGGAAACGTGTTGCTGTGGGCTGCGACGATATGc                                                                                >  2:1386193/1‑150 (MQ=255)
                                                                                    tgaCCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATg                                                                   <  2:3340801/150‑1 (MQ=255)
                                                                                    tgaCCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCATCGTCGAGCACTTATTGCAACAGTTCAAAAAAAAGTTTCATGGTTACTAAGGAAACGTGTTGCTGGGGGCTGCGACGATATGCCCAGAAcatcaag                                                                   >  2:1473783/1‑148 (MQ=255)
                                                                                      aCCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGAt                                                                 <  1:3810372/150‑1 (MQ=255)
                                                                                      aCCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGAt                                                                 <  1:1658148/150‑1 (MQ=255)
                                                                                       ccAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATc                                                                <  2:1807767/150‑1 (MQ=255)
                                                                                        cAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATca                                                               <  2:2544863/150‑1 (MQ=255)
                                                                                          aaaaTGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATcaca                                                             <  1:2498997/150‑1 (MQ=255)
                                                                                                 aaaGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGa                                                      <  1:3162887/150‑1 (MQ=255)
                                                                                                 aaaGCGTTCCGGTTCCGGTACGGTCGTGTTTCTGTGTGCCGTCGTCGAGCACTTTTTGCAACAGGTCTAAATAAATTTTCATGGTTGCTCACGAAACCTGTTGCTGTGGGCGGCGACCATATGCCCACACCATCAAGATCACAccccccc                                                      >  2:800196/1‑146 (MQ=255)
                                                                                                     cGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAat                                                  <  1:104193/150‑1 (MQ=255)
                                                                                                        tccggttccggtACGGTCGTTTTTATTTGTTCCTTCGTCGAGCACATTTTTGATCAGGTCTAAAAAAATTTTGATGGTTCCTCAGGAAACGTGGTTCTGTGGGCTGCGACGATATTCCCATACCATCAAGAACACACCCCCCACAaacaa                                               >  2:3760647/1‑149 (MQ=255)
                                                                                                           ggttccggtACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCgg                                            <  1:320033/150‑1 (MQ=255)
                                                                                                           ggttccggtACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCgg                                            <  2:1937373/150‑1 (MQ=255)
                                                                                                           ggttccggtACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCgg                                            <  1:1571863/150‑1 (MQ=255)
                                                                                                           ggttccggtACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCgg                                            <  1:76220/150‑1 (MQ=255)
                                                                                                            gttccggtACGGTCGTTTTTCAGTGTGCCGTCGTCGAGCACATTTTTCATCATTTCTCAAAACATTTTGATGGTTCCGCAGGAAACGTTTTGCTGTGGGCGGCGACGATATGCCCAGACCAACATGATCCAAACCCCCACCATCATCggg                                           >  2:3543931/1‑150 (MQ=255)
                                                                                                             gtccggtACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGa                                          <  1:39711/149‑1 (MQ=255)
                                                                                                              tccggtACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGAt                                         <  2:2922391/150‑1 (MQ=255)
                                                                                                                   tACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGaaa                                    <  1:1872170/150‑1 (MQ=255)
                                                                                                                   tACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGaaa                                    <  2:488371/150‑1 (MQ=255)
                                                                                                                   tACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGaaa                                    <  2:783343/150‑1 (MQ=255)
                                                                                                                    aCGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAg                                   <  1:4184189/150‑1 (MQ=255)
                                                                                                                      ggTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGaa                                 <  1:3579159/150‑1 (MQ=255)
                                                                                                                           tttttCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGc                            <  2:4094278/150‑1 (MQ=255)
                                                                                                                           tttttCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGc                            >  1:2467460/1‑150 (MQ=255)
                                                                                                                             tttCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGccc                          <  2:1039302/150‑1 (MQ=255)
                                                                                                                             tttCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGccc                          <  1:2612200/150‑1 (MQ=255)
                                                                                                                              ttCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGcccc                         >  2:1543712/1‑150 (MQ=255)
                                                                                                                                 tgtgtgCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATg                      <  1:3459270/150‑1 (MQ=255)
                                                                                                                                 tgtgtgCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATg                      >  2:1158857/1‑150 (MQ=255)
                                                                                                                                 tgtgtgCCTTCGTCGAGCACTTTGTGCATGAGGTCTAAATACTGGTTCATGGTTCGTCCGGAACCGTGTTGCTGTGGGCTGCGACGATCTGCCCAGACCCTCATGCACACACCGGCGACACTCATCGGGAGGGACCGACTATGCCCCagg                      >  2:2618218/1‑148 (MQ=255)
                                                                                                                                   tgtgCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATGCt                    >  2:1882274/1‑150 (MQ=255)
                                                                                                                                      gCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATGCTGAt                 <  2:979472/150‑1 (MQ=255)
                                                                                                                                      gCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATGCTGAt                 <  2:3696254/150‑1 (MQ=255)
                                                                                                                                      gCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATGCTGAt                 <  1:3693440/150‑1 (MQ=255)
                                                                                                                                      gCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATGCTGAt                 <  2:2034550/150‑1 (MQ=255)
                                                                                                                                        cTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCAGCGACAATCATCGGGATGGAAAGAATTTGCCCCATGCTGAAGt               >  2:2810261/1‑150 (MQ=255)
                                                                                                                                         ttCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACACTCATCGGGATGGAAAGAATTTGCCCCATGCTGATGTa              >  2:2492954/1‑150 (MQ=255)
                                                                                                                                         ttCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATGCTGATGTa              >  1:896360/1‑150 (MQ=255)
                                                                                                                                         ttCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATGCTGATGTa              >  1:4339520/1‑150 (MQ=255)
                                                                                                                                         ttCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATGCTGATGTa              <  1:1092561/150‑1 (MQ=255)
                                                                                                                                            gtcgAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATGCTGATGTACTg           <  1:734230/150‑1 (MQ=255)
                                                                                                                                            gtcgAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATGCTGATGTACTg           >  1:613821/1‑150 (MQ=255)
                                                                                                                                                   aCTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATGCTGATGTACTGCACCCAg    <  1:2070282/150‑1 (MQ=255)
                                                                                                                                                     tttttGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATGCTGATGTACTGCACCCAGGc  >  2:4323659/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCAGTTCATGGATGATGGAACGCAGGTGGCAACGTTTAGTTGTCACCAGCGGGAATCCATCTTGCAGGTTAAAACGCATCTGATGACCAAAAATGGAAAGCGTTCCGGTTCCGGTACGGTCGTTTTTCTGTGTGCCTTCGTCGAGCACTTTTTGCATCAGTTCTAAATACTGTTTCATGGTTCCTCAGGAAACGTGTTGCTGTGGGCTGCGACGATATGCCCAGACCATCATGATCACACCCGCGACAATCATCGGGATGGAAAGAATTTGCCCCATGCTGATGTACTGCACCCAGGC  >  NC_000913/2964977‑2965275

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: