Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 83,579 (A)7→8 intergenic (‑50/+43) leuA ← / ← leuL 2‑isopropylmalate synthase/leu operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091383,5791.A85.2% 176.5 / 22.9 54intergenic (‑50/+43)leuA/leuL2‑isopropylmalate synthase/leu operon leader peptide
Reads supporting (aligned to +/- strand):  ref base . (5/3);  new base A (21/25);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

GTTTTTAACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAA‑C‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGTAA  >  NC_000913/83329‑83819
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   
gTTTTTAACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCTCGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGTCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATaaaaaaa                                                                                                                                                                                                                                                      <  2:591947/251‑1 (MQ=255)
      aaCCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTGTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTGCGGGCATaaaaaaaacccg                                                                                                                                                                                                                                                 <  1:609345/250‑6 (MQ=255)
           ggcgggcgATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑t                                                                                                                                                                                                                                         <  2:154868/251‑1 (MQ=255)
            gcgggcgATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑Tg                                                                                                                                                                                                                                        <  1:194798/251‑1 (MQ=255)
                  gATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCgg                                                                                                                                                                                                                                  <  1:610573/251‑1 (MQ=255)
                      gTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGttt                                                                                                                                                                                                                              <  1:489226/251‑1 (MQ=255)
                            aCCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTtttgtt                                                                                                                                                                                                                        >  2:337456/1‑251 (MQ=255)
                                    aaaaTCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCg                                                                                                                                                                                                                <  1:375908/251‑1 (MQ=255)
                                          gCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTgg                                                                                                                                                                                                          >  1:370381/1‑251 (MQ=255)
                                            ccGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCt                                                                                                                                                                                                        >  1:143038/1‑251 (MQ=255)
                                                 gAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGGTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAA‑C‑CCGCGCCA‑TGGCGCGGGGTTTTTGTTTGACTGCGTGCTGGCTTAATgg                                                                                                                                                                                                  >  2:4951/1‑250 (MQ=255)
                                                                     ttCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGACGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTc                                                                                                                                                                               >  2:442905/1‑251 (MQ=255)
                                                                        cATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAA‑C‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTAcc                                                                                                                                                                           <  1:473161/251‑1 (MQ=255)
                                                                                     cccATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaaga                                                                                                                                                                  >  1:691938/1‑248 (MQ=255)
                                                                                        aTACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑t                                                                                                                                                                                                                                         <  1:645488/174‑1 (MQ=255)
                                                                                        aTACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑t                                                                                                                                                                                                                                         >  2:645488/1‑174 (MQ=255)
                                                                                        aTACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGc                                                                                                                                                            <  1:402493/251‑1 (MQ=255)
                                                                                        aTACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGc                                                                                                                                                            <  1:511925/251‑1 (MQ=255)
                                                                                                     ccAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagt                                                                                                                                               <  1:135640/251‑1 (MQ=255)
                                                                                                         cgcAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGAc                                                                                                                                           <  1:659325/251‑1 (MQ=255)
                                                                                                         cgcAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGAc                                                                                                                                           <  2:287454/251‑1 (MQ=255)
                                                                                                           naaTTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAA‑C‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACt                                                                                                                                                                                                                   >  2:305844/2‑176 (MQ=255)
                                                                                                           cAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAA‑C‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACt                                                                                                                                                                                                                   <  1:305844/176‑1 (MQ=255)
                                                                                                               ttGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATa                                                                                                                                      <  2:636220/250‑1 (MQ=255)
                                                                                                                   aGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAGAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCg                                                                                                                                 >  1:173908/1‑251 (MQ=255)
                                                                                                                    gTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCCTGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTAc                                                                                                                                                                            >  2:592960/1‑207 (MQ=255)
                                                                                                                    gTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGTGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTAc                                                                                                                                                                            <  1:592960/207‑1 (MQ=255)
                                                                                                                       tttctttcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAAc                                                                                                                              <  2:697556/250‑1 (MQ=255)
                                                                                                                            ttcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAA‑C‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtg                                                                                                                          <  1:4951/248‑1 (MQ=255)
                                                                                                                              cacaCTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCgn                                                                                                                                                                                        <  2:325854/185‑2 (MQ=255)
                                                                                                                              cacaCTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGc                                                                                                                                                                                        >  1:325854/1‑185 (MQ=255)
                                                                                                                              cacaCTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAg                                                                                                                      <  1:417579/251‑1 (MQ=255)
                                                                                                                                cacTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTc                                                                                                                    >  1:428820/1‑251 (MQ=255)
                                                                                                                                       gCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAAt                                                                                                             <  2:396872/251‑1 (MQ=255)
                                                                                                                                           gCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAg                                                                                                          <  1:207163/250‑1 (MQ=255)
                                                                                                                                                  aCGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTcc                                                                                                      >  1:436464/1‑247 (MQ=255)
                                                                                                                                                    gCCTGTTCACCGTCGCGCACTGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑AATg                                                                                                >  2:187181/1‑250 (MQ=255)
                                                                                                                                                      ntGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAA‑c‑cc                                                                                                                                                                                                                                                 >  2:328484/2‑104 (MQ=255)
                                                                                                                                                      cTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAA‑c‑cc                                                                                                                                                                                                                                                 <  1:328484/104‑1 (MQ=255)
                                                                                                                                                        gTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAAt                                                                                             <  2:110721/249‑1 (MQ=255)
                                                                                                                                                              cGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATg                                                                                     >  2:167025/1‑251 (MQ=255)
                                                                                                                                                               gTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGc                                                                                    >  2:598321/1‑251 (MQ=255)
                                                                                                                                                               gTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTT‑AGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCt                                                                                   <  1:469756/251‑1 (MQ=255)
                                                                                                                                                                      aaTGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCtttt                                                                                >  2:408387/1‑248 (MQ=255)
                                                                                                                                                                       aTGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAA‑AACCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCT‑GATGAATGCGATATGCTTTTag                                                                              >  1:326268/1‑249 (MQ=255)
                                                                                                                                                                                   aaaTAATGACTTCCTGGCTCATGGATTGCGTCCTTGTCTCTTTTAGACCGACTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTACTGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGACGATGCGTTTAGTAGTAGTCGACCGATAATGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCTTTTAGAGTTACTGGATAc                                                                 >  2:608446/1‑250 (MQ=255)
                                                                                                                                                                                    aataatGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGAAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCTTTTAGAGTTACTGGATACa                                                                >  2:232756/1‑250 (MQ=255)
                                                                                                                                                                                                  ggCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCa                                                   <  2:691938/249‑1 (MQ=255)
                                                                                                                                                                                                  ggCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAA‑C‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAAcc                                                >  1:477463/1‑251 (MQ=255)
                                                                                                                                                                                                     tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCt                                              >  1:170250/1‑251 (MQ=255)
                                                                                                                                                                                                      cATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCt                                              <  1:527466/250‑1 (MQ=255)
                                                                                                                                                                                                                gTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATaa                                                                                                                                     >  2:224621/1‑154 (MQ=255)
                                                                                                                                                                                                                gTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGAGGCGTTTAGTAGTAGTAGACCGATaa                                                                                                                                     <  1:224621/154‑1 (MQ=255)
                                                                                                                                                                                                                  ccTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATaaaa                                  >  2:146537/1‑250 (MQ=255)
                                                                                                                                                                                                                             ttAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCa                       <  1:337456/250‑1 (MQ=255)
                                                                                                                                                                                                                              tAGAGCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGc                     >  2:689064/1‑251 (MQ=255)
                                                                                                                                                                                                                                gagCGCCTCGCTTCGGGCATAAAAAAAAC‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCCCTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCACCCCTGACGCAATAAAAACGTCCCGCCAGCGt                   <  1:598321/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cgggcataaaaaaaaa‑a‑CCGCGCAAATGG‑GCGGGTTTTTTGTTTGGCTGCGTGATGGCATAATGCTGGATGCCGCTCCCTCGTCTAACCCGCAAAGGAGATGGGGTTAGGAGGAGGAGACCGATAAAGCGGACGATGTGAGGCATTAAATCAGCTCCAAGA‑GAATGCGATATGCTTTTAGAGTTACTGGATACAAAAAAGGATGTCAACCCTGACGCAATAAAAAAGTCCCCCCAACGTGAGTTCTGCATc       >  2:163879/10‑251 (MQ=255)
                                                                                                                                                                                                                                               ggCATAAAAAAA‑C‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGt    >  1:564858/1‑250 (MQ=255)
                                                                                                                                                                                                                                                cataaaaaaaa‑C‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGTaa  <  2:135452/247‑1 (MQ=255)
                                                                                                                                                                                                                                                  taaaaaaaa‑C‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGTATGCCGCTCACTCGTCTACCGCGCATAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCtt                                                                                  <  2:515970/167‑1 (MQ=255)
                                                                                                                                                                                                                                                  taaaaaaaa‑C‑CCGCGCAA‑TGGCGCGGGTTTTTTGATTGACTGCGTGCTGGCTTAATGCTGGATGCCGCACACTCGTCTACTGCGCAAAGAAGATGCGTTTAGTCGTAGCAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCtt                                                                                  >  1:515970/3‑169 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   
GTTTTTAACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAA‑C‑CCGCGCAA‑TGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCA‑GATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCAATAAAAACGTCCCGCCAGCGTGAGTTCTGCATCCGTAA  >  NC_000913/83329‑83819

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