Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,768,615 G→T 9.0% intergenic (+42/‑50) 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,6150GT 9.0% 179.5 / 10.0 85intergenic (+42/‑50)yfjP/yfjQCP4‑57 prophage; 50S ribosome‑binding GTPase family protein/CP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (35/42);  new base T (5/3);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

GTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCC  >  NC_000913/2768469‑2768754
                                                                                                                                                  |                                                                                                                                           
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:1234470/1‑149 (MQ=255)
                                                                        gTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTg                                                                   >  2:1143362/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)
                                                                                                                                                  |                                                                                                                                           
GTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCC  >  NC_000913/2768469‑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: