Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 774,637 C→G 5.4% intergenic (+35/‑115) 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,6370CG 5.4% 343.6 / 9.4 111intergenic (+35/‑115)ybgE/ybgCputative inner membrane protein in cydABX‑ybgE operon/acyl‑CoA thioester hydrolase
Reads supporting (aligned to +/- strand):  ref base C (51/52);  new base G (4/2);  total (57/54)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-01

Only 100 of 112 total aligned reads displayed.

TGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCG  >  NC_000913/774488‑774783
                                                                                                                                                     |                                                                                                                                                  
gggTGTGATTCACGGCGGGGGCTGTCGTCGGCAGAAGGTTCTTTGGCCAGGGACTTTGTGCCCATTGCTTGCCGATCTTGTTCTCCATGTCGGGCTGATTTTCTTCTTCTTTTAACTCAGAATTCTCTTTAAAAAATTATGGGCCGCTcc                                                                                                                                                    <  1:266977/149‑1 (MQ=255)
   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)
        ttCACGGCGTGGGATTTCGTACGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTGTCATTGTCGGGCTGATTTTCTTCTTCTTTTATATCAGAATTCTCTTTAAAAAATTATGTGCCGATCCAGGaccct                                                                                                                                            >  2:121103/1‑148 (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)
                                         tttGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTc                                                                                                           <  1:981797/150‑1 (MQ=255)
                                          ttGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTcc                                                                                                          <  2:1677173/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                                                                                        <  1:472297/150‑1 (MQ=255)
                                                            cccATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATcgcgcgc                                                                                        >  2:1500610/1‑150 (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                                                                         >  2:1712400/1‑150 (MQ=255)
                                                                           tatTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGACCTCATTCCCAAACATCAATCGCGAGCGTATAGTAGCAGCGt                                                                         >  1:1192036/1‑150 (MQ=255)
                                                                               gTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaa                                                                     <  1:417852/150‑1 (MQ=255)
                                                                               gTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaa                                                                     <  1:1510127/150‑1 (MQ=255)
                                                                               gTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaa                                                                     >  1:1585297/1‑150 (MQ=255)
                                                                               gTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaa                                                                     <  1:1580601/150‑1 (MQ=255)
                                                                               gTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaa                                                                     >  2:24412/1‑150 (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)
                                                                                        gTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGTGCCGCTCCAGGCCCATAAATATTTACGCTCACTTAACTTGCAATCATTCCCAAACCTAAATCGCGAGAGTATAGTAGAATCGTTTAAAAGCTCTaa                                                            >  2:444815/1‑150 (MQ=255)
                                                                                         tCGGGCTTATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCTGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAAc                                                           <  2:100784/150‑1 (MQ=255)
                                                                                            ggCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCATGCCCATAAATTTTTACGATACCTTAAATTGCCCTCATTCACAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTATAACtat                                                        >  1:1191236/1‑148 (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:283230/1‑150 (MQ=255)
                                                                                                                   aTCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAg                                 <  1:1578974/150‑1 (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)
                                                                                                                          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: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)
                                                                                                                                  aaaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGt                  >  2:1263873/1‑150 (MQ=255)
                                                                                                                                   aaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTATAGTGAATACAACGCTGtt                 >  2:200163/1‑150 (MQ=255)
                                                                                                                                   aaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGtt                 <  1:1329908/150‑1 (MQ=255)
                                                                                                                                   aaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGtt                 <  1:428989/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)
                                                                                                                                                     |                                                                                                                                                  
TGGTGTGATTCACGGCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCG  >  NC_000913/774488‑774783

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: