Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,121,217 A→C 7.9% intergenic (‑262/+25) bssS ← / ← dinI regulator of biofilm formation/DNA damage‑inducible protein I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,121,2170AC 7.9% 239.7 / 9.4 85intergenic (‑262/+25)bssS/dinIregulator of biofilm formation/DNA damage‑inducible protein I
Reads supporting (aligned to +/- strand):  ref base A (38/40);  new base C (4/3);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAAACAGATAAATTTCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACC  >  NC_000913/1121069‑1121365
                                                                                                                                                    |                                                                                                                                                    
ggAAACAGATAAATTTCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGAc                                                                                                                                                     >  2:1248223/1‑150 (MQ=255)
    aCAGATAAATTTCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCgg                                                                                                                                                 <  1:829274/150‑1 (MQ=255)
     cAGATAAATTTCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGc                                                                                                                                                >  2:59062/1‑150 (MQ=255)
       gATAAATTTCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGTCGCCCCGGctt                                                                                                                                              <  1:1001323/150‑3 (MQ=255)
       gATAAATTTCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCaa                                                                                                                                              >  2:1169229/1‑150 (MQ=255)
       gATAAATTTCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTCGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCaa                                                                                                                                              >  2:152374/1‑150 (MQ=255)
         tAAATTTCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCaaaa                                                                                                                                            >  2:1192430/1‑150 (MQ=255)
         tAAATTTCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCaaaa                                                                                                                                            >  1:1141301/1‑150 (MQ=255)
          aaaTTTCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGTCGCCCCGGcttttt                                                                                                                                           <  1:847389/150‑6 (MQ=255)
           aaTTTCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGTCGCCCCGGctttttt                                                                                                                                          <  1:526797/150‑7 (MQ=255)
              ttCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAAt                                                                                                                                       <  1:1169229/150‑1 (MQ=255)
               tCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATc                                                                                                                                      >  2:706694/1‑150 (MQ=255)
               tCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATc                                                                                                                                      <  1:439112/150‑1 (MQ=255)
               tCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATATTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATc                                                                                                                                      <  1:930858/150‑1 (MQ=255)
                    ttttAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGc                                                                                                                                 <  1:1207826/150‑1 (MQ=255)
                    ttttAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTAATCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGc                                                                                                                                 >  1:922148/1‑150 (MQ=255)
                       tAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCata                                                                                                                              >  2:135199/1‑150 (MQ=255)
                           aTTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCAtattat                                                                                                                          >  1:47663/1‑150 (MQ=255)
                              tataTGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCg                                                                                                                       >  1:865891/1‑150 (MQ=255)
                               atatGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCGAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGc                                                                                                                      >  2:831743/1‑150 (MQ=255)
                               ataaGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATAATGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCAGCGGCGACCCGGCAAAAAAAATCAATGCATattattaga                                                                                                                      >  2:1333107/1‑147 (MQ=255)
                                 atGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTg                                                                                                                    <  2:1330720/150‑1 (MQ=255)
                                   gAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGAc                                                                                                                  >  2:19819/1‑150 (MQ=255)
                                    aaCAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACa                                                                                                                 >  1:626619/1‑150 (MQ=255)
                                      cAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACaaa                                                                                                               >  1:947158/1‑150 (MQ=255)
                                      cAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACaaa                                                                                                               >  1:656872/1‑150 (MQ=255)
                                        aTAAAATAGCATGTCATTCATATTTTTTTAGCATATTCTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTAGCTGACAAAcc                                                                                                             >  1:1134949/1‑150 (MQ=255)
                                        aTAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAAcc                                                                                                             <  1:911646/150‑1 (MQ=255)
                                         tAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCa                                                                                                            >  2:804645/1‑150 (MQ=255)
                                          aaaaCAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAg                                                                                                           <  2:1180636/150‑1 (MQ=255)
                                              cAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCAt                                                                                                       >  2:261371/1‑150 (MQ=255)
                                               aGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTCGATGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATc                                                                                                      <  1:152374/150‑1 (MQ=255)
                                                   tGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGgcg                                                                                                  <  2:1141301/150‑1 (MQ=255)
                                                   tGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGgcg                                                                                                  <  1:1248223/150‑1 (MQ=255)
                                                    gTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGgcgc                                                                                                 >  2:140175/1‑150 (MQ=255)
                                                     tcattcatATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCt                                                                                                <  2:1214614/150‑1 (MQ=255)
                                                       attcatATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCttt                                                                                              <  1:549540/150‑1 (MQ=255)
                                                           atatTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGAGACACGGAAAAAAAAATCACTGCATATTATTCGCTGACAAACCACTCATCGGCGAtttaaca                                                                                          >  2:657854/1‑146 (MQ=255)
                                                            tatTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCAc                                                                                         >  1:327368/1‑150 (MQ=255)
                                                               ttttttAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGtt                                                                                      >  1:33635/1‑150 (MQ=255)
                                                                    tAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTg                                                                                 <  2:1431267/150‑1 (MQ=255)
                                                                     aGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTgg                                                                                >  1:655507/1‑150 (MQ=255)
                                                                      gCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGga                                                                               <  1:804645/150‑1 (MQ=255)
                                                                       cATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGgag                                                                              <  2:922148/150‑1 (MQ=255)
                                                                        atatTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTTTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGgaga                                                                             <  1:1347691/150‑1 (MQ=255)
                                                                             gtgCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTc                                                                        >  2:1063304/1‑150 (MQ=255)
                                                                             gtgCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTc                                                                        <  1:140175/150‑1 (MQ=255)
                                                                              tgCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCa                                                                       <  1:135199/150‑1 (MQ=255)
                                                                               gCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCAc                                                                      <  1:1192430/150‑1 (MQ=255)
                                                                               gCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCAc                                                                      <  2:626619/150‑1 (MQ=255)
                                                                               gCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCAc                                                                      <  2:1389462/150‑1 (MQ=255)
                                                                                cAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACt                                                                     <  1:1132599/150‑1 (MQ=255)
                                                                                cAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACt                                                                     >  1:1202843/1‑150 (MQ=255)
                                                                                     attTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATgc                                                                >  2:1259765/1‑150 (MQ=255)
                                                                                      ttttGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCAc                                                                      >  1:1145005/1‑143 (MQ=255)
                                                                                      ttttGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCAc                                                                      <  2:1145005/143‑1 (MQ=255)
                                                                                                gtgtAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTAtct                                                     <  1:971664/150‑1 (MQ=255)
                                                                                                gtgtAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTAtct                                                     <  1:1333107/150‑1 (MQ=255)
                                                                                                   tAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGCCAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTAtgttct                                                  >  2:178755/1‑150 (MQ=255)
                                                                                                     gAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCttt                                                <  1:59062/150‑1 (MQ=255)
                                                                                                        aTTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGTGCGACCCGCCAAAAAAAATAACTGCATATTATTCGCTGACAAAACACTCATCGGCGCTTTCACAAGTTTCCTGGAGAATTACACTAATGCGCTGTTTATCTTCTTTTGt                                             >  2:22321/1‑150 (MQ=255)
                                                                                                        aTTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGt                                             >  2:468424/1‑150 (MQ=255)
                                                                                                        aTTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGt                                             <  2:327368/150‑1 (MQ=255)
                                                                                                            tGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTcgcg                                         <  1:1103899/150‑1 (MQ=255)
                                                                                                            tGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTcgcg                                         <  1:765099/150‑1 (MQ=255)
                                                                                                            tGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTcgcg                                         <  2:656872/150‑1 (MQ=255)
                                                                                                            tGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTcgcg                                         >  1:758280/1‑150 (MQ=255)
                                                                                                              tACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGcc                                       <  2:33635/150‑1 (MQ=255)
                                                                                                              tACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGcc                                       <  1:711627/150‑1 (MQ=255)
                                                                                                               actttaaaTTCGTAAGTAAAATAAAAAAGGCGGGGCGACGCGGCAAAAATAATCACTGCATATTATTGGCTGACAAACCAGTCATCGTCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCa                                      <  2:655507/145‑1 (MQ=255)
                                                                                                                cTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCaa                                     >  1:138887/1‑150 (MQ=255)
                                                                                                                  cgaaaTTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATa                                   <  1:831743/148‑1 (MQ=255)
                                                                                                                      aTTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCg                               <  1:276534/150‑1 (MQ=255)
                                                                                                                            aagtaaAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAat                         >  2:825474/1‑150 (MQ=255)
                                                                                                                              gtaaAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAatta                       >  2:149822/1‑150 (MQ=255)
                                                                                                                                aaaataaaaAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAattatt                     <  1:349936/150‑1 (MQ=255)
                                                                                                                                aaaataaaaAAGCAGAGCCTCCACGACAAAAAAAATCACAGCATATTATTCGATGACAAACCAGTCATCGACGATATCCAACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATATTATTTTGTCGCGCCAATAACCGATAAattatt                     >  2:1227062/1‑150 (MQ=255)
                                                                                                                                      tttttgCCGGGTCGCCCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGACGCGCCAATAACCGATAAATTATTCGCTGc               >  1:941799/6‑150 (MQ=255)
                                                                                                                                       aaaaGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCg              >  2:1224938/1‑150 (MQ=255)
                                                                                                                                        tttgCCGGGTCGCCCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCgg             >  1:1186558/4‑150 (MQ=255)
                                                                                                                                        tttgCCGGGTCGCCCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCgg             >  2:76538/4‑150 (MQ=255)
                                                                                                                                        aaaGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCgg             <  1:19819/150‑1 (MQ=255)
                                                                                                                                         aaGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACTAACCAGTCATCGGCGCTTTCGCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGc            <  1:655169/150‑1 (MQ=255)
                                                                                                                                               gggCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACg      >  2:30969/1‑150 (MQ=255)
                                                                                                                                                   gACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTAcc  <  2:511018/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGAAACAGATAAATTTCTGCTTTTAAGATTTATATGAACAATAAAACAGCATGTCATTCATATTTTTTTAGCATATTGTGCAATTATTTTGAGGAAGTGTAGAAATTTTGTACTCAAAATTCGTAAGTAAAATAAAAAAGCCGGGGCGACCCGGCAAAAAAAATCACTGCATATTATTCGCTGACAAACCAGTCATCGGCGCTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACC  >  NC_000913/1121069‑1121365

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