Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,322,396 A→T 5.0% I184N (ATT→AAT)  trpE ← anthranilate synthase subunit TrpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,322,3960AT 5.0% 322.9 / 3.2 97I184N (ATT→AAT) trpEanthranilate synthase subunit TrpE
Reads supporting (aligned to +/- strand):  ref base A (43/49);  new base T (2/3);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTC  >  NC_000913/1322248‑1322545
                                                                                                                                                     |                                                                                                                                                     
atGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAAt                                                                                                                                                      <  2:947980/150‑1 (MQ=255)
 tGCGGCACGGAAACCACTGGCAGCGGCGGCGCGACTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATc                                                                                                                                                     <  2:180109/150‑1 (MQ=255)
    ggCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTTCTGATGGTCAATCAc                                                                                                                                                   >  2:1388697/1‑150 (MQ=255)
         ggAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAg                                                                                                                                             >  2:110020/1‑150 (MQ=255)
         ggAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAg                                                                                                                                             >  1:4345/1‑150 (MQ=255)
          gAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTGATTGACCATCAga                                                                                                                                            <  2:342035/150‑2 (MQ=255)
              ccACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGATGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGttt                                                                                                                                        <  2:12070/150‑1 (MQ=255)
               cACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTGATTGACCATCAgaaaaaa                                                                                                                                       <  1:1384752/150‑7 (MQ=255)
                aCTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTGATTGACCATCAgaaaaaaa                                                                                                                                      <  2:402814/150‑8 (MQ=255)
                aCTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCa                                                                                                                                      <  2:1236885/150‑1 (MQ=255)
                 cTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAg                                                                                                                                     <  1:449979/150‑1 (MQ=255)
                     cAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGGGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCgag                                                                                                                                 <  1:384362/150‑1 (MQ=255)
                      aGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCgaga                                                                                                                                <  1:823842/150‑1 (MQ=255)
                       gcggcggcgCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGAt                                                                                                                               <  1:771958/150‑1 (MQ=255)
                        cggcggcgCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATa                                                                                                                              >  1:1296524/1‑150 (MQ=255)
                        cggcggcgCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATa                                                                                                                              <  1:932849/150‑1 (MQ=255)
                         ggcggcgCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATaa                                                                                                                             <  1:110020/150‑1 (MQ=255)
                            ggcgcggcTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATaaaaa                                                                                                                          <  1:1158886/150‑1 (MQ=255)
                              cgcggcTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACa                                                                                                                        >  2:1497281/1‑150 (MQ=255)
                               gcggcTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAg                                                                                                                       <  2:316392/150‑1 (MQ=255)
                               gcggcTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAg                                                                                                                       <  2:665523/150‑1 (MQ=255)
                                cggcTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGa                                                                                                                      <  1:1082505/150‑1 (MQ=255)
                                     tCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCa                                                                                                                 <  2:50183/150‑1 (MQ=255)
                                      cGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAg                                                                                                                <  2:1296524/150‑1 (MQ=255)
                                       ggTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAgg                                                                                                               <  1:212408/150‑1 (MQ=255)
                                               gCTGACGTAGTTCGTTCAGTCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGTAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTAt                                                                                                       <  2:846085/150‑1 (MQ=255)
                                               gCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTAt                                                                                                       >  2:1375648/1‑150 (MQ=255)
                                                cTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTAtt                                                                                                      >  2:921378/1‑150 (MQ=255)
                                                  gACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTAtttt                                                                                                    >  1:3358/1‑150 (MQ=255)
                                                  gACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTAtttt                                                                                                    >  1:96440/1‑150 (MQ=255)
                                                  gACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTAtttt                                                                                                    >  1:1326781/1‑150 (MQ=255)
                                                   aCGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTc                                                                                                   >  1:327063/1‑150 (MQ=255)
                                                    cGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGTTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTcc                                                                                                  >  1:1525845/1‑150 (MQ=255)
                                                    cGTAGTGCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTcc                                                                                                  <  2:555810/150‑1 (MQ=255)
                                                      tAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCgt                                                                                                <  1:949159/150‑2 (MQ=255)
                                                            gttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAg                                                                                          >  2:1559903/1‑150 (MQ=255)
                                                              tcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGtt                                                                                        >  1:1548313/1‑150 (MQ=255)
                                                               cAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGTTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTg                                                                                       <  2:1525845/150‑1 (MQ=255)
                                                                 ggCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCg                                                                                     >  1:203759/1‑150 (MQ=255)
                                                                      gCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTaaa                                                                                <  1:1559903/150‑1 (MQ=255)
                                                                       cAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTaaa                                                                                <  1:1388697/150‑1 (MQ=255)
                                                                        aGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATc                                                                              >  2:1485332/1‑150 (MQ=255)
                                                                        aGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATc                                                                              >  1:723619/1‑150 (MQ=255)
                                                                         gTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCt                                                                             >  1:1065050/1‑150 (MQ=255)
                                                                         gTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTAGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGTGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCt                                                                             <  2:1602970/150‑1 (MQ=255)
                                                                              aCGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCaaa                                                                        <  1:1020655/150‑1 (MQ=255)
                                                                               cGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAAt                                                                       <  1:879633/150‑1 (MQ=255)
                                                                                gttgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATc                                                                      >  1:474762/1‑150 (MQ=255)
                                                                                  tgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATTAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTc                                                                           <  2:575528/143‑1 (MQ=255)
                                                                                  tgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATGAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTc                                                                           >  1:575528/1‑143 (MQ=255)
                                                                                   gttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCg                                                                   >  1:1348050/1‑150 (MQ=255)
                                                                                    tttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGc                                                                  >  1:1587519/1‑150 (MQ=255)
                                                                                    tttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATAACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCAGCTGACAGTTGCGGTAAATCTTCAAAtaacgc                                                                  >  1:875448/1‑145 (MQ=255)
                                                                                       ttcttcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac                                                               >  2:7483/1‑150 (MQ=255)
                                                                                            cttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGt                                                          <  2:327063/150‑1 (MQ=255)
                                                                                             ttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTc                                                         >  1:374645/1‑150 (MQ=255)
                                                                                              tcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCa                                                        >  1:437609/1‑150 (MQ=255)
                                                                                               cATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCAt                                                       >  2:287563/1‑150 (MQ=255)
                                                                                               cATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCAt                                                       >  2:338611/1‑150 (MQ=255)
                                                                                                aTTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATa                                                      >  1:306571/1‑150 (MQ=255)
                                                                                                  tCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAag                                                    <  1:1375648/150‑1 (MQ=255)
                                                                                                    ggAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAagag                                                  <  2:96440/150‑1 (MQ=255)
                                                                                                       gCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAAc                                               <  2:3358/150‑1 (MQ=255)
                                                                                                        cAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACa                                              <  2:1587519/150‑1 (MQ=255)
                                                                                                              ggCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGgccgc                                        >  2:316705/1‑150 (MQ=255)
                                                                                                               gCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGgccgcc                                       <  2:474762/150‑1 (MQ=255)
                                                                                                                cTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGccgccg                                      >  1:1348363/1‑150 (MQ=255)
                                                                                                                cTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGccgccg                                      <  1:1497281/150‑1 (MQ=255)
                                                                                                                 tGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCga                                     <  1:1104787/150‑1 (MQ=255)
                                                                                                                  ggCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCgaa                                    <  1:1486306/150‑1 (MQ=255)
                                                                                                                   gCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCgaag                                   >  1:375450/1‑150 (MQ=255)
                                                                                                                     cTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCgaagaa                                 >  2:996306/1‑150 (MQ=255)
                                                                                                                     cTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCgaagaa                                 <  2:203759/150‑1 (MQ=255)
                                                                                                                      tGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAAc                                <  2:311146/150‑1 (MQ=255)
                                                                                                                      tGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAAc                                >  2:590601/1‑150 (MQ=255)
                                                                                                                        aaTACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACAt                              <  2:1065050/150‑1 (MQ=255)
                                                                                                                          tACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATgg                            >  1:1608655/1‑150 (MQ=255)
                                                                                                                              ggTGC‑TTTTTTTCTGATGGTCAATCACCATCATCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGATAACAGACCGCCGAAGAACATGGCTtc                        <  2:688825/150‑1 (MQ=255)
                                                                                                                               gTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTtct                       <  2:1326781/150‑1 (MQ=255)
                                                                                                                               gTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTtct                       <  2:1163284/150‑1 (MQ=255)
                                                                                                                                  c‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGt                    <  1:859226/150‑1 (MQ=255)
                                                                                                                                  c‑TTTTTTTCTGATGGTCAAACACCATCAGCGTTTCAGCGAGATAAAAACAGAAATAAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGt                    >  1:818236/1‑150 (MQ=255)
                                                                                                                                         cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttcc              >  2:1385097/3‑150 (MQ=255)
                                                                                                                                          tCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttcct             >  1:915008/1‑150 (MQ=255)
                                                                                                                                          gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAGCATGGCTTCTCGTTcttcct             >  1:911498/2‑150 (MQ=255)
                                                                                                                                           cTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttcctt            >  1:959011/1‑150 (MQ=255)
                                                                                                                                            tGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttccttc           >  2:180025/1‑150 (MQ=255)
                                                                                                                                            tGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttccttc           <  1:608770/150‑1 (MQ=255)
                                                                                                                                            tGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttccttc           >  2:1215873/1‑150 (MQ=255)
                                                                                                                                             gATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCg          <  2:374645/150‑1 (MQ=255)
                                                                                                                                              aTGATCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGc                           <  1:1633307/132‑1 (MQ=255)
                                                                                                                                              aTGATCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGc                           >  2:1633307/1‑132 (MQ=255)
                                                                                                                                               tGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGt        <  1:963724/150‑1 (MQ=255)
                                                                                                                                               tGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGt        >  2:1164017/1‑150 (MQ=255)
                                                                                                                                               tGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTAAAATCCCCCCACAAGGTCATAAGAGAACACGCCGCCGAAGAACATGGCTACTCGTTCTTCCTTCgga        >  1:86365/1‑149 (MQ=255)
                                                                                                                                                ggTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAGGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTa       <  2:4345/150‑1 (MQ=255)
                                                                                                                                                 gTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTAc      <  1:1485332/150‑1 (MQ=255)
                                                                                                                                                     aTCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTc  >  2:741688/1‑150 (MQ=255)
                                                                                                                                                     aTCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTc  >  1:1285888/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGCGGCACGGAAACCACTGGCAGCGGCGGCGCGGCTTCGGTCAGTTGCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGC‑TTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTC  >  NC_000913/1322248‑1322545

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