Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,350,678 Δ1 bp 7.3% noncoding (102/121 nt) arcZ → sRNA positive antisense regulator of rpoS; binds Hfq

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,350,6780T. 7.3% 398.2 / 24.5 137noncoding (102/121 nt)arcZsRNA positive antisense regulator of rpoS; binds Hfq
Reads supporting (aligned to +/- strand):  ref base T (52/73);  new base . (5/5);  total (58/79)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.73e-01

Only 100 of 137 total aligned reads displayed.

AACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACC  >  NC_000913/3350531‑3350823
                                                                                                                                                   |                                                                                                                                                  
    gATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGcc                                                                                                                                             >  2:1543676/1‑150 (MQ=255)
     aTTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCg                                                                                                                                            <  2:1312246/150‑1 (MQ=255)
       tCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCggg                                                                                                                                          <  1:1429009/150‑1 (MQ=255)
        cATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCgggg                                                                                                                                         <  2:134754/150‑1 (MQ=255)
        cATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGCGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCgggg                                                                                                                                         <  1:1095811/150‑1 (MQ=255)
        cATGTAACAAATCATTTAAGTCTTGCTATCTTTACTCCGTGCGGCCTGAAAACCAGTGCTGTGCCCTTGTAACTCATCATAATATTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCgggg                                                                                                                                         <  1:268874/150‑1 (MQ=255)
          tGTATCAAATGATTTAAGTTTTGCTAGCTTAACTGCGTGCGGCCTGTAAATCAGTGCTGTGCCCTTGTAACTCATCTTATTAATTTACGGCGCAGCCAAGATTTCCCTGGGGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtg                                                                                                                                       <  2:1331119/150‑2 (MQ=255)
          tGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTc                                                                                                                                       <  2:554668/150‑1 (MQ=255)
          tGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtg                                                                                                                                       <  1:1018695/150‑2 (MQ=255)
          tGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtg                                                                                                                                       <  2:455168/150‑2 (MQ=255)
          tGTAACAAATCATTTAAGTTTTGCAATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCCTAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGCTAGACCGGGGtg                                                                                                                                       <  1:1138336/150‑2 (MQ=255)
            tAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAt                                                                                                                                     >  1:1561251/1‑150 (MQ=255)
               cAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAtttt                                                                                                                                  >  1:1356517/1‑150 (MQ=255)
                 aaTCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAtttttt                                                                                                                                <  1:905037/150‑1 (MQ=255)
                   tCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAg                                                                                                                              <  2:1452214/150‑1 (MQ=255)
                    cATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGt                                                                                                                             <  1:463701/150‑1 (MQ=255)
                      tttAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTgg                                                                                                                           <  2:245587/150‑1 (MQ=255)
                        tAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCt                                                                                                                         <  2:1578914/150‑1 (MQ=255)
                          aGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCttt                                                                                                                       <  1:415747/150‑1 (MQ=255)
                          aGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCttt                                                                                                                       <  1:679714/150‑1 (MQ=255)
                          aGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGATGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCttt                                                                                                                       <  2:514057/150‑1 (MQ=255)
                           gTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCCCCCCGGTCTAGCCGGGGTCATTTTTTCGTGGCtttt                                                                                                                      >  2:881981/1‑150 (MQ=255)
                           gTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCtttt                                                                                                                      <  1:409008/150‑1 (MQ=255)
                            ttttGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTg                                                                                                                     <  1:73909/150‑1 (MQ=255)
                            ttttGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTg                                                                                                                     <  2:477153/150‑1 (MQ=255)
                            ttttGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTg                                                                                                                     <  2:49510/150‑1 (MQ=255)
                             tttGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGc                                                                                                                    <  2:452987/150‑1 (MQ=255)
                               tGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGcca                                                                                                                  <  1:691932/150‑1 (MQ=255)
                                 cTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccacc                                                                                                                >  2:907773/1‑150 (MQ=255)
                                  aaTCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccaccc                                                                                                               <  1:123509/149‑1 (MQ=255)
                                     cTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                                                                                            >  2:1121728/1‑150 (MQ=255)
                                     cTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                                                                                            >  2:1127609/1‑150 (MQ=255)
                                       tAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCt                                                                                                          >  2:1603061/1‑150 (MQ=255)
                                           tgcgtgcgGCCTGAAAAACAGTGCTTTGCCCTTGTAACTCATCATAATAAATTACGGCGCAGCCAAGATTTCCCTGGTGATGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCa                                                                                                      <  1:1227483/150‑1 (MQ=255)
                                           tgcgtgcgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATTATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCa                                                                                                      <  1:1130790/150‑1 (MQ=255)
                                            gcgtgcgGCCTGAAAAACAGTGCTGGGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCACACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                                                                                     <  1:1468971/150‑1 (MQ=255)
                                               tgcgGCCTGAAAATCAGTGCTGTGCCGTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGAGTTGGCGCAGTATTCGCGCTCCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCAc                                                                                                  <  2:146417/150‑1 (MQ=255)
                                                tcgGCCTGAAAAAGAGTGCTGTGCCGTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACt                                                                                                 <  2:615455/149‑1 (MQ=255)
                                                 cgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACtt                                                                                                >  1:934353/1‑150 (MQ=255)
                                                    ccTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTa                                                                                             <  2:756089/150‑1 (MQ=255)
                                                     cTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAc                                                                                            >  1:1711557/1‑150 (MQ=255)
                                                         aaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcg                                                                                        <  1:1121728/150‑1 (MQ=255)
                                                         aaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcg                                                                                        <  1:1127609/150‑1 (MQ=255)
                                                         aaaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcg                                                                                        >  1:887712/1‑150 (MQ=255)
                                                          aaaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGACAAGATTTCGCTGGTGTTGGCGCAGTATTCGCGCACCCCGGACTAGAAGGGGTCATTTTTTAGTGGCTTTTGCAACCAACGCTTTAAGCTCTTCTCcctcgt                                                                                       >  1:1224311/1‑150 (MQ=255)
                                                           aaCAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTg                                                                                      <  1:1024802/150‑1 (MQ=255)
                                                             cAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGAc                                                                                    <  2:1661476/150‑1 (MQ=255)
                                                             cAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGACACCCACGCTTTCAGCACTTCTACTTCGTGAc                                                                                    >  1:1611253/1‑150 (MQ=255)
                                                              aGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACg                                                                                   >  1:1138106/1‑150 (MQ=255)
                                                                tGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGcc                                                                                 >  2:1164471/1‑150 (MQ=255)
                                                                 gCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCa                                                                                >  1:1717534/1‑150 (MQ=255)
                                                                   tgtgCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCAct                                                                              >  1:1062996/1‑150 (MQ=255)
                                                                   tgtgCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCAct                                                                              >  1:1496657/1‑150 (MQ=255)
                                                                    gtgCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCTGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCCCCCACGCTTTCAGCACTTCTAAGTCGTGACTCCAatc                                                                             >  1:1033515/1‑150 (MQ=255)
                                                                    gtgCCATTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCGACCCACGCTTTCAGCACTTCTACGTCGTGACGCCActc                                                                             <  1:296292/150‑1 (MQ=255)
                                                                       cccTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACtcttc                                                                          >  1:1091206/1‑150 (MQ=255)
                                                                         cTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCACTCTTCttt                                                                       >  1:1521964/1‑150 (MQ=255)
                                                                          ttGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCttt                                                                       <  2:1561251/150‑1 (MQ=255)
                                                                           tGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTCCGTCGTGACGCCACTCTTCTTTc                                                                      >  1:1579112/1‑150 (MQ=255)
                                                                             tAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAt                                                                    <  1:959319/150‑1 (MQ=255)
                                                                                  catcatAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCtt                                                               <  1:855544/150‑1 (MQ=255)
                                                                                      ataataATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAAt                                                           >  1:503262/1‑150 (MQ=255)
                                                                                       taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTc                                                              <  1:1166700/146‑1 (MQ=255)
                                                                                           aatTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCAtt                                                      <  1:551518/150‑1 (MQ=255)
                                                                                               tACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCAcc                                                  <  2:1356517/150‑1 (MQ=255)
                                                                                               tACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCAcc                                                  >  2:1039583/1‑150 (MQ=255)
                                                                                               tACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCGCTTCAATCCATTCAcc                                                  >  2:1247764/1‑150 (MQ=255)
                                                                                                aCGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCg                                                 <  2:55206/150‑1 (MQ=255)
                                                                                                  ggCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGGCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGAc                                               <  2:1510028/150‑1 (MQ=255)
                                                                                                     gcACCCAAGATATCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGCGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGAAGACACTCTTCATTCATCTCGTCAATCCATTCACCGACGtt                                            >  1:193569/1‑150 (MQ=255)
                                                                                                        gCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATc                                         >  2:537666/1‑150 (MQ=255)
                                                                                                        gCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATc                                         <  1:1353901/150‑1 (MQ=255)
                                                                                                         ccAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCt                                        <  2:1496657/150‑1 (MQ=255)
                                                                                                          cAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCtt                                       >  2:943786/1‑150 (MQ=255)
                                                                                                                ttCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAgg                                 <  2:849027/150‑1 (MQ=255)
                                                                                                                  cccTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGcc                               >  2:1419759/1‑150 (MQ=255)
                                                                                                                  cccTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGcc                               >  2:222520/1‑150 (MQ=255)
                                                                                                                    cTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGATGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCgg                             <  1:1193290/150‑1 (MQ=255)
                                                                                                                       gtgtTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAg                          >  2:1089828/1‑150 (MQ=255)
                                                                                                                         gtTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGTTCATTTTTTAGTTGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGt                        <  2:1224311/150‑1 (MQ=255)
                                                                                                                         gtTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGt                        <  1:907773/150‑1 (MQ=255)
                                                                                                                               gcAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGgag                  <  1:1429555/150‑1 (MQ=255)
                                                                                                                               gcAGTATTCGCGCACCCAGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAATGTCAGgag                  >  1:927216/1‑150 (MQ=255)
                                                                                                                               gcAGTAATCGCGCACCCCGGTCTAGCAGGGGTCATTTTTTAGGTGCTTATGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGgag                  <  2:528836/150‑1 (MQ=255)
                                                                                                                                 aGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGAc                <  2:1056601/150‑1 (MQ=255)
                                                                                                                                       aaatgaCCCCGGCTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCCCTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGTTCAGGAGACTGAAtt          >  2:192454/6‑150 (MQ=255)
                                                                                                                                        aatgaCCCCGGCTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAAttt         >  1:513169/5‑150 (MQ=255)
                                                                                                                                         cgcACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTg        <  2:1711557/150‑1 (MQ=255)
                                                                                                                                         cgcACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTg        >  2:1006149/1‑150 (MQ=255)
                                                                                                                                          gcACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGc       <  2:1212345/150‑1 (MQ=255)
                                                                                                                                          gcACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCt      <  2:1521964/150‑1 (MQ=255)
                                                                                                                                          gcACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGGGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGc       <  1:1486450/150‑1 (MQ=255)
                                                                                                                                           gaCCCCGGCTAGACCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACATCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCt      >  2:1087126/2‑150 (MQ=255)
                                                                                                                                           cACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCt      <  1:1405349/150‑1 (MQ=255)
                                                                                                                                             ccccGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGa    >  1:365009/1‑150 (MQ=255)
                                                                                                                                             ccccGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGa    <  2:440059/150‑1 (MQ=255)
                                                                                                                                             ccccGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGa    >  2:182614/1‑150 (MQ=255)
                                                                                                                                              cccGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGAc   <  1:1164471/150‑1 (MQ=255)
                                                                                                                                                  gTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGa    >  1:1626112/1‑145 (MQ=255)
                                                                                                                                                  gTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGa    <  2:1626112/145‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACC  >  NC_000913/3350531‑3350823

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