Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,275,152 T→A 5.4% intergenic (‑111/‑319) yjcB ← / → yjcC putative inner membrane protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,275,1520TA 5.4% 411.5 / 5.3 137intergenic (‑111/‑319)yjcB/yjcCputative inner membrane protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (53/75);  new base A (2/6);  total (56/81)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

CCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAA  >  NC_000913/4275003‑4275301
                                                                                                                                                     |                                                                                                                                                         
    cGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgt                                                                                                                                                       >  1:1439730/1‑150 (MQ=255)
     gaagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgtt                                                                                                                                                      >  2:453689/1‑150 (MQ=255)
     gaagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgtt                                                                                                                                                      <  2:706607/150‑1 (MQ=255)
     gaagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGGGCCCGGGAGGCGTGGTACCCTCTGGCGATCGCGGGGCGTGCAGAACACAAAAATGCAACACTCtgtt                                                                                                                                                      >  2:87748/1‑150 (MQ=255)
      aagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACCCTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgttt                                                                                                                                                     >  1:17110/1‑150 (MQ=255)
       agaagaaCCACGCCAGTGGTGAGGGTAGCCATCACACAATTCCTCTGAACCGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgtttg                                                                                                                                                    <  2:1272310/150‑1 (MQ=255)
        gaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAGGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgtttgt                                                                                                                                                   >  1:539392/1‑150 (MQ=255)
          agaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGCGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAACCACTCTgtttgttt                                                                                                                                                 >  1:1063933/1‑150 (MQ=255)
           gaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTc                                                                                                                                                <  1:633592/150‑1 (MQ=255)
             aCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGCCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCAt                                                                                                                                              >  2:1388747/1‑150 (MQ=255)
             aCCACGCCAGTGGTGAGGGTAGACATCAACCAATTCCTCTGAAAACAGCCCATGCCCTGGCCACCGCGGCGTTGTTCACGGGAGGCGTGGTACCCTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACGCAGTTTGTTTAAt                                                                                                                                              >  2:570112/1‑150 (MQ=255)
              ccACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTTGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCAtt                                                                                                                                             <  2:254170/150‑1 (MQ=255)
              ccACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACAAACAGAGTGTTTCAtt                                                                                                                                             <  1:364620/150‑1 (MQ=255)
               cACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTa                                                                                                                                            >  2:619306/1‑150 (MQ=255)
                  gCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACAAACAGAGTGTTTCAtttttg                                                                                                                                         <  1:980497/150‑5 (MQ=255)
                  gCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGAGGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACAAACAGAGTGTTTCAtttttg                                                                                                                                         <  1:1377193/150‑5 (MQ=255)
                     aGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTtg                                                                                                                                      >  1:155625/1‑150 (MQ=255)
                       tggggATGGTAGCCATCAAACAATACCTCTGAAAATAGCCGATGCCCGTGACAGCGCGGCGTTGTTCACGGGAGGCGTTGTAAACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTtgtg                                                                                                                                    <  2:354919/150‑1 (MQ=255)
                        ggtgttcgTAGCCATCAAACAATTCCTCTGACAAGAGCCGATGCCCTGGACCGCGCGGCGTTGTGCACGGGAGGCGTGTTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGa                                                                                                                                   <  2:285975/143‑1 (MQ=255)
                        ggtgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCATGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGa                                                                                                                                   >  1:1048915/1‑150 (MQ=255)
                           gAGGGTCGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGTACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACt                                                                                                                                <  1:847730/150‑1 (MQ=255)
                                tAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATc                                                                                                                           <  2:745544/150‑1 (MQ=255)
                                 aGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATca                                                                                                                          >  1:166687/1‑150 (MQ=255)
                                  gCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATcac                                                                                                                         <  2:194160/150‑1 (MQ=255)
                                  cccATCAAACAATACCTCTGAAAAGAGCCGATGCCCTGGACAGCGCTGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATcac                                                                                                                         <  1:1328577/149‑1 (MQ=255)
                                   ccATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATcaca                                                                                                                        <  1:649768/150‑1 (MQ=255)
                                    cATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACa‑‑‑‑a                                                                                                                   >  2:342109/1‑150 (MQ=255)
                                       cAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑attg                                                                                                                >  2:162376/1‑150 (MQ=255)
                                         aaCAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑attgat                                                                                                              <  1:1008514/150‑1 (MQ=255)
                                          aCAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑attgatt                                                                                                             <  1:564758/150‑1 (MQ=255)
                                              ttCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGAttgttt                                                                                                         <  2:1106376/150‑1 (MQ=255)
                                               tCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGAttgtttg                                                                                                        <  1:375744/150‑1 (MQ=255)
                                                ccTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGAttgtttgt                                                                                                       >  2:379171/1‑150 (MQ=255)
                                                ccTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGAGTGGTACACTCTGGCTATCGCGGTGCTTGCAGAACACAAAAATGAAACACTCTGTGTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGAttgtttgt                                                                                                       >  1:444299/1‑150 (MQ=255)
                                                 ctctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATtgtttgtt                                                                                                      <  1:38505/150‑1 (MQ=255)
                                                 ctctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATtgtttgtt                                                                                                      >  2:613471/1‑150 (MQ=255)
                                                 ctctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATtgtttgtt                                                                                                      <  2:587223/150‑1 (MQ=255)
                                                 ctctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATAGATtgtttgtt                                                                                                      >  1:762859/1‑150 (MQ=255)
                                                     gAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAGCACTCTGTTTGTTTGATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGcc                                                                                                  >  1:492063/1‑150 (MQ=255)
                                                     gAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGcc                                                                                                  >  1:530428/1‑150 (MQ=255)
                                                     caacagagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGcc                                                                                                  <  2:1467885/146‑1 (MQ=255)
                                                       aaaGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCat                                                                                                <  1:453689/150‑1 (MQ=255)
                                                          gagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATAtt                                                                                             <  2:917966/150‑1 (MQ=255)
                                                            gCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGCGGCTTGCAGAACACAAAAATGAAACACTCTGCTGGTTTCATAAATATTGTGAACTAAATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAg                                                                                           >  2:333306/1‑150 (MQ=255)
                                                              cGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGc                                                                                         <  1:493568/150‑1 (MQ=255)
                                                                  gCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTg                                                                                     <  2:607226/150‑1 (MQ=255)
                                                                  gCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTg                                                                                     <  1:1061312/150‑1 (MQ=255)
                                                                    ccTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGAc                                                                                   >  2:965279/1‑150 (MQ=255)
                                                                    ccTGCACAGCGCGGCGTTGTTCACTGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGACACACTCTGTTTGTTTCAGTAAGTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGAc                                                                                   <  1:279208/150‑1 (MQ=255)
                                                                     cTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACt                                                                                  >  1:937540/1‑150 (MQ=255)
                                                                     cTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGGTAGCCATATTAGGCCGTGACt                                                                                  <  2:1069728/150‑1 (MQ=255)
                                                                     cTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGATAGCCATATTAGGCCGTGACt                                                                                  <  2:527780/150‑1 (MQ=255)
                                                                          cAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGTGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTAt                                                                             <  1:1596732/150‑1 (MQ=255)
                                                                           aGCGCGCCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTAtt                                                                            <  1:1014348/150‑1 (MQ=255)
                                                                                  cgTTGTTCACGCGAGGCGTGGTACACTCGGGCAATCGGGGGGCTTGAAGAACAGAAAAATGAAACACACTGTTTGTTGCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTAc                                                                     <  2:111820/150‑1 (MQ=255)
                                                                                    ttgttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAg                                                                   <  1:1297826/150‑1 (MQ=255)
                                                                                     tgttCACGGGAGGCGTGGTACACTCTGGCTGTCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGa                                                                  >  1:354836/1‑150 (MQ=255)
                                                                                      gttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGAt                                                                 >  1:1325494/1‑150 (MQ=255)
                                                                                       ttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGAtt                                                                <  2:804899/150‑1 (MQ=255)
                                                                                        tCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACCCTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTa                                                               <  2:1039839/150‑1 (MQ=255)
                                                                                           cGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACT‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtgt                                                            >  2:1592868/1‑150 (MQ=255)
                                                                                            gggAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtgtg                                                           <  2:539392/150‑1 (MQ=255)
                                                                                               aGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGtt                                                        <  2:914241/150‑1 (MQ=255)
                                                                                                   gTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTc                                                    <  1:1411029/150‑1 (MQ=255)
                                                                                                     ggTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGATGTACAGATTATGTGGTTTTTCAg                                                  >  1:1096319/1‑150 (MQ=255)
                                                                                                      gTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGt                                                 >  1:332145/1‑150 (MQ=255)
                                                                                                       tACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTg                                                <  2:1633583/150‑1 (MQ=255)
                                                                                                         cacTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAg                                                  <  2:680464/150‑1 (MQ=255)
                                                                                                         cacTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGc                                              <  2:1554938/150‑1 (MQ=255)
                                                                                                            tctGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGGTTCATTAATTTTGTGAACTATTTCACA‑‑‑‑ATTGATTGTTTGTTAGCCATTTTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCAtt                                           >  2:346659/1‑150 (MQ=255)
                                                                                                               ggCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGCCTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAg                                        >  1:331100/1‑150 (MQ=255)
                                                                                                                  tATCGCGGTGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGc                                     <  1:208774/150‑1 (MQ=255)
                                                                                                                  tATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGc                                     <  1:357788/150‑1 (MQ=255)
                                                                                                                  tATCGCGGGGCTTGCAGAACAATAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGc                                     <  1:119345/150‑1 (MQ=255)
                                                                                                                   aTCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCa                                    <  2:189430/150‑1 (MQ=255)
                                                                                                                     cgcgGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCata                                  <  1:619306/150‑1 (MQ=255)
                                                                                                                      gcgGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCatat                                 <  2:1214032/150‑1 (MQ=255)
                                                                                                                       cgGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATc                                <  2:1387257/150‑1 (MQ=255)
                                                                                                                          ggCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGTGCATATCtta                             >  1:1157637/1‑149 (MQ=255)
                                                                                                                          ggCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTc                             >  1:998958/1‑150 (MQ=255)
                                                                                                                          ggCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTc                             >  1:664274/1‑150 (MQ=255)
                                                                                                                           gCTTGCATTAAAAAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGTACTTTATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTcc                            <  2:937540/150‑1 (MQ=255)
                                                                                                                           gCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTcc                            <  2:155625/150‑1 (MQ=255)
                                                                                                                           gCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTcc                            <  1:560526/150‑1 (MQ=255)
                                                                                                                               gCAGAACACAAAAATGAAACTCTCTGTTTTTTTCATTAATTTTGTGAAATATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGACGTGACTTTTATTGCTGTACAGATTATTTGGTTTTTCAGTGGCATTAAGAGCATATCTTcaagac                        >  2:837495/1‑145 (MQ=255)
                                                                                                                                  gAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCgccgcc                     <  2:39393/150‑1 (MQ=255)
                                                                                                                                       caAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGc                <  1:965066/150‑1 (MQ=255)
                                                                                                                                        attaaTGAAACAAACAGAGTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCAGTGACTTTTATTGCTGTACAGATTATGTGAATTTTAAGTGGCATTAAGTGCATATCTTCCCGCCGCCTCTGCa               >  2:261029/4‑150 (MQ=255)
                                                                                                                                           aaTGAAACAAACAGAGTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGc                                     <  2:67932/125‑1 (MQ=255)
                                                                                                                                           aaTGAAACAAACAGAGTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGc                                     >  1:67932/1‑125 (MQ=255)
                                                                                                                                            aTGATACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTcc           <  1:1368840/150‑1 (MQ=255)
                                                                                                                                            aTGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTcc           <  1:1258471/150‑1 (MQ=255)
                                                                                                                                            aTGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCGCTGCATTcc           <  1:1561939/150‑1 (MQ=255)
                                                                                                                                              gAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTg         <  2:1486134/150‑1 (MQ=255)
                                                                                                                                              gAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTg         <  1:379171/150‑1 (MQ=255)
                                                                                                                                               aaaCACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGt        >  1:207545/1‑150 (MQ=255)
                                                                                                                                                aaCACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCGTGTa       >  2:1533013/1‑150 (MQ=255)
                                                                                                                                                aaCACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTAAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTa       <  1:162376/150‑1 (MQ=255)
                                                                                                                                                   acTCTGTTTGTTTCATTAATTTTGTGAACTATATCAAA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGa    >  1:1566344/1‑150 (MQ=255)
                                                                                                                                                     tctGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCAAGCCGCCTCTGCATTCCTGTAGGaaa  >  1:144464/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                         
CCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAA  >  NC_000913/4275003‑4275301

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