Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 774,648 A→T 5.4% intergenic (+46/‑104) ybgE → / → ybgC putative inner membrane protein in cydABX‑ybgE operon/acyl‑CoA thioester hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913774,6480AT 5.4% 314.1 / 4.3 107intergenic (+46/‑104)ybgE/ybgCputative inner membrane protein in cydABX‑ybgE operon/acyl‑CoA thioester hydrolase
Reads supporting (aligned to +/- strand):  ref base A (51/49);  new base T (4/2);  total (56/51)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

Only 100 of 108 total aligned reads displayed.

CGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAG  >  NC_000913/774500‑774797
                                                                                                                                                    |                                                                                                                                                     
  gCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAAttt                                                                                                                                                    <  1:1432213/150‑1 (MQ=255)
  gCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTACAGGCCCATAAAttt                                                                                                                                                    >  2:361364/1‑150 (MQ=255)
   cGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTTTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAAtttt                                                                                                                                                   >  2:1578974/1‑150 (MQ=255)
     tGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTTACCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTa                                                                                                                                                 <  2:1168008/150‑1 (MQ=255)
     tGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCTGGATGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTa                                                                                                                                                 <  2:41776/150‑1 (MQ=255)
     gtggCTTTCGACCGCATAAGGTGCTATGGCAAGGGATTTTTTGCCCATTGCTTGCCGAAATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTCAATCAGAATTCTCTTTAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTa                                                                                                                                                 <  2:258098/148‑1 (MQ=25)
       ggCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATTAATTTTTACg                                                                                                                                               >  1:124471/1‑150 (MQ=255)
         cctgcggccGCAGAAGGTGCCCTGGCACGGGACTTATTGCCCACTGCTTGCCGATATTGATCTCAATCTCGGGCTGATTTTCATCTTCTTTTACATCAGAATTCTCTTTAAAAACTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCt                                                                                                                                             <  2:1286422/143‑1 (MQ=255)
             cGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCt                                                                                                                                         >  2:609961/1‑150 (MQ=255)
              gTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCtt                                                                                                                                        <  2:1213525/150‑1 (MQ=255)
                ccGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTaa                                                                                                                                      >  1:1564101/1‑150 (MQ=255)
                  gCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACt                                                                                                                                    >  1:1658951/1‑150 (MQ=255)
                  gCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACt                                                                                                                                    <  1:908279/150‑1 (MQ=255)
                     gAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGc                                                                                                                                 <  2:617754/150‑1 (MQ=255)
                       aGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGccc                                                                                                                               >  1:1672647/1‑150 (MQ=255)
                        ggTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCt                                                                                                                              <  2:1592523/150‑1 (MQ=255)
                            cTTTGGCCAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGCGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCAtt                                                                                                                          <  2:51953/150‑1 (MQ=255)
                             tttGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTc                                                                                                                         <  1:981797/150‑1 (MQ=255)
                              ttGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTcc                                                                                                                        <  2:1677173/150‑1 (MQ=255)
                               tGGCATGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTAATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTATAAAAAATTATGGGCCTCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTccc                                                                                                                       <  2:1452588/150‑1 (MQ=255)
                               tGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTccc                                                                                                                       <  2:661225/150‑1 (MQ=255)
                                ggCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCa                                                                                                                      >  1:1563261/1‑150 (MQ=255)
                                    aGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCTGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTACCTTAACTTGCCCTCATTCCCAAAca                                                                                                                  >  1:26016/1‑149 (MQ=255)
                                    aGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCGCAAAcc                                                                                                                  <  2:1138518/150‑1 (MQ=255)
                                      ggATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTc                                                                                                                <  1:480414/150‑1 (MQ=255)
                                      ggATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTACCTTATCTTGCCCTCATTACCAAACCTc                                                                                                                >  1:49331/1‑150 (MQ=255)
                                       gATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGATCAAGGCCCATAAATTTTTACGCTCAATTAACTTGCCATCATTACCAAACCTCa                                                                                                               >  1:1041528/1‑150 (MQ=255)
                                         ttttttGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAAt                                                                                                             <  2:124471/150‑1 (MQ=255)
                                          tttttGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATc                                                                                                            >  1:1567842/1‑150 (MQ=255)
                                               gCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATcgcgcg                                                                                                       <  2:744547/150‑1 (MQ=255)
                                                cccATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATcgcgcgc                                                                                                      >  2:1500610/1‑150 (MQ=255)
                                                cccATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATcgcgcgc                                                                                                      <  1:472297/150‑1 (MQ=255)
                                                cccATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATcgcgcgc                                                                                                      >  2:428989/1‑150 (MQ=255)
                                                     tgcttgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGtata                                                                                                 <  1:615517/150‑1 (MQ=255)
                                                         tgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTAtagtag                                                                                             >  2:1329908/1‑150 (MQ=255)
                                                         tgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTAtagtag                                                                                             >  2:1059527/1‑150 (MQ=255)
                                                          gcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTAACGCTCCCTTAACTAGCCATCATTCCCAAACATCAATCGCGCGCGTAtattaga                                                                                            >  2:817363/1‑149 (MQ=255)
                                                            cGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCAAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTagcag                                                                                          >  1:1175398/1‑150 (MQ=255)
                                                             gATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTagcagc                                                                                         <  1:1574091/150‑1 (MQ=255)
                                                               tatTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGt                                                                                       >  1:412592/1‑150 (MQ=255)
                                                               tatTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGt                                                                                       >  2:1712400/1‑150 (MQ=255)
                                                               tatTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGACCTCATTCCCAAACATCAATCGCGAGCGTATAGTAGCAGCGt                                                                                       >  1:1192036/1‑150 (MQ=255)
                                                                  tGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTa                                                                                    >  2:1195944/1‑150 (MQ=255)
                                                                   gTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaa                                                                                   >  2:24412/1‑150 (MQ=255)
                                                                   gTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaa                                                                                   >  1:1585297/1‑150 (MQ=255)
                                                                   gTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaa                                                                                   <  1:417852/150‑1 (MQ=255)
                                                                   gTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaa                                                                                   <  1:1580601/150‑1 (MQ=255)
                                                                   gTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaa                                                                                   <  1:1510127/150‑1 (MQ=255)
                                                                     tctcTTTGTCGGGCTGATTTTCTTCTTCTTTTTAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTATCTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaaaa                                                                                 <  1:844091/150‑1 (MQ=255)
                                                                      ctcATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAg                                                                                >  2:1481971/1‑150 (MQ=255)
                                                                      ctcATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCACAAACCTCAAACCCGCGCTTATAGTAGAATCGTTTAAAAg                                                                                >  1:535526/1‑150 (MQ=255)
                                                                       tcATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGACCTAATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGc                                                                               >  1:54506/1‑150 (MQ=255)
                                                                         aTTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGctc                                                                             <  1:282941/150‑1 (MQ=255)
                                                                          ttGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGctct                                                                            <  2:1018482/150‑1 (MQ=255)
                                                                             tCGGGCTTATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCTGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAAc                                                                         <  2:100784/150‑1 (MQ=255)
                                                                                ggCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCATGCCCATAAATTTTTACGATACCTTAAATTGCCCTCATTCACAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTATAACtat                                                                      >  1:1191236/1‑148 (MQ=255)
                                                                                ggCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACttt                                                                      <  1:609961/150‑1 (MQ=255)
                                                                                ggCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACttt                                                                      >  1:1363819/1‑150 (MQ=255)
                                                                                 gCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACtttt                                                                     <  2:461689/150‑1 (MQ=255)
                                                                                 gCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCCCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACtttt                                                                     >  2:1273287/1‑150 (MQ=255)
                                                                                  gtGATTTTCTTCTTCTATTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGACCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCATACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTttg                                                                    <  1:160026/149‑1 (MQ=255)
                                                                                         tcttcttctTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCagt                                                             >  1:155289/1‑148 (MQ=255)
                                                                                              ttcttTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCgg                                                        >  2:493954/1‑150 (MQ=255)
                                                                                               tattttAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTAATCTCCCTTAACTTGCCCTCATTCCCAAAACTCAATAGAGCGCGTATAGTAGCAGAGTTTAAAAGATCTAACTATTGTTGCATTACAggg                                                       >  2:535749/3‑150 (MQ=255)
                                                                                                cttTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGa                                                      <  1:361364/150‑1 (MQ=255)
                                                                                                 ttttAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGAt                                                     <  1:1712400/150‑1 (MQ=255)
                                                                                                      aaTCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTACCAAACCTCAATCGCGCGCGTATAGGAGCAGCGTTTAAACGCTCTAACTTTTGTTGCATTACCGGGATGTaga                                                >  1:293566/1‑148 (MQ=255)
                                                                                                       aTCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAACTTTTTACGCGCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATGACCGGGATGTAAAg                                               >  1:849827/1‑150 (MQ=255)
                                                                                                       aTCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAg                                               <  1:1578974/150‑1 (MQ=255)
                                                                                                       aTCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAg                                               >  1:283230/1‑150 (MQ=255)
                                                                                                        tCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGt                                              >  2:383456/1‑150 (MQ=255)
                                                                                                        tCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGt                                              >  2:747928/1‑150 (MQ=255)
                                                                                                        acataaTTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGt                                              <  2:1072062/146‑1 (MQ=255)
                                                                                                         cAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACATACCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTg                                             >  2:18663/1‑150 (MQ=255)
                                                                                                           gAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCCCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGaa                                           <  1:550918/150‑1 (MQ=255)
                                                                                                             aTTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTGGCATTACCGGGATGTAAAGTGAATa                                         <  1:952659/150‑1 (MQ=255)
                                                                                                              ttCTCTTTAAATAATTTTTGTAGTCTCCAGTCCTTTAAATTTTTACGCTTCCTTAAATTGCCCTCATTCCCAAATCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAATGCACTAACTTTTGTTTCATTACATGTATGTAATGTGCATAc                                        <  2:619448/150‑1 (MQ=255)
                                                                                                                 tctTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAAc                                     <  1:747928/150‑1 (MQ=255)
                                                                                                                 tctTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAAc                                     <  1:383456/150‑1 (MQ=255)
                                                                                                                  ctTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACg                                    >  1:1186235/1‑150 (MQ=255)
                                                                                                                    gtAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCt                                  >  2:1409922/2‑150 (MQ=38)
                                                                                                                    gtAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCt                                  >  2:379315/2‑150 (MQ=38)
                                                                                                                    gtAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATCTAAAGTGAAaacaaccca                                  >  2:1191978/2‑147 (MQ=37)
                                                                                                                    gtAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTACGCTCCCTTAACTTGCCCACATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAAATTTTGTTTCATTACCGGGATGTAAAGTGAATACAACGCt                                  >  2:331369/2‑150 (MQ=37)
                                                                                                                      aaaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGt                                >  1:1500848/1‑150 (MQ=255)
                                                                                                                       aaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTATAGTGAATACAACGCTGtt                               >  2:200163/1‑150 (MQ=255)
                                                                                                                       aaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGtt                               <  1:1329908/150‑1 (MQ=255)
                                                                                                                       aaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCACAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTATAACTTTTGTTGCATTTCCGGGATGTAATGTGAATACATCGCTGtt                               >  2:188566/1‑150 (MQ=255)
                                                                                                                        aaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGttt                              <  2:1658951/150‑1 (MQ=255)
                                                                                                                            ttATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGAt                          <  2:663828/150‑1 (MQ=255)
                                                                                                                               tGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCTTTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGc                       >  1:139412/1‑150 (MQ=255)
                                                                                                                                 ggCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCg                     >  1:1376771/1‑150 (MQ=255)
                                                                                                                                  gCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAACCCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGAGGTAAAGTGAATACAACGCTGTTTCGATGGCCgg                    >  1:1729895/1‑150 (MQ=255)
                                                                                                                                      cTCCAGGCAGATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTTTTGCATTACCGGGATGTTAAGTGAATACTACGCTGTTTCGATGGCCGGTTcg                <  1:1151541/150‑1 (MQ=255)
                                                                                                                                          aGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTc            <  2:1500848/150‑1 (MQ=255)
                                                                                                                                           ggCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTACCTAACCTGAATCTCGTTCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTATGTTGCATTACCGGGATGTAAAGTGAATACAACGCGGTTTCGATGGCCGGTTCGCGTct           <  1:1077130/150‑1 (MQ=255)
                                                                                                                                            gCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAACGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTcta          >  1:1041978/1‑150 (MQ=255)
                                                                                                                                            gCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTcta          >  2:1021983/1‑150 (MQ=255)
                                                                                                                                                aTAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTAt      >  2:1574049/1‑150 (MQ=255)
                                                                                                                                                    aTTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGGATAGTAGCATCGTTTAAAAGCTCTAAGTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAg  <  1:345527/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAG  >  NC_000913/774500‑774797

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: