Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,078,306:1 +G 5.4% intergenic (+172/‑269) 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,3061.G 5.4% 464.7 / 21.1 112intergenic (+172/‑269)yeeW/yoeFCP4‑44 prophage; putative uncharacterized protein YeeW/CP4‑44 prophage; putative uncharacterized protein YoeF
Reads supporting (aligned to +/- strand):  ref base . (57/49);  new base G (2/4);  total (59/53)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

Only 100 of 113 total aligned reads displayed.

TAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACC  >  NC_000913/2078157‑2078446
                                                                                                                                                      |                                                                                                                                            
taGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACa                                                                                                                                               >  2:878442/1‑150 (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)
             cacaGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCa                                                                                                                                 <  1:1285239/150‑1 (MQ=255)
             cacaGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCa                                                                                                                                 <  1:1312229/150‑1 (MQ=255)
                aGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCACCCTGAATACAGCTGGCGTTTTCAttt                                                                                                                              >  1:803236/1‑150 (MQ=255)
                 gAAACAGCACTGTTGTCCTGTTTTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTta                                                                                                                             <  1:1049196/150‑1 (MQ=255)
                 gAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTta                                                                                                                             >  1:43026/1‑150 (MQ=255)
                  aaaCAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTtat                                                                                                                            <  1:506256/150‑1 (MQ=255)
                  aaaCAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTtat                                                                                                                            <  2:908617/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)
                            gttgtCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCATGCGCCATCTACAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAg                                                                                                                  >  2:1160900/1‑150 (MQ=255)
                                   tGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAAcc                                                                                                           <  1:1381559/150‑1 (MQ=255)
                                     tATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGa                                                                                                         >  1:1390645/1‑150 (MQ=255)
                                      aTTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGat                                                                                                        <  2:1196785/150‑1 (MQ=255)
                                      aTTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGat                                                                                                        <  1:539120/150‑1 (MQ=255)
                                       ttAAGCAGGCAAGAGAAGAGGTCCTTACCCTGCGCCAGCTGCAGCGTCTGCACCAGAACCGCGGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCGGGCGGTTTCATTTATAAACAGAAAGGACAACCGata                                                                                                       <  2:864110/150‑1 (MQ=255)
                                        tAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataa                                                                                                      <  2:1390645/150‑1 (MQ=255)
                                        tAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataa                                                                                                      >  2:79668/1‑150 (MQ=255)
                                        tAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataa                                                                                                      >  2:441518/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:967098/1‑150 (MQ=255)
                                                            tCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCAATTATAAACAGAAATGAAAACCGATAATAATGGAATGGACGCCCCGa                                                                                  >  2:1292661/1‑150 (MQ=255)
                                                                  cccTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGa                                                                            >  2:47956/1‑150 (MQ=255)
                                                                    cTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAg                                                                          >  1:1163803/1‑150 (MQ=255)
                                                                    cTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAg                                                                          <  1:773664/150‑1 (MQ=255)
                                                                     tGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGc                                                                         <  2:16448/150‑1 (MQ=255)
                                                                         ccATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAACGGAATGGACGCCCCGACCATGAAGCAACg                                                                     >  1:574498/1‑150 (MQ=255)
                                                                          cATCTGCAGCTTCTGCACCAGACCCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGAAAATAATGGAATGGACGCCCCGACCATGAAGCAACgg                                                                    <  2:1365943/150‑1 (MQ=255)
                                                                          cATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACgg                                                                    <  1:617002/150‑1 (MQ=255)
                                                                              tGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAgg                                                                >  1:1546918/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)
                                                                                   cTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATAtt                                                           <  1:79668/150‑1 (MQ=255)
                                                                                       tGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTg                                                       <  2:1216763/150‑1 (MQ=255)
                                                                                        gCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGa                                                      >  1:900383/1‑150 (MQ=255)
                                                                                        gCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGa                                                      >  1:1382456/1‑150 (MQ=255)
                                                                                          aCCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATa                                                    <  1:1495728/150‑1 (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)
                                                                                                    gTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAg                                          <  1:895399/150‑1 (MQ=255)
                                                                                                        aTCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGa                                      <  1:480157/150‑1 (MQ=255)
                                                                                                        aTCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGa                                      >  2:818906/1‑150 (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)
                                                                                                                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)
                                                                                                                 aaCAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTa                             >  2:149481/1‑150 (MQ=255)
                                                                                                                     ggCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGc                         >  1:1386429/1‑150 (MQ=255)
                                                                                                                     ggCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGc                         >  1:753926/1‑150 (MQ=255)
                                                                                                                      gCTTCATTTCTCCTGACGCCAGCCTGAATACAGCGGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGcc                        <  2:699657/150‑1 (MQ=255)
                                                                                                                       cTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAAACGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCg                       >  2:1227359/1‑150 (MQ=255)
                                                                                                                        ttCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGt                      >  1:1129620/1‑150 (MQ=255)
                                                                                                                          cATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGttt                    <  2:309231/150‑1 (MQ=255)
                                                                                                                          cATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGttt                    <  2:682321/150‑1 (MQ=255)
                                                                                                                          cATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGttt                    >  1:1442382/1‑150 (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)
                                                                                                                               ctcCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACg               <  1:875999/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)
                                                                                                                                            ccTGAATACAGCTGGCGTTTTCATTAATAAACACAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGAcc  <  1:1165812/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                            
TAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACC  >  NC_000913/2078157‑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: