Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,454 2 bp→AA coding (97‑98/993 nt) rpoS ← RNA polymerase sigma factor RpoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,4540TA100.0% 168.9 / NA 47Q33M (CAG→TTG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (14/33);  total (14/33)
*NC_0009132,867,4550GA95.7% 158.5 / NA 47Q33M (CAG→TTG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/33);  minor base C (2/0);  total (14/33)
Fisher's exact test for biased strand distribution p-value = 8.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAAC  >  NC_000913/2867211‑2867690
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                           
aGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctac                                                                                                                                                                                                                                       <  2:142599/251‑1 (MQ=255)
       gATTGCCATAACCGCGGGCAATGTTTACCACCAGACGCAAGTCACTCTCGATCATCCGGCGGAGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGcc                                                                                                                                                                                                                                <  2:322585/251‑1 (MQ=255)
                 aCGGCGGGCAATTTTTACCACCAGACGAAAGGTACGCTCGATCAGCCGGCGGCGAGAGGCGACAACTCCACGCAGTGCGCGACGCGCAAAATAAACTTCCTCTTCGGCCGGTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACCCACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCaa                                                                                                                                                                                                                      <  2:244797/251‑1 (MQ=255)
                    gcggGCAATTTTTACCACCAGGCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGGGCGCGATGCGCAAAATAAACTTCTTCTACGGCCGTTAACAGTGGTGAATAACCATTCTCACCATGGTAAAGCTGAGGCGCGTCCAACACTCGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCTCTGGGTTCCAATTCTACTAAGGCCTTGTCGTCaaaaa                                                                                                                                                                                                                   <  1:86940/251‑1 (MQ=255)
                     cggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACTCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAAc                                                                                                                                                                                                                  <  1:326432/251‑1 (MQ=255)
                       ggCAAGTTTTACCACCAGACGCAAGTTACTCTCGATCAGCCGGCGGCGAGAGGCGACATCTCCTCGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCCATCGCACCAAGGGAAAGCTGAGGCGCGTCCAGCACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCt                                                                                                                                                                                                                <  2:452381/251‑1 (MQ=255)
                                       aGACGAAAGTGACTCTCGAGCATAAGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatc                                                                                                                                                                                                <  1:355023/251‑1 (MQ=255)
                                        gACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatca                                                                                                                                                                                               <  2:334022/251‑1 (MQ=255)
                                        gACGCAAGTTACTCGCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatca                                                                                                                                                                                               <  1:242379/251‑1 (MQ=255)
                                             aaGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTc                                                                                                                                                                                          >  1:332013/1‑251 (MQ=255)
                                                  aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAt                                                                                                                                                                                     <  1:329229/251‑1 (MQ=255)
                                                           catcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACGCCGTTCTCATCAAATTCCGCATCTTCAttt                                                                                                                                                                             <  2:407429/250‑1 (MQ=255)
                                                               cggcggCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                                                                                                        <  1:56278/251‑1 (MQ=255)
                                                                            cGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTCACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGc                                                                                                                                                           >  2:329884/1‑251 (MQ=255)
                                                                               cATCTGCACGCAGTGCGCGACGCGCAAAATAAACTTTTTTCTCGGCCGCTACCTGTGGTGAATAACCAATCTCCCCAAGGCAAAGCTGAGTCGCGTCCAACACACTCTGTGTGGCTCCCTGCGATAACAGCTCCTCTTCGGCCAAATCGCTATCCCCCGGTACCACTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTCATCATCAATTTCAGCCTCTTAATTTAAATCATGACCTTTCCGCGt                                                                                                                                                         >  2:162848/1‑250 (MQ=255)
                                                                               cATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGt                                                                                                                                                         >  1:442712/1‑250 (MQ=255)
                                                                               cATCCCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTa                                                                                                                                                        <  2:271816/251‑1 (MQ=255)
                                                                                   cccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCt                                                                                                                                                    <  1:25788/250‑1 (MQ=255)
                                                                                    ccccGCAGTGCGCGAGTCTCCAAAAAAACTTCTTCTTCGGCCGTTAACAGTGGTGATTACCCAATCTCACCTAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGTCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTc                                                                                                                                                     <  1:50640/246‑1 (MQ=255)
                                                                                      aCGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCACATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGAc                                                                                                                                                 >  1:222136/1‑251 (MQ=255)
                                                                                                          aaTAAACTTCTTCTTAGGCCGTTAACAGCGGTGACTAACCAATCTCCCCAAGGTAAAGCTGAGCCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCACATCGTTATCACTTGGTTCCAATTCTACTACTGCCTTTTCGTCAAAACCCTCACCCCCGTTCTCATCAAATTCCGCCTCTTCTTTTAAATCATGAACTTTCAGCGTATTCTAACTCATAAGTTGGCTCCTAccc                                                                                                                             >  2:182877/1‑251 (MQ=255)
                                                                                                             aaaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTg                                                                                                                          <  2:291241/251‑1 (MQ=255)
                                                                                                               aCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGAt                                                                                                                        <  2:408613/251‑1 (MQ=255)
                                                                                                                   cttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTAccc                                                                                                                             <  1:10121/242‑1 (MQ=255)
                                                                                                                   cttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTAccc                                                                                                                             >  2:10121/1‑242 (MQ=255)
                                                                                                                     tcttcGGCCGTAAACAGGGGTTAATAACCAAACTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCGGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCAATGGGATCCAATTCTACTAAGGCCTTTTCGTCAAAAACCGCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCtt                                                                                                                   <  2:254449/250‑1 (MQ=255)
                                                                                                                        tcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTTGGTGGCTCCCTGCGATAACAGTTCCTCGGCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAATCCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGAc                                                                                                                <  2:104412/250‑1 (MQ=255)
                                                                                                                         cGGCCGGTAACAGGGGTGAATAACCAATCTNCCCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACg                                                                                                               <  2:204116/250‑1 (MQ=255)
                                                                                                                              gTTAAAAGTGGTGAATAACCAATCTCTCCAAGGTAAAGCTGATTCGCGTCCAACACACGCTGTGTTGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAAGCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCGCAACTCCGTTCTCATCAAATTCCTCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGGGATCCCTTGACGCAAc                                                                                                           <  2:252622/249‑1 (MQ=255)
                                                                                                                                tAACAGTGGTAAATACACAATCACACAAAGGTAAAGCTGAGTCTAGTCCAACACACGCTGTCTGGCTCACTGCGAAAACAGTTCCTCTTCGGCCACATCGTTATCACTGGGGTCCACTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATAAACTTTCAGCGCATTCTGACTCATAAAGCGGCTCCGACCCGTGCTCCCTTGACGGAACATTc                                                                                                       >  2:287663/1‑251 (MQ=255)
                                                                                                                                    aGTGGTGAATATCCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTCCTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACGCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGc                                                                                                   <  2:386723/251‑1 (MQ=255)
                                                                                                                                       ggtgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCaa                                                                                                 <  1:335566/250‑1 (MQ=255)
                                                                                                                                       ggtgAATAAACAATCTCACCAAGGTAAAGCTGAGGTGCGTACAACCCACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTAATAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCac                                                                                                 <  2:85870/250‑2 (MQ=255)
                                                                                                                                              aaCCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTg                                                                                          <  1:157448/250‑1 (MQ=255)
                                                                                                                                                  aaTCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTcc                                                                                     >  1:2161/1‑251 (MQ=255)
                                                                                                                                                   aTCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTNGACGGAACATTCAAGCAAAAGCCTGGTTccg                                                                                    >  2:203246/1‑251 (MQ=255)
                                                                                                                                                                      gAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCa                                                                 <  2:68964/251‑1 (MQ=255)
                                                                                                                                                                            cgTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAAc                                                           <  2:98386/251‑1 (MQ=255)
                                                                                                                                                                                   cacacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCgg                                                    >  1:228827/1‑251 (MQ=255)
                                                                                                                                                                                    acacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGCTCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCggg                                                   >  1:233013/1‑251 (MQ=255)
                                                                                                                                                                                      ccGCTGTGTGGCTCACAGCAATAACAGTTCCTCTCCGGCCAAAGCATGTTCACTGGGTTCTAATCCTACTAAGGCCTTATCGACAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTTTTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTCGACGGAACATGCCAGGAAAAGCCGTGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCgggct                                                 <  2:179355/250‑3 (MQ=255)
                                                                                                                                                                                           gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTAc                                              <  1:358154/249‑1 (MQ=255)
                                                                                                                                                                                           gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTAc                                              <  1:302920/249‑1 (MQ=255)
                                                                                                                                                                                           gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACgg                                            >  1:367650/1‑251 (MQ=255)
                                                                                                                                                                                                         gATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTa                               >  1:106608/1‑250 (MQ=255)
                                                                                                                                                                                                                   ccTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCCGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCaa                    <  1:324379/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aTCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAAc  <  1:357962/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                           
AGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAAC  >  NC_000913/2867211‑2867690

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