Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,078,297 Δ1 bp 5.3% intergenic (+163/‑278) yeeW → / → yoeF CP4‑44 prophage; putative uncharacterized protein YeeW/CP4‑44 prophage; putative uncharacterized protein YoeF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,078,2970C. 5.3% 382.4 / 13.3 114intergenic (+163/‑278)yeeW/yoeFCP4‑44 prophage; putative uncharacterized protein YeeW/CP4‑44 prophage; putative uncharacterized protein YoeF
Reads supporting (aligned to +/- strand):  ref base C (58/48);  new base . (2/4);  total (61/53)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

Only 100 of 116 total aligned reads displayed.

AGGGTATAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACC  >  NC_000913/2078151‑2078446
                                                                                                                                                  |                                                                                                                                                      
aGGGTATAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCATATTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTg                                                                                                                                                     >  1:555901/1‑150 (MQ=255)
aGGGTATAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTg                                                                                                                                                     <  1:567027/150‑1 (MQ=255)
        gTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGGTGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGc                                                                                                                                            >  1:50474/1‑150 (MQ=255)
         tGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCt                                                                                                                                           >  1:111302/1‑150 (MQ=255)
            cATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGc                                                                                                                                        >  2:946322/1‑150 (MQ=255)
            cATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGc                                                                                                                                        >  2:940555/1‑150 (MQ=255)
               tAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCTGTATTCAGGCTGGCGtc                                                                                                                                     <  1:558954/150‑2 (MQ=255)
               tAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCTGTATTCAGGCTGGCGtc                                                                                                                                     <  1:567917/150‑2 (MQ=255)
               tAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCTGTATTCAGGCTGGCGtc                                                                                                                                     <  2:1572588/150‑2 (MQ=255)
                aaCCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGttt                                                                                                                                    >  2:1381559/1‑150 (MQ=255)
                  ccACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCTGTATTCAGGCTGGCGtcagg                                                                                                                                  <  2:873833/150‑5 (MQ=255)
                   cacaGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCa                                                                                                                                 <  1:1285239/150‑1 (MQ=255)
                   cacaGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCa                                                                                                                                 <  1:1312229/150‑1 (MQ=255)
                      aGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCACCCTGAATACAGCTGGCGTTTTCAttt                                                                                                                              >  1:803236/1‑150 (MQ=255)
                      aGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCCCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCAttt                                                                                                                              <  1:841543/150‑1 (MQ=255)
                       gAAACAGCACTGTTGTCCTGTTTTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTta                                                                                                                             <  1:1049196/150‑1 (MQ=255)
                        aaaCAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTtat                                                                                                                            <  2:908617/150‑1 (MQ=255)
                        aaaCAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTtat                                                                                                                            <  1:506256/150‑1 (MQ=255)
                             gCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACa                                                                                                                       >  2:1550818/1‑150 (MQ=255)
                             gCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACa                                                                                                                       >  1:312647/1‑150 (MQ=255)
                             gCACTGTTGTCCTGTATTAAGCACGCAAGAGAAGATGACCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGTACCGGTGAATCTCACTCAACCCGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACa                                                                                                                       <  2:1356505/150‑1 (MQ=255)
                              cACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAg                                                                                                                      >  1:1248121/1‑150 (MQ=255)
                                         tGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAAcc                                                                                                           <  1:1381559/150‑1 (MQ=255)
                                         tGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAAcc                                                                                                           >  1:1196785/1‑150 (MQ=255)
                                           tATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGa                                                                                                         >  1:1390645/1‑150 (MQ=255)
                                            aTTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGat                                                                                                        <  2:1196785/150‑1 (MQ=255)
                                             ttAAGCAGGCAAGAGAAGAGGTCCTTACCCTGCGCCAGCTGCAGCGTCTGCACCAGAACCGCGGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCGGGCGGTTTCATTTATAAACAGAAAGGACAACCGata                                                                                                       <  2:864110/150‑1 (MQ=255)
                                              tAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataa                                                                                                      >  2:441518/1‑150 (MQ=255)
                                              tAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataa                                                                                                      <  2:1390645/150‑1 (MQ=255)
                                              tAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataa                                                                                                      >  2:79668/1‑150 (MQ=255)
                                               aaGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataat                                                                                                     >  2:763627/1‑150 (MQ=255)
                                                   aGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatg                                                                                                 <  2:1624685/150‑1 (MQ=255)
                                                   aGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatg                                                                                                 <  2:312647/150‑1 (MQ=255)
                                                    ggCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGATAGGAAAACCGATAATaatgg                                                                                                >  2:1165812/1‑150 (MQ=255)
                                                    ggCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatgg                                                                                                <  1:441518/150‑1 (MQ=255)
                                                    ggCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatgg                                                                                                <  1:813365/150‑1 (MQ=255)
                                                        agagAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatggaatg                                                                                            >  2:1540464/1‑150 (MQ=255)
                                                        agagAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatggaatg                                                                                            >  1:1405328/1‑150 (MQ=255)
                                                        agagAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatggaatg                                                                                            >  1:967098/1‑150 (MQ=255)
                                                                  tCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCAATTATAAACAGAAATGAAAACCGATAATAATGGAATGGACGCCCCGa                                                                                  >  2:1292661/1‑150 (MQ=255)
                                                                          cTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAg                                                                          >  1:1163803/1‑150 (MQ=255)
                                                                               ccATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAACGGAATGGACGCCCCGACCATGAAGCAACg                                                                     >  1:574498/1‑150 (MQ=255)
                                                                                cATCTGCAGCTTCTGCACCAGACCCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGAAAATAATGGAATGGACGCCCCGACCATGAAGCAACgg                                                                    <  2:1365943/150‑1 (MQ=255)
                                                                                cATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACgg                                                                    <  1:617002/150‑1 (MQ=255)
                                                                                 aTCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGt                                                                   <  1:509160/150‑1 (MQ=255)
                                                                                  tCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTa                                                                  >  2:1041408/1‑150 (MQ=255)
                                                                                    tGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAgg                                                                >  2:895399/1‑150 (MQ=255)
                                                                                    tGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAgg                                                                >  2:841760/1‑150 (MQ=255)
                                                                                     gCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGa                                                               >  2:1495728/1‑150 (MQ=255)
                                                                                         cTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATATTGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATAtt                                                           <  1:64464/150‑1 (MQ=255)
                                                                                         cTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATAtt                                                           <  1:79668/150‑1 (MQ=255)
                                                                                             tGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTg                                                       <  2:1216763/150‑1 (MQ=255)
                                                                                              gCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGa                                                      >  1:900383/1‑150 (MQ=255)
                                                                                                aCCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATa                                                    <  1:1495728/150‑1 (MQ=255)
                                                                                                 ccAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATaa                                                   >  1:36401/1‑150 (MQ=255)
                                                                                                 ccAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATaa                                                   >  2:577710/1‑150 (MQ=255)
                                                                                                  cAGAACCGGTGAATCTCACTCAACAGGCTTCATTTATCCTGACGCCAGACTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACTCCCCGACCATGAAGCAACGGTAGGATATTACTGATaaa                                                  <  2:555901/150‑1 (MQ=255)
                                                                                                   aGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAAc                                                 <  1:1265815/150‑1 (MQ=255)
                                                                                                      aCCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCtt                                              >  1:229347/1‑150 (MQ=255)
                                                                                                        cGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTcc                                            <  2:50474/150‑1 (MQ=255)
                                                                                                         ggTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCa                                           >  2:377433/1‑150 (MQ=255)
                                                                                                         ggTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACAGATAATAATGGAATGGACGCCCCGACCATGAAGCAACCGTAGGATATTACTGATAAACCTTCCa                                           >  2:336494/1‑150 (MQ=255)
                                                                                                          gTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAg                                          <  1:895399/150‑1 (MQ=255)
                                                                                                              aTCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGa                                      >  2:818906/1‑150 (MQ=255)
                                                                                                              aTCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGa                                      <  1:480157/150‑1 (MQ=255)
                                                                                                                 tcactcaACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAggg                                   >  1:406022/1‑149 (MQ=255)
                                                                                                                  cactcaACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTg                                  <  2:464683/150‑1 (MQ=255)
                                                                                                                  cactcaACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTg                                  >  2:959357/1‑150 (MQ=255)
                                                                                                                   actcaACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGc                                 <  1:1540464/150‑1 (MQ=255)
                                                                                                                    ctcaACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCg                                <  1:841760/150‑1 (MQ=255)
                                                                                                                      caACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATTAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATTTTACTGATAAACCTTCCAGTGGAGGTGCGtt                              >  2:1502539/1‑150 (MQ=255)
                                                                                                                      caACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGtt                              >  2:1191935/1‑150 (MQ=255)
                                                                                                                      caACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGtt                              >  1:625483/1‑150 (MQ=255)
                                                                                                                      caACAGGCTTCATTTCTCCTGACGCCAGACTGAATACAGATGGCGTTTTCATTTATAAACAGAAAGGAAAAGCGATAATAAGGGAATGGACGCCCGGACCATGAAGCAACGGTAGGATATTACTGATAAACATTCGAGTTGATGTGCggt                              >  2:1030330/1‑148 (MQ=255)
                                                                                                                      caACAGGCTTCATTTCTCCTGACGCAAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCgat                              >  2:1124259/1‑148 (MQ=255)
                                                                                                                       aaCAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCTATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTa                             <  1:1550818/150‑1 (MQ=255)
                                                                                                                           ggCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGc                         >  1:1386429/1‑150 (MQ=255)
                                                                                                                           ggCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGc                         >  1:753926/1‑150 (MQ=255)
                                                                                                                             cTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAAACGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCg                       >  2:1227359/1‑150 (MQ=255)
                                                                                                                              ttCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGt                      >  1:1129620/1‑150 (MQ=255)
                                                                                                                                cATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGttt                    <  2:309231/150‑1 (MQ=255)
                                                                                                                                cATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGttt                    >  1:1442382/1‑150 (MQ=255)
                                                                                                                                cATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGttt                    <  2:682321/150‑1 (MQ=255)
                                                                                                                                cATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAACGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGttt                    <  2:574498/150‑1 (MQ=255)
                                                                                                                                    tctcCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAAc                <  1:47956/150‑1 (MQ=255)
                                                                                                                                    aatgaaaaCGCCAGCTGTATTCAGGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAAc                >  2:1203640/8‑150 (MQ=255)
                                                                                                                                     ctcCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACg               <  1:1041408/150‑1 (MQ=255)
                                                                                                                                      tgaaaaCGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTGTTAACGa              >  1:517274/6‑150 (MQ=255)
                                                                                                                                      tgaaaaCGCCAGCTGTATTCAGGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGa              >  2:1009537/6‑150 (MQ=255)
                                                                                                                                       ccTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGAt             >  2:887893/1‑150 (MQ=255)
                                                                                                                                       ccTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGAt             >  1:1022546/1‑150 (MQ=255)
                                                                                                                                        cTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGAtt            <  1:763627/150‑1 (MQ=255)
                                                                                                                                          gACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAg          <  1:56728/150‑1 (MQ=255)
                                                                                                                                           aCGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGc         >  1:179654/1‑150 (MQ=255)
                                                                                                                                            cGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCa        >  1:723828/1‑150 (MQ=255)
                                                                                                                                             gCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTTCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCAt       >  2:678980/1‑150 (MQ=255)
                                                                                                                                               cAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGTAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTg     <  2:1442382/150‑1 (MQ=255)
                                                                                                                                                 gCCTGGATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGAc   >  1:1005769/1‑150 (MQ=255)
                                                                                                                                                 gCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGAc   <  2:36401/150‑1 (MQ=255)
                                                                                                                                                  ccTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGAcc  <  1:818906/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
AGGGTATAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACC  >  NC_000913/2078151‑2078446

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