Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 32,916 G→A 5.2% M700I (ATG→ATA carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091332,9160GA 5.2% 494.4 / 10.2 156M700I (ATG→ATAcarBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base G (74/73);  new base A (3/5);  total (77/79)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

Only 100 of 157 total aligned reads displayed.

GCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGT  >  NC_000913/32769‑33065
                                                                                                                                                   |                                                                                                                                                     
gcgcgcgcTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGc                                                                                                                                                     <  2:1511760/150‑1 (MQ=255)
  gcgcgcTGGAAGCTGCGGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCgg                                                                                                                                                   <  1:708843/150‑1 (MQ=255)
     cgcTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTag                                                                                                                                                <  1:1866027/150‑1 (MQ=255)
     cgcTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTag                                                                                                                                                >  2:869375/1‑150 (MQ=255)
       gtGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTagag                                                                                                                                              <  2:543693/149‑1 (MQ=255)
       cTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTagag                                                                                                                                              >  1:1022400/1‑150 (MQ=255)
       cTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTagag                                                                                                                                              >  1:896321/1‑150 (MQ=255)
       cTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTagag                                                                                                                                              <  1:66436/150‑1 (MQ=255)
         ggAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGaa                                                                                                                                            <  2:1683276/150‑1 (MQ=255)
         ggAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGaa                                                                                                                                            >  1:1278199/1‑150 (MQ=255)
         ggAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCCATTTCAATAGCGGTaacggt                                                                                                                                            <  1:819562/150‑6 (MQ=255)
         ggAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCCATTTCAATAGCGGTaacggt                                                                                                                                            <  1:1489593/150‑6 (MQ=255)
            aGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGc                                                                                                                                         <  1:246966/150‑1 (MQ=255)
              cggctgGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCCGAAGACCGGGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGTAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGa                                                                                                                                       <  2:64106/148‑1 (MQ=255)
                gctgGCGTACCGGTAATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGaaa                                                                                                                                     <  1:442584/150‑1 (MQ=255)
                 ctgGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAag                                                                                                                                    >  1:447228/1‑150 (MQ=255)
                  tgGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAaga                                                                                                                                   <  2:1279824/150‑1 (MQ=255)
                   ggCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAagag                                                                                                                                  >  1:1706309/1‑150 (MQ=255)
                     cGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGAGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGCGCGTCGGAACCTGAAACAACCGGCGCACGCCCCCGTGACCGCTATTGAAATGGCGGTAGAGACGGCGATAGAGAt                                                                                                                                >  1:298106/1‑150 (MQ=255)
                       tACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTg                                                                                                                              <  2:1305914/150‑1 (MQ=255)
                       tACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTg                                                                                                                              <  1:1306051/150‑1 (MQ=255)
                         ccGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCGTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGc                                                                                                                            <  2:171462/150‑1 (MQ=255)
                          cGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCt                                                                                                                           <  2:903822/150‑1 (MQ=255)
                          cGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCt                                                                                                                           <  2:413928/150‑1 (MQ=255)
                               aTCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGGACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACCACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCg                                                                                                                      >  1:1158837/1‑150 (MQ=255)
                                  ggCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtg                                                                                                                   >  2:1573787/1‑150 (MQ=255)
                                  ggCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtg                                                                                                                   >  2:1205512/1‑150 (MQ=255)
                                     aCCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtggtg                                                                                                                <  2:531131/150‑1 (MQ=255)
                                      ccAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtggtgg                                                                                                               >  1:1077578/1‑150 (MQ=255)
                                        aGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTa                                                                                                             <  1:1413232/150‑1 (MQ=255)
                                              gATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTAcgtccg                                                                                                       <  2:757593/150‑1 (MQ=255)
                                              gATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTAcgtccg                                                                                                       <  1:1580043/150‑1 (MQ=255)
                                                 gCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGATTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCt                                                                                                    >  2:33908/1‑150 (MQ=255)
                                                  cTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCtt                                                                                                   <  2:118034/150‑1 (MQ=255)
                                                  cTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCtt                                                                                                   >  2:610686/1‑150 (MQ=255)
                                                  cTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCtt                                                                                                   >  1:1064049/1‑150 (MQ=255)
                                                   tATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGACACACCCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTa                                                                                                  >  2:513034/1‑150 (MQ=255)
                                                    aTCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTAc                                                                                                 <  2:1757888/150‑1 (MQ=255)
                                                    aTCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTAc                                                                                                 >  1:612710/1‑150 (MQ=255)
                                                      cGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGt                                                                                               >  2:1646173/1‑150 (MQ=255)
                                                       gACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGtt                                                                                              >  2:1571902/1‑150 (MQ=255)
                                                           gTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcg                                                                                          >  2:1641738/1‑150 (MQ=255)
                                                           gTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcg                                                                                          >  2:1535362/1‑150 (MQ=255)
                                                            tGCAGAAGACCGTGAACGCTTCCTGCATGCGGTTGCGCGTCTGAAACTGAAACAACCGGCGAACGCCCCCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGATACCCGCTGGGGTTACGTCCGTCTTACGTTCTcgg                                                                                         >  2:1127326/1‑150 (MQ=255)
                                                            tGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcgg                                                                                         >  1:1413957/1‑150 (MQ=255)
                                                             gCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcggc                                                                                        >  2:695462/1‑150 (MQ=255)
                                                             gCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcggc                                                                                        >  1:1732004/1‑150 (MQ=255)
                                                               agaagaCCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcggcgg                                                                                      <  1:1802106/150‑1 (MQ=255)
                                                               agaagaCCGAGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcggcgg                                                                                      <  1:1809317/150‑1 (MQ=255)
                                                                  agaCCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCg                                                                                   <  2:1300713/150‑1 (MQ=255)
                                                                  agaCCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCGATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCg                                                                                   >  1:952352/1‑150 (MQ=255)
                                                                     ccGTGAACGCTTCCAGCATGCGGTTGAGCGGCTGAAACTGAACCAACCGGCGAACGCCACCGTTACCGCTATTGACATCGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGc                                                                                <  1:347005/150‑1 (MQ=255)
                                                                      cGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCg                                                                               <  2:1205569/150‑1 (MQ=255)
                                                                         gAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATg                                                                            <  1:529355/150‑1 (MQ=255)
                                                                         gAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATg                                                                            <  2:896321/150‑1 (MQ=255)
                                                                          aaCGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATgg                                                                           <  2:234245/150‑1 (MQ=255)
                                                                               ttCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATc                                                                      >  1:1634705/1‑150 (MQ=255)
                                                                                tCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCg                                                                     >  1:1464028/1‑150 (MQ=255)
                                                                                tCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCg                                                                     >  1:582276/1‑150 (MQ=255)
                                                                                  cAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTc                                                                   <  2:746099/150‑1 (MQ=255)
                                                                                  cAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTc                                                                   >  1:615135/1‑150 (MQ=255)
                                                                                     cATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTAt                                                                >  2:1729873/1‑150 (MQ=255)
                                                                                          ggTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGa                                                           >  2:203112/1‑150 (MQ=255)
                                                                                            ttGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGACATTGGCTACGCGATGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATTGAAATCGTCTATGACGAAg                                                         <  2:512530/150‑1 (MQ=255)
                                                                                              gAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGCAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCt                                                       >  1:491614/1‑150 (MQ=255)
                                                                                                     tgaaactgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGc                                                >  1:1304197/1‑150 (MQ=255)
                                                                                                      gaaactgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCg                                               >  2:1049416/1‑150 (MQ=255)
                                                                                                        atctgaaacACCCGGCGAAGGCCACCGTTACCGCTAGTGTAATGGCGGTAGTGAAGGCGTAAGTGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAATTCGTCTATGACGATGCTGACCTGCGTc                                             <  1:1247066/148‑1 (MQ=255)
                                                                                                        aactgaaacAACCGGCGATCGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTc                                             <  2:1022400/150‑1 (MQ=255)
                                                                                                        aactgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTc                                             >  1:725272/1‑150 (MQ=255)
                                                                                                         actgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCg                                            >  1:1181829/1‑150 (MQ=255)
                                                                                                          ctgtaacAACCGGCGAACGCCTCCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGTTAGAGATTGGATACCCGCTGGTGGTACTTACGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGc                                           <  2:400120/150‑1 (MQ=255)
                                                                                                           tgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCt                                          <  1:1099857/150‑1 (MQ=255)
                                                                                                            gaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTa                                         >  2:124298/1‑150 (MQ=255)
                                                                                                             aaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTAc                                        >  2:736439/1‑150 (MQ=255)
                                                                                                                caacCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTc                                     >  2:1461282/1‑150 (MQ=255)
                                                                                                                caacCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGGGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTc                                     <  1:857936/150‑1 (MQ=255)
                                                                                                                 aacCGGCGTACGCCACCGTTACCGCTATTGAATTGGCGGTAGAGACTGCGAAAGAGATTGGCTCCCCGCTGGTGGTACGTCCGTCTTACTGTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTcc                                    <  2:277171/150‑1 (MQ=255)
                                                                                                                 aacCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTcc                                    <  2:464975/150‑1 (MQ=255)
                                                                                                                   ccGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGAAGAAGCAGACCTGCGTCGCTACTTCCAg                                  >  2:304967/1‑150 (MQ=255)
                                                                                                                   ccGGCGAACGCCACCGTTACAGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTAAGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGCCCTGAGTCGCTACTTCCAg                                  >  2:413063/1‑150 (MQ=255)
                                                                                                                    cGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGa                                 >  1:243599/1‑150 (MQ=255)
                                                                                                                       cGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGAcgg                              >  1:521232/1‑150 (MQ=255)
                                                                                                                        gAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGAcggc                             <  2:1722305/150‑1 (MQ=255)
                                                                                                                        gAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGAcggc                             <  1:1535362/150‑1 (MQ=255)
                                                                                                                                tctctACCGCCATTTCAATAGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCg                     >  2:146585/5‑150 (MQ=255)
                                                                                                                                 ctctACCGCCATTTCAATAGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCgt                    >  2:976802/4‑150 (MQ=255)
                                                                                                                                  gTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCgtg                   <  1:1049994/150‑1 (MQ=255)
                                                                                                                                      ccGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTa               <  1:560097/150‑1 (MQ=255)
                                                                                                                                      ccGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTa               <  1:1620288/150‑1 (MQ=255)
                                                                                                                                      ccGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTa               >  1:1028307/1‑150 (MQ=255)
                                                                                                                                        gCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGCCCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAAc             >  2:1066673/1‑150 (MQ=255)
                                                                                                                                        gCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAAc             >  1:1126578/1‑150 (MQ=255)
                                                                                                                                        gCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGATGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCTATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAAc             >  1:1132414/1‑150 (MQ=255)
                                                                                                                                           aTTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGCCCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGAt          >  1:1638932/1‑150 (MQ=255)
                                                                                                                                              gAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATgcg       >  2:1290057/1‑150 (MQ=255)
                                                                                                                                               aaaTGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCGCGGCTGTCGGGCGATGGAAATCGTCTATGACGCAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGCTgcgc      >  2:488800/1‑150 (MQ=255)
                                                                                                                                               aaaTGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGGCTGGCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATgcgc      <  1:1205512/150‑1 (MQ=255)
                                                                                                                                                  tGGCGGTAGAGAAGGCGAACGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAg   <  2:1426555/150‑1 (MQ=255)
                                                                                                                                                   ggCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAgt  >  1:1395951/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGT  >  NC_000913/32769‑33065

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