Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,768,609 T→C 9.7% intergenic (+36/‑56) yfjP → / → yfjQ CP4‑57 prophage; 50S ribosome‑binding GTPase family protein/CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,768,6090TC 9.7% 179.8 / 11.9 85intergenic (+36/‑56)yfjP/yfjQCP4‑57 prophage; 50S ribosome‑binding GTPase family protein/CP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (34/36);  new base C (5/3);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

TCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCC  >  NC_000913/2768461‑2768754
                                                                                                                                                    |                                                                                                                                                 
tCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGGCAGTGGCCCGCGCTGTCTGGGATTGCTGCGTCTGAGTGGTTAAGTCTGCCCTATTATTCTCCGAGCCCGGt                                                                                                                                                   >  1:1112107/1‑149 (MQ=255)
        gTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTg                                                                                                                                           >  1:965812/1‑149 (MQ=255)
               ccTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAgg                                                                                                                                    <  1:367764/149‑1 (MQ=255)
                    ccccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCAAGCCCTGCCGCACTTGCGACAGGGctcg                                                                                                                               <  1:245812/149‑3 (MQ=25)
                     cccTTATTCCCGCCCCGGTGCGGACCATTATTCATGCCGTGCGTTCGTCGGTGGTGGCAGTGGCCCGCGATGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTTCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttc                                                                                                                              <  2:1541434/149‑1 (MQ=37)
                     cccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttc                                                                                                                              <  2:1535332/149‑1 (MQ=255)
                      ccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTATATTCTGCCGTATTATTCTCCGGGCCCTGTCGCTAGTGCGGCAGGTctttct                                                                                                                             >  1:1522547/1‑149 (MQ=25)
                       cTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttctt                                                                                                                            <  1:614393/149‑1 (MQ=255)
                       cTTAATCCCGCCCAGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGGGTTAGTGGCCCGCGCTGTCTGGGTTTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttctt                                                                                                                            <  2:904101/149‑1 (MQ=255)
                              ccGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGt                                                                                                                     >  2:1098971/1‑149 (MQ=255)
                              ccGCCCCGGAGCGGACCATTATTCAGGCCGTGCGGTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGt                                                                                                                     <  1:413485/149‑1 (MQ=255)
                                   ccGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGGCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAc                                                                                                                <  1:144027/149‑1 (MQ=255)
                                    cGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCGTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTAGTTTTAcc                                                                                                               >  1:1386565/1‑149 (MQ=255)
                                    cGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCCAGTGCGGCAGGGCTTTCTTTATTTGTTTTAcc                                                                                                               >  1:1380962/1‑149 (MQ=255)
                                    cGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAcc                                                                                                               <  2:76165/149‑1 (MQ=255)
                                    cGGGGCGGTCCATTATTGAGGCCGGGCGTTGCTGGGTGGTGTCAGTGGCCCGCGCTGTCTGTGATTTCTTCTTGTGTGAGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAcc                                                                                                               <  1:663402/149‑1 (MQ=255)
                                      gagCGGACCATTATACCGCCCGTGCGGTCCGCGGGGGTGGCAGGGGCCCGCGCGGTCTGGGATTGCTTCTTCTGAGGGTTTAATTCTGCCCTATTATTCTCCGCGCCCTGGCGCAAGTGCGGCAGGGCGTTCGTTATTTGTTTTACCCa                                                                                                             <  2:1392790/147‑1 (MQ=255)
                                          ggACCATTATTCAGGACGGGCGTTCCTCGGTGGTGTCTGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACt                                                                                                         <  2:1513941/149‑1 (MQ=255)
                                           gACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTg                                                                                                        >  2:733078/1‑149 (MQ=255)
                                            aCCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTga                                                                                                       >  2:84992/1‑149 (MQ=255)
                                             ccATTATTCAGGGCGTGCGTTCCTCGGTGGTGTCAGTGGCCCACGCTGTCTGGGATTTCTTGTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgag                                                                                                      <  2:1551083/149‑1 (MQ=255)
                                             ccATTAATCAGGCCGGGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCGGTCTGGGATTTCTTCTTCTGAGGGTTTAATTCTGCCCTATTCTGCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgag                                                                                                      <  2:1125659/149‑1 (MQ=38)
                                              cATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGCGCCCTGTCGCAAGTGCGGCAGGGCGTTGTTTATTTGTTTTACCCATACTgagg                                                                                                     >  2:154420/1‑149 (MQ=255)
                                              cATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCACTATTATTATCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACACCTACTgagg                                                                                                     >  2:1389352/1‑149 (MQ=255)
                                                ttattCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgaggag                                                                                                   >  1:792870/1‑149 (MQ=255)
                                                    gcAGGCCGTGCGCTCCTCGGTGGTGGCAGCGGCACTCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTg                                                                                               <  2:175860/148‑1 (MQ=38)
                                                     cAGGCCGTGCGTGCCTCGGGGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGc                                                                                              <  1:720636/149‑1 (MQ=255)
                                                      tgtccGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTTGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCt                                                                                             <  1:1121365/146‑1 (MQ=255)
                                                        gCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGAGTTCTTCTTCTGAGTGTTTAATTCTGCGCGATTTTTCTCCGAGCCCTGGCGCAAGTGCGGCAGGGCTGTCTTTTTTTGTTTTACGCATCCTGAGGAGTCTGCttt                                                                                           >  1:848441/1‑148 (MQ=37)
                                                           gTGCGTTCCTCGGTGGTGGCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGa                                                                                        <  2:1391625/149‑1 (MQ=255)
                                                            tGCGTGCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGAc                                                                                       <  2:1345767/149‑1 (MQ=255)
                                                            ggCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGGCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGAc                                                                                       <  2:381169/148‑1 (MQ=255)
                                                              cGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCAGAGTGTTTAATTCTGCCCTATTAATCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGAccc                                                                                     <  2:170143/149‑1 (MQ=255)
                                                                ttCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGt                                                                                   >  1:67053/1‑149 (MQ=255)
                                                                  ccTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCt                                                                                 >  1:38660/1‑149 (MQ=255)
                                                                        gggggtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGACTGTTAAATTCTGCCCTATTATTCTCCGAGCCCTGGCGCAAGGGCGGCAGGGCTTTCTTTATTTGTTTTCCCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTc                                                                           <  2:505290/145‑1 (MQ=37)
                                                                          ggcgcAGTGGCCCGCGCTGGGTGGGATTTCTTCTTCTGAGTGGTAAATTCTGCCCTATTTTTCTCCGAGCCCGCTCGCACGTGCTGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTcgc                                                                         <  1:715347/145‑1 (MQ=38)
                                                                             tCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCt                                                                      >  2:815060/1‑149 (MQ=255)
                                                                             tCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCt                                                                      >  2:478281/1‑149 (MQ=255)
                                                                               aGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCttt                                                                    >  2:5970/1‑149 (MQ=255)
                                                                               aGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGGGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCttt                                                                    <  1:1129951/149‑1 (MQ=255)
                                                                                gTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTg                                                                   >  2:1143362/1‑149 (MQ=255)
                                                                                gTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTg                                                                   >  2:1234470/1‑149 (MQ=255)
                                                                                  ggCCCGCGCTGGCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGgc                                                                 <  1:1556180/149‑1 (MQ=255)
                                                                                     ccGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGgcgcg                                                              >  2:371266/1‑148 (MQ=255)
                                                                                     ccGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCt                                                              >  1:995757/1‑149 (MQ=255)
                                                                                     ccGCGCTGTCTGGGATTTCTTCTTCTGAGTGAGAAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTGCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCt                                                              <  1:1036588/149‑1 (MQ=255)
                                                                                      cgcgCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCATATTATTCTCCGAGCCCTGGCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTTGCGCTTTGGCGCTg                                                             <  1:1484451/149‑1 (MQ=255)
                                                                                         gcTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCCTTGGCGCTGCaa                                                          <  1:587331/149‑1 (MQ=255)
                                                                                           tGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAAt                                                        <  1:84992/149‑1 (MQ=255)
                                                                                            gTCTTGGTTTTCTTCTTTTGAGTGGTTAATTCTGCCCTTTAATTGTCCTATCCCTTTCTCAAGTGAGGCAGTGCTTTCTTTATTTGTTTTACCCATACTTAGGAGTCTTATTAAGACCCGTCTGGCTTCGCTCTTTGGCGCTGCATATc                                                       <  2:1111080/149‑1 (MQ=255)
                                                                                              cTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATACTCCGAGCCCTGTCTCAAGTGAGGCAGGGCTTTCTTTATTTGTTTTAAACATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTTGCGCTGCAAATCtt                                                     >  1:1266614/1‑149 (MQ=38)
                                                                                               tGGGATTTCTTCTTCTTAGTGTTTAATTCTGACCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTa                                                    <  2:929428/149‑1 (MQ=255)
                                                                                               tGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTa                                                    >  1:1293112/1‑149 (MQ=255)
                                                                                               tGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTATCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGTAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTa                                                    >  2:625864/1‑149 (MQ=255)
                                                                                               ggggATTGCTTCTTCTGAGGGTTTAATTCTGCCCTATTATTCTCCGATCCCTTGCGCAAGTGCGGCAGGGCTTTCTTTATTTGGTTTACCCTTACTGAGGAGTCTGCTTATGAGCCGTCTGGCTTGGCGCTTTGGCGCTGCAAATCTTa                                                    <  1:1152057/148‑1 (MQ=255)
                                                                                                gggAGTTCTCCTTATGAGTGTTAAATTCTGCCCTATTATTCTCCGCGCCCTGTCGCACGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGCGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTAt                                                   <  2:1418081/149‑1 (MQ=37)
                                                                                                   aTTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTATCCGAGCCCTGTCGGAAGTGCGGCAGGGCTTTCTTTATTTGTTATACCCAGTCTGAGGAGTCTGCTTATGACCCCTCTGGCTTCGCGCTTTGTCGCTGCAAATCTTATTCg                                                >  2:541663/1‑149 (MQ=38)
                                                                                                    tttCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGc                                               >  2:783623/1‑149 (MQ=255)
                                                                                                        ttcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCAGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCAAACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTCTCGCTGCAAATCTTATTAGCCGTg                                           >  1:350924/1‑149 (MQ=38)
                                                                                                         tcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGa                                          <  2:866675/149‑1 (MQ=255)
                                                                                                                aggctttAATTCTTCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGCCCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtccg                                   <  2:1522547/145‑1 (MQ=255)
                                                                                                                    tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTaa                               >  1:332298/1‑149 (MQ=255)
                                                                                                                         ttCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGt                          >  1:1142453/1‑149 (MQ=255)
                                                                                                                              cccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAaga                     >  1:337972/1‑149 (MQ=255)
                                                                                                                                cTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTCCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGt                   <  1:554467/149‑1 (MQ=255)
                                                                                                                                 tattatTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGtt                  <  1:825725/149‑1 (MQ=255)
                                                                                                                                 tattatTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGAGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGATAACCCGTGAAGAGtt                  >  2:128123/1‑149 (MQ=255)
                                                                                                                                 tattatTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACACATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGtt                  >  2:361145/1‑149 (MQ=255)
                                                                                                                                  attattCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTACCCCGTGAAGAGTTa                 >  1:727286/1‑149 (MQ=255)
                                                                                                                                  aatattCTCCGAGCCCTTCCGCACTTGCGGCAGGGCTTTCTTTATTTGTTTCACCCATTCTGTGGACTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTa                 <  2:1450408/147‑1 (MQ=21)
                                                                                                                                   ttattCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTAt                <  1:733078/149‑1 (MQ=255)
                                                                                                                                     actcgccGAGCCCTGGCGCCAGTGCGGCAGGGCTGTCTTTATTTGTTTTACCCCTACTGAGGAGGCTGCTTATGCCCCGTCTGGCTTCGCGCGTTGGCGCTGCAAATCTTATTCGCCGTGACCGGCCGTAAACCCGTGAAGAGTTAttt              <  2:1112107/144‑1 (MQ=255)
                                                                                                                                      ttCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTc             <  1:1206912/149‑1 (MQ=255)
                                                                                                                                         tcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTAATTGTTTTACCCATACAGAGGAGTCTGCTTATGACCCGTCTGGCTTCGGGCTTTGGCGATGCAAATCTTATTCGCCGTGCCCGTCCGTTAACCCGTGAAGAGGTATTTcgcg          >  1:1340795/1‑149 (MQ=25)
                                                                                                                                          gaaaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgt         >  2:412816/4‑149 (MQ=21)
                                                                                                                                          ccGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgt         <  2:1275196/149‑1 (MQ=255)
                                                                                                                                           cGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgtg        <  1:1480258/149‑1 (MQ=255)
                                                                                                                                            aaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgtga       >  1:922971/2‑148 (MQ=255)
                                                                                                                                             aGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGATTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgtggt      <  1:712799/149‑1 (MQ=255)
                                                                                                                                             aGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTGGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgtggt      <  2:898470/149‑1 (MQ=255)
                                                                                                                                               cccTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTTGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGc    >  2:795694/1‑149 (MQ=21)
                                                                                                                                                ccTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCAGTTAACCCGTGAAGAGTTATTTCGCGTGGTGcc   >  1:1109380/1‑149 (MQ=255)
                                                                                                                                                ccTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCATATCTTATTCGCGGTGCCCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGcc   >  1:1180338/1‑149 (MQ=21)
                                                                                                                                                ccTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGcc   >  1:184626/1‑149 (MQ=255)
                                                                                                                                                 cTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGccc  >  2:255638/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCC  >  NC_000913/2768461‑2768754

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