Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 774,638 A→C 5.0% intergenic (+36/‑114) 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,6380AC 5.0% 326.3 / 4.8 111intergenic (+36/‑114)ybgE/ybgCputative inner membrane protein in cydABX‑ybgE operon/acyl‑CoA thioester hydrolase
Reads supporting (aligned to +/- strand):  ref base A (53/51);  new base C (4/2);  total (57/54)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 112 total aligned reads displayed.

TGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTC  >  NC_000913/774491‑774787
                                                                                                                                                   |                                                                                                                                                     
tgtgATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAgg                                                                                                                                                     >  2:282941/1‑150 (MQ=255)
 gtgATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGc                                                                                                                                                    >  1:617754/1‑150 (MQ=255)
  tgAGTCACTGCGGGGTCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGATCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGcc                                                                                                                                                   <  1:1307667/150‑1 (MQ=255)
   gATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGccc                                                                                                                                                  <  2:4483/150‑1 (MQ=255)
   gATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGAGTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCCGGCGGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGccc                                                                                                                                                  <  2:1160931/150‑1 (MQ=255)
   gATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGAATTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGccc                                                                                                                                                  >  1:661225/1‑150 (MQ=255)
    aTTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGCATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCa                                                                                                                                                 <  2:894944/150‑1 (MQ=255)
         cGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAAt                                                                                                                                            <  1:1658800/150‑1 (MQ=255)
           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)
              gtggCTTTCGACCGCATAAGGTGCTATGGCAAGGGATTTTTTGCCCATTGCTTGCCGAAATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTCAATCAGAATTCTCTTTAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTa                                                                                                                                       <  2:258098/148‑1 (MQ=25)
                ggCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTaaa                                                                                                                                     <  1:1405954/150‑3 (MQ=38)
                ggCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATTAATTTTTACg                                                                                                                                     >  1:124471/1‑150 (MQ=255)
                  cctgcggccGCAGAAGGTGCCCTGGCACGGGACTTATTGCCCACTGCTTGCCGATATTGATCTCAATCTCGGGCTGATTTTCATCTTCTTTTACATCAGAATTCTCTTTAAAAACTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCt                                                                                                                                   <  2:1286422/143‑1 (MQ=255)
                    ttCGGCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTcc                                                                                                                                 <  1:937265/150‑1 (MQ=255)
                      cGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCt                                                                                                                               >  2:609961/1‑150 (MQ=255)
                       gTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCtt                                                                                                                              <  2:1213525/150‑1 (MQ=255)
                           gCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACt                                                                                                                          <  1:908279/150‑1 (MQ=255)
                           gCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACt                                                                                                                          >  1:1658951/1‑150 (MQ=255)
                              gAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGc                                                                                                                       <  2:617754/150‑1 (MQ=255)
                                 ggTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCt                                                                                                                    <  2:1592523/150‑1 (MQ=255)
                                     cTTTGGCCAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGCGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCAtt                                                                                                                <  2:51953/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)
                                                                  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:1580601/150‑1 (MQ=255)
                                                                            gTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaa                                                                         <  1:1510127/150‑1 (MQ=255)
                                                                              tctcTTTGTCGGGCTGATTTTCTTCTTCTTTTTAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTATCTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaaaa                                                                       <  1:844091/150‑1 (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)
                                                                                     gTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGTGCCGCTCCAGGCCCATAAATATTTACGCTCACTTAACTTGCAATCATTCCCAAACCTAAATCGCGAGAGTATAGTAGAATCGTTTAAAAGCTCTaa                                                                >  2:444815/1‑150 (MQ=255)
                                                                                      tCGGGCTTATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCTGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAAc                                                               <  2:100784/150‑1 (MQ=255)
                                                                                         ggCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCATGCCCATAAATTTTTACGATACCTTAAATTGCCCTCATTCACAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTATAACtat                                                            >  1:1191236/1‑148 (MQ=255)
                                                                                         ggCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACttt                                                            <  1:609961/150‑1 (MQ=255)
                                                                                          gCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACtttt                                                           <  2:461689/150‑1 (MQ=255)
                                                                                          gCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCCCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACtttt                                                           >  2:1273287/1‑150 (MQ=255)
                                                                                           gtGATTTTCTTCTTCTATTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGACCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCATACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTttg                                                          <  1:160026/149‑1 (MQ=255)
                                                                                             gATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTttgtt                                                        >  2:1661966/1‑150 (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:747928/1‑150 (MQ=255)
                                                                                                                 tCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGt                                    >  2:383456/1‑150 (MQ=255)
                                                                                                                 acataaTTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGt                                    <  2:1072062/146‑1 (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:383456/150‑1 (MQ=255)
                                                                                                                          tctTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAAc                           <  1:747928/150‑1 (MQ=255)
                                                                                                                           ctTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACg                          >  1:1186235/1‑150 (MQ=255)
                                                                                                                             gtAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCt                        >  2:379315/2‑150 (MQ=38)
                                                                                                                             gtAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCt                        >  2:1409922/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)
                                                                                                                               aaaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGt                      >  2:1263873/1‑150 (MQ=255)
                                                                                                                                aaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTATAGTGAATACAACGCTGtt                     >  2:200163/1‑150 (MQ=255)
                                                                                                                                aaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGtt                     <  1:428989/150‑1 (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)
                                                                                                                                                   |                                                                                                                                                     
TGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTC  >  NC_000913/774491‑774787

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: