Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 580,697 T→C 5.8% intergenic (‑252/‑137) ylcI ← / → nohD DUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,6970TC 5.8% 182.6 / 3.3 69intergenic (‑252/‑137)ylcI/nohDDUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base T (38/26);  new base C (2/2);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

TATTTGACGTGGTTTGATGGCGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACA  >  NC_000913/580548‑580846
                                                                                                                                                     |                                                                                                                                                     
tattTGACGTGGTTTGATGGCGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGgt                                                                                                                                                       <  1:184045/150‑1 (MQ=255)
 attTGACGTGGTTTGATGGCGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCAGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGgtg                                                                                                                                                      >  2:706455/1‑150 (MQ=255)
   ttGACGTGGTTTGATGGCGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGtt                                                                                                                                                    <  2:1205459/150‑1 (MQ=255)
    tGACGTGGTTTGATGGCGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGttt                                                                                                                                                   >  1:703063/1‑150 (MQ=255)
     gACGTGGTTTGATGGCGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTc                                                                                                                                                  >  1:421102/1‑150 (MQ=255)
         tGGTTTGATGGCGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGCGTTTCCGtt                                                                                                                                              >  2:96816/1‑150 (MQ=255)
         tGGTTTGATGGCGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCCATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGCGTTTCCttt                                                                                                                                              >  2:96818/1‑150 (MQ=255)
          ggTTTGATGGCGTAGATGCACGTTGTGACATGTAGATGATAATTATTAGCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGttc                                                                                                                                             <  1:98232/150‑1 (MQ=255)
            tttGATGGCGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGttctt                                                                                                                                           <  1:575540/150‑1 (MQ=255)
              tGATGGCGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGttcttct                                                                                                                                         >  2:1113845/1‑150 (MQ=255)
              tGATGGCGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCAGttcttct                                                                                                                                         >  1:1180972/1‑150 (MQ=255)
               gATGGCGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCTGTTTAATTTGTTTCCTttcttctt                                                                                                                                        >  2:810963/1‑150 (MQ=255)
               gATGGCCTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGttcttctt                                                                                                                                        <  2:232858/150‑1 (MQ=255)
                aTGGCGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGTCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGttcttcttc                                                                                                                                       <  1:170071/150‑1 (MQ=255)
                    cGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTtcgtcg                                                                                                                                   <  2:336450/150‑1 (MQ=255)
                     gTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTcgtcgt                                                                                                                                  >  1:621144/1‑150 (MQ=255)
                        gATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCATAAc                                                                                                                               >  1:705496/1‑150 (MQ=255)
                            cACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTaa                                                                                                                           >  1:260859/1‑150 (MQ=255)
                            cACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTaa                                                                                                                           >  2:552915/1‑150 (MQ=255)
                            cACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCAGATCTTTTTCGTCGTAACTTaa                                                                                                                           >  2:541673/1‑150 (MQ=255)
                                  gtgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTAACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtatt                                                                                                                     <  2:31887/150‑1 (MQ=255)
                                      cctGTAGCTGATAATTATTATGATTTTGAGGGTCCTTTCAGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGTGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtatttatt                                                                                                                 <  2:1313169/148‑1 (MQ=14)
                                       aTGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTatttattt                                                                                                                >  1:799286/1‑150 (MQ=255)
                                        tGTAGTTGATAATTATTATGATTTTTTGCGTCCTTTCCGGAGATCCGACAGTTTACTGGGCGGAGTCCTCGCTGGTTTTCGCTATTTTTGAAAATTTTCAGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTAtttattta                                                                                                               <  1:121796/150‑1 (MQ=14)
                                         gTAGATGATAATTATTATCATTTTGTGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaaa                                                                                                             >  2:689066/1‑150 (MQ=255)
                                           aGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaaaa                                                                                                            >  2:50892/1‑150 (MQ=255)
                                           aGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaaaa                                                                                                            >  2:1268300/1‑150 (MQ=255)
                                            gattatAATTATTATCATTTTGTGGGTCCTTTCAGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTGTTCGCTATTTATGAAAATTTTCAGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAAt                                                                                                           >  2:560536/1‑150 (MQ=255)
                                                 tAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGTGGCGTCGACCTCGCGGGTTTTATCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTAGTCGTAACTTAATGTATTTATTTAAATTAcccc                                                                                                      >  1:94726/1‑150 (MQ=17)
                                                   attattATCATTTTGGGGGTCCTTTCCGGCGATCCGACAGGGTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGGTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCt                                                                                                    <  2:985415/150‑1 (MQ=255)
                                                         atCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaag                                                                                              >  1:643259/1‑150 (MQ=255)
                                                          tCATGTTGGGGGTACTTACCGGCGATCCGACAGGTTACGGGGCGGAGACCTCGCGGGATATCGCTATTTATGAAAATATTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTAAATGTATTTATTTAAAATACCCCCTGAaaaga                                                                                             <  1:711390/150‑1 (MQ=14)
                                                              tttGTGGGTCCTTCCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAgg                                                                                         <  2:127396/150‑1 (MQ=255)
                                                                tgtgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAAAGTATTTATTTAAAATGCACCCTGAAAAGAAAGGaa                                                                                       >  1:1136357/1‑150 (MQ=255)
                                                                 gtgGGTCCTTTCCGGCGATCCGACAGGGTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGaaa                                                                                      <  1:87203/150‑1 (MQ=255)
                                                                        cTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAgg                                                                               <  1:746017/150‑1 (MQ=39)
                                                                          ttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTg                                                                             >  2:665043/1‑150 (MQ=39)
                                                                           tCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGc                                                                            >  1:390888/1‑150 (MQ=39)
                                                                           tCCGGCGATCCGAAAGGTTACGGGGCGGCTACCTCGCGGGTTTTCGATATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGc                                                                            <  1:116351/150‑1 (MQ=255)
                                                                               gCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGaa                                                                        >  2:679293/1‑150 (MQ=38)
                                                                                     cGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGAGAGGTGCTGAAAGCGAg                                                                  >  1:819798/1‑150 (MQ=255)
                                                                                     cGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCGATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACGTAATGTATTTATTTAAAATACCCCCTCAAAAGACAGGCAACCACCGGTGCTGAACGCGAg                                                                  >  2:625674/1‑150 (MQ=14)
                                                                                      gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                 >  2:444766/1‑150 (MQ=39)
                                                                                           gTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttttt                                                            >  1:870580/1‑150 (MQ=40)
                                                                                              aCGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGc                                                         <  2:643259/150‑1 (MQ=255)
                                                                                               cGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGcc                                                        >  1:127689/1‑150 (MQ=255)
                                                                                                 ggccggcgACCTCGCGGGATTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCATCGTCGAAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCctc                                                      <  2:698890/147‑1 (MQ=255)
                                                                                                  ggcggcGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCctct                                                     >  2:1262035/1‑150 (MQ=255)
                                                                                                     ggcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTc                                                  >  1:257182/1‑150 (MQ=255)
                                                                                                      gcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGACAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCg                                                 <  2:844704/150‑1 (MQ=255)
                                                                                                       ctaCCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGGGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGt                                                <  2:1071026/148‑1 (MQ=255)
                                                                                                        gACCTCGCGGGTTCTCGCTATTTCTGAAAATTTTCCGGTTTAAGGTGTTTCCGTGCTTCTTCGTCGTAACTTACTGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGtt                                               <  1:919701/150‑1 (MQ=255)
                                                                                                         aCCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGttt                                              >  2:514404/1‑150 (MQ=255)
                                                                                                          ccTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTc                                             >  2:621872/1‑150 (MQ=255)
                                                                                                           cTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTcc                                            <  1:242991/150‑1 (MQ=255)
                                                                                                            tCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCt                                           >  1:666447/1‑150 (MQ=255)
                                                                                                             cgcgGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGATACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCtt                                          >  2:1029848/1‑150 (MQ=255)
                                                                                                             cgcgGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCATGATAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTGATTGGCCTATGTCGTTTccgt                                          >  2:1248295/1‑148 (MQ=14)
                                                                                                                gggTTTTCTCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCTTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                       <  1:417925/150‑1 (MQ=255)
                                                                                                                gggTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                       >  1:41931/1‑150 (MQ=255)
                                                                                                                 tgTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCACCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGATCTTTTTGGCCTCTGTCGTTTCCTTtctc                                      >  1:209977/2‑150 (MQ=255)
                                                                                                                          tATTTATGAAAATTTTCCGGCTTAAGGTGTTTCCGTGATTCTTCGTCTTAACTTAATGTATTTATTTAAAAAACCCCCTGAAACGCAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGt                             <  2:193221/150‑1 (MQ=14)
                                                                                                                           aTTTATGACACATTTCCGGTTTACGGCGTTTCCGGTCTTCTTCGGCGTAACTTATTGTATTTATTTATAGTACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTc                            <  1:96818/150‑1 (MQ=11)
                                                                                                                           aTTTATGAAAATTTTCCGGTTTAAGGCGTTTCCGTTCTTGTTCGTCGTAACTTAATGTATTTATTTAAAGTACCCCCTGAAAAGAAAGGAAACGTCAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGtg                            <  1:96816/150‑2 (MQ=255)
                                                                                                                            tttATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTcc                           >  2:78048/1‑150 (MQ=255)
                                                                                                                              tATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACAACCTGAAAAGAAAGGAAACGACTGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGt                         >  1:203944/1‑150 (MQ=255)
                                                                                                                                tGAAATTTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTTGTAACTTAATGTATTTATTTAAATTACCGCCTGAAAAGATAGGAAACGACAGGTGCTGATAGGGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTgg                       <  1:202121/150‑1 (MQ=14)
                                                                                                                                           cGGTTTAAGGTGTTTCCGTTCTTCTTCGGCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATg            <  1:98911/150‑1 (MQ=255)
                                                                                                                                             gTTTAAGGTGTTTCCTTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATTCCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGa          <  1:810963/150‑1 (MQ=255)
                                                                                                                                               ttAAGGTTTTTCAGTTCTTGTTCTTCTTACATTATTTTTTTTATTTAAAATACACCGTGAAACCCAAGTACACGACAGGTGCTGAATGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAg        <  1:1263270/150‑1 (MQ=14)
                                                                                                                                                tAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGt       >  1:587870/1‑150 (MQ=255)
                                                                                                                                                     tgtTTCCGTTATTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaaca  <  2:308063/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATTTGACGTGGTTTGATGGCGTAGATGCACGTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACA  >  NC_000913/580548‑580846

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: