Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,275,151 C→A 5.1% intergenic (‑110/‑320) 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,1510CA 5.1% 406.8 / 8.2 136intergenic (‑110/‑320)yjcB/yjcCputative inner membrane protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (53/76);  new base A (2/5);  total (55/81)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

Only 100 of 140 total aligned reads displayed.

CCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGA  >  NC_000913/4275003‑4275299
                                                                                                                                                    |                                                                                                                                                    
ccAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACAct                                                                                                                                                     >  2:344597/1‑150 (MQ=255)
 cAGCGAAGCAGAACCACGCCAGTGGTGAGGGTATCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGCGGCTGGCAGAACTCAAAAATGCAACActc                                                                                                                                                    >  1:936393/1‑150 (MQ=255)
 cAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACCGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACActc                                                                                                                                                    >  2:698191/1‑150 (MQ=255)
    cGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgt                                                                                                                                                 >  1:1439730/1‑150 (MQ=255)
     gaagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgtt                                                                                                                                                >  2:453689/1‑150 (MQ=255)
     gaagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgtt                                                                                                                                                <  2:706607/150‑1 (MQ=255)
      aagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgttt                                                                                                                                               <  2:1208253/150‑1 (MQ=255)
       agaagaaCCACGCCAGTGGTGAGGGTAGCCATCACACAATTCCTCTGAACCGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgtttg                                                                                                                                              <  2:1272310/150‑1 (MQ=255)
       agaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACCGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgttac                                                                                                                                              >  2:263817/1‑148 (MQ=255)
       agaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgtttg                                                                                                                                              >  1:548839/1‑150 (MQ=255)
          agaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTgtttgttt                                                                                                                                           >  1:659133/1‑150 (MQ=255)
          agaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGCGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAACCACTCTgtttgttt                                                                                                                                           >  1:1063933/1‑150 (MQ=255)
           gaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTc                                                                                                                                          <  1:633592/150‑1 (MQ=255)
            aaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCa                                                                                                                                         <  2:1390129/150‑1 (MQ=255)
             aCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGCCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCAt                                                                                                                                        >  2:1388747/1‑150 (MQ=255)
             aCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCAt                                                                                                                                        <  1:995825/150‑1 (MQ=255)
             aCCACGCCAGTGGTGAGGGTAGACATCAACCAATTCCTCTGAAAACAGCCCATGCCCTGGCCACCGCGGCGTTGTTCACGGGAGGCGTGGTACCCTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACGCAGTTTGTTTAAt                                                                                                                                        >  2:570112/1‑150 (MQ=255)
              ccACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTTGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCAtt                                                                                                                                       <  2:254170/150‑1 (MQ=255)
              ccACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACAAACAGAGTGTTTCAtt                                                                                                                                       <  1:364620/150‑1 (MQ=255)
                  gCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACAAACAGAGTGTTTCAtttttg                                                                                                                                   <  1:980497/150‑5 (MQ=255)
                  gCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACAAACAGAGTGTTTCAtttttg                                                                                                                                   <  2:184608/150‑5 (MQ=255)
                  gCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGAGGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACAAACAGAGTGTTTCAtttttg                                                                                                                                   <  1:1377193/150‑5 (MQ=255)
                    cAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAAtttt                                                                                                                                 >  2:560526/1‑150 (MQ=255)
                     aGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTtg                                                                                                                                >  1:155625/1‑150 (MQ=255)
                        ggtgttcgTAGCCATCAAACAATTCCTCTGACAAGAGCCGATGCCCTGGACCGCGCGGCGTTGTGCACGGGAGGCGTGTTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGa                                                                                                                             <  2:285975/143‑1 (MQ=255)
                        ggtgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCATGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGa                                                                                                                             >  1:1048915/1‑150 (MQ=255)
                           gAGGGTCGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGTACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACt                                                                                                                          <  1:847730/150‑1 (MQ=255)
                                tAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATc                                                                                                                     >  1:361407/1‑150 (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)
                                    cATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACaa                                                                                                                 >  2:342109/1‑150 (MQ=255)
                                      tCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAatt                                                                                                               <  2:65243/150‑1 (MQ=255)
                                         aaCAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAattgat                                                                                                            <  1:1008514/150‑1 (MQ=255)
                                         aaCAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAattgat                                                                                                            <  2:361407/150‑1 (MQ=255)
                                          aCAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAattgatt                                                                                                           <  1:564758/150‑1 (MQ=255)
                                              ttCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGAttgttt                                                                                                       >  1:1651260/1‑150 (MQ=255)
                                              ttCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGAttgttt                                                                                                       <  2:1106376/150‑1 (MQ=255)
                                               tCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGAttgtttg                                                                                                      <  1:375744/150‑1 (MQ=255)
                                                ccTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGAttgtttgt                                                                                                     >  2:379171/1‑150 (MQ=255)
                                                 ctctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATtgtttgtt                                                                                                    <  1:38505/150‑1 (MQ=255)
                                                 ctctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATtgtttgtt                                                                                                    >  2:613471/1‑150 (MQ=255)
                                                     gAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAGCACTCTGTTTGTTTGATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGcc                                                                                                >  1:492063/1‑150 (MQ=255)
                                                     gAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGcc                                                                                                >  1:530428/1‑150 (MQ=255)
                                                     caacagagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGcc                                                                                                <  2:1467885/146‑1 (MQ=255)
                                                       aaaGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCat                                                                                              <  1:453689/150‑1 (MQ=255)
                                                         agagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCatat                                                                                            >  2:1561939/1‑150 (MQ=255)
                                                         agagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCatat                                                                                            >  1:189430/1‑150 (MQ=255)
                                                          gagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATAtt                                                                                           <  2:917966/150‑1 (MQ=255)
                                                             ccGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAgg                                                                                        <  1:451199/150‑1 (MQ=255)
                                                                  gCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTg                                                                                   <  2:607226/150‑1 (MQ=255)
                                                                  gCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTg                                                                                   <  1:1061312/150‑1 (MQ=255)
                                                                    ccTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGAc                                                                                 >  2:965279/1‑150 (MQ=255)
                                                                    ccTGCACAGCGCGGCGTTGTTCACTGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGACACACTCTGTTTGTTTCAGTAAGTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGAc                                                                                 <  1:279208/150‑1 (MQ=255)
                                                                     cTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGATAGCCATATTAGGCCGTGACt                                                                                <  2:527780/150‑1 (MQ=255)
                                                                           aGCGCGCCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTAtt                                                                          <  1:1014348/150‑1 (MQ=255)
                                                                                  cgTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTAc                                                                   <  1:344597/150‑1 (MQ=255)
                                                                                  cgTTGTTCACGCGAGGCGTGGTACACTCGGGCAATCGGGGGGCTTGAAGAACAGAAAAATGAAACACACTGTTTGTTGCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTAc                                                                   <  2:111820/150‑1 (MQ=255)
                                                                                    ttgttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAg                                                                 <  1:1297826/150‑1 (MQ=255)
                                                                                     tgttCACGGGAGGCGTGGTACACTCTGGCTGTCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGa                                                                >  1:354836/1‑150 (MQ=255)
                                                                                      gttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGAt                                                               >  1:1325494/1‑150 (MQ=255)
                                                                                       ttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGAtt                                                              <  2:804899/150‑1 (MQ=255)
                                                                                        tCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACCCTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTa                                                             <  2:1039839/150‑1 (MQ=255)
                                                                                         cACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAt                                                            >  1:1449338/1‑150 (MQ=255)
                                                                                          aCGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtg                                                           <  1:518409/150‑1 (MQ=255)
                                                                                           cGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACTATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtgt                                                          >  2:1592868/1‑150 (MQ=255)
                                                                                           cGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtgt                                                          >  2:1058384/1‑150 (MQ=255)
                                                                                           cGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtgt                                                          >  2:223559/1‑150 (MQ=255)
                                                                                               aGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGtt                                                      <  2:914241/150‑1 (MQ=255)
                                                                                                   gTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTc                                                  <  1:1411029/150‑1 (MQ=255)
                                                                                                     ggTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAg                                                <  1:384545/150‑1 (MQ=255)
                                                                                                     ggTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGATGTACAGATTATGTGGTTTTTCAg                                                >  1:1096319/1‑150 (MQ=255)
                                                                                                      gTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGATGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTg                                              <  2:453489/150‑1 (MQ=255)
                                                                                                       tACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTg                                              <  2:1633583/150‑1 (MQ=255)
                                                                                                       tACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTg                                              <  1:3051/150‑1 (MQ=255)
                                                                                                        acacTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTgg                                             <  2:1546631/150‑1 (MQ=255)
                                                                                                           ctctGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCAt                                          <  2:1285747/150‑1 (MQ=255)
                                                                                                            tctGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGGTTCATTAATTTTGTGAACTATTTCACAATTGATTGTTTGTTAGCCATTTTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCAtt                                         >  2:346659/1‑150 (MQ=255)
                                                                                                             ctGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAAAACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGAATTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTa                                        >  1:1584537/1‑150 (MQ=255)
                                                                                                               ggCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGCCTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAg                                      >  1:331100/1‑150 (MQ=255)
                                                                                                                 cTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAggg                                    <  1:1347231/150‑1 (MQ=255)
                                                                                                                  tATCGCGGGGCTTGCAGAACAATAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGc                                   <  1:119345/150‑1 (MQ=255)
                                                                                                                     cgcgGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCata                                <  1:619306/150‑1 (MQ=255)
                                                                                                                        ggggCTTGCAGAACACAAAACTGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCt                             <  2:587740/150‑1 (MQ=255)
                                                                                                                          ggCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGTGCATATCtta                           >  1:1157637/1‑149 (MQ=255)
                                                                                                                          ggCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTc                           >  1:664274/1‑150 (MQ=255)
                                                                                                                           gCTTGCATTAAAAAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGTACTTTATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTcc                          <  2:937540/150‑1 (MQ=255)
                                                                                                                           gCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTcc                          <  2:155625/150‑1 (MQ=255)
                                                                                                                           gCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTcc                          <  1:560526/150‑1 (MQ=255)
                                                                                                                               gCAGAACACAAAAATGAAACTCTCTGTTTTTTTCATTAATTTTGTGAAATATATCACAATTGATTGTTTGTTAGCCATATTAGGACGTGACTTTTATTGCTGTACAGATTATTTGGTTTTTCAGTGGCATTAAGAGCATATCTTcaagac                      >  2:837495/1‑145 (MQ=255)
                                                                                                                                       caAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGc              <  1:965066/150‑1 (MQ=255)
                                                                                                                                           aaTGAAACAAACAGAGTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGc                                   <  2:67932/125‑1 (MQ=255)
                                                                                                                                           aaTGAAACAAACAGAGTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGc                                   >  1:67932/1‑125 (MQ=255)
                                                                                                                                            aTGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTcc         <  1:1258471/150‑1 (MQ=255)
                                                                                                                                              gAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTTCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTg       <  2:902367/150‑1 (MQ=255)
                                                                                                                                               aaaCACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGt      >  1:207545/1‑150 (MQ=255)
                                                                                                                                                aaCACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCGTGTa     >  2:1533013/1‑150 (MQ=255)
                                                                                                                                                aaCACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTAAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTa     <  1:162376/150‑1 (MQ=255)
                                                                                                                                                   acTCTGTTTGTTTCATTAATTTTGTGAACTATATCAAAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGa  >  1:1566344/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGA  >  NC_000913/4275003‑4275299

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