Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,275,156 T→A 6.2% intergenic (‑115/‑315) 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,1560TA 6.2% 367.5 / 9.3 128intergenic (‑115/‑315)yjcB/yjcCputative inner membrane protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (46/73);  new base A (2/5);  total (49/79)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

Only 100 of 141 total aligned reads displayed.

CGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTA  >  NC_000913/4275007‑4275304
                                                                                                                                                     |                                                                                                                                                        
cGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgt                                                                                                                                                          >  1:1439730/1‑150 (MQ=255)
 gaagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgtt                                                                                                                                                         >  2:453689/1‑150 (MQ=255)
 gaagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGGGCCCGGGAGGCGTGGTACCCTCTGGCGATCGCGGGGCGTGCAGAACACAAAAATGCAACACTCtgtt                                                                                                                                                         >  2:87748/1‑150 (MQ=255)
  aagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACCCTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgttt                                                                                                                                                        >  1:17110/1‑150 (MQ=255)
  aagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgttt                                                                                                                                                        <  2:1208253/150‑1 (MQ=255)
   agaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACCGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgttac                                                                                                                                                       >  2:263817/1‑148 (MQ=255)
   agaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgtttg                                                                                                                                                       >  1:548839/1‑150 (MQ=255)
    gaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAGGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgtttgt                                                                                                                                                      >  1:539392/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)
         aCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCAt                                                                                                                                                 <  1:995825/150‑1 (MQ=255)
         aCCACGCCAGTGGTGAGGGTAGACATCAACCAATTCCTCTGAAAACAGCCCATGCCCTGGCCACCGCGGCGTTGTTCACGGGAGGCGTGGTACCCTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACGCAGTTTGTTTAAt                                                                                                                                                 >  2:570112/1‑150 (MQ=255)
          ccACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTTGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCAtt                                                                                                                                                <  2:254170/150‑1 (MQ=255)
              gCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACAAACAGAGTGTTTCAtttttg                                                                                                                                            <  2:184608/150‑5 (MQ=255)
              gCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGAGGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACAAACAGAGTGTTTCAtttttg                                                                                                                                            <  1:1377193/150‑5 (MQ=255)
                cAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAAtttt                                                                                                                                          >  2:560526/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)
                            tAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATc                                                                                                                              >  1:361407/1‑150 (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)
                                  tCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑att                                                                                                                    <  2:65243/150‑1 (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                                                                                                         >  2:613471/1‑150 (MQ=255)
                                             ctctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATtgtttgtt                                                                                                         <  2:587223/150‑1 (MQ=255)
                                                 gAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAGCACTCTGTTTGTTTGATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGcc                                                                                                     >  1:492063/1‑150 (MQ=255)
                                                 caacagagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGcc                                                                                                     <  2:1467885/146‑1 (MQ=255)
                                                     agagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCatat                                                                                                 >  2:1561939/1‑150 (MQ=255)
                                                     agagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCatat                                                                                                 >  1:189430/1‑150 (MQ=255)
                                                      gagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATAtt                                                                                                <  2:917966/150‑1 (MQ=255)
                                                         ccGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAgg                                                                                             <  1:451199/150‑1 (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)
                                                                              cgTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTAc                                                                        <  1:344597/150‑1 (MQ=255)
                                                                                ttgttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAg                                                                      <  1:1297826/150‑1 (MQ=255)
                                                                                  gttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGAt                                                                    >  1:1325494/1‑150 (MQ=255)
                                                                                     cACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAt                                                                 >  1:1449338/1‑150 (MQ=255)
                                                                                      aCGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtg                                                                <  1:518409/150‑1 (MQ=255)
                                                                                       cGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACT‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtgt                                                               >  2:1592868/1‑150 (MQ=255)
                                                                                       cGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtgt                                                               >  2:223559/1‑150 (MQ=255)
                                                                                       cGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtgt                                                               >  2:1058384/1‑150 (MQ=255)
                                                                                           aGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGtt                                                           <  2:914241/150‑1 (MQ=255)
                                                                                            ggCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGttt                                                          >  1:608225/1‑150 (MQ=255)
                                                                                               gTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTc                                                       <  1:1411029/150‑1 (MQ=255)
                                                                                                 ggTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAg                                                     <  1:384545/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)
                                                                                                    acacTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTgg                                                  <  2:1546631/150‑1 (MQ=255)
                                                                                                     cacTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAg                                                     <  2:680464/150‑1 (MQ=255)
                                                                                                       ctctGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCAt                                               <  2:1285747/150‑1 (MQ=255)
                                                                                                        tctGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGGTTCATTAATTTTGTGAACTATTTCACA‑‑‑‑ATTGATTGTTTGTTAGCCATTTTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCAtt                                              >  2:346659/1‑150 (MQ=255)
                                                                                                           ggCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGCCTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAg                                           >  1:331100/1‑150 (MQ=255)
                                                                                                             cTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAggg                                         <  1:1347231/150‑1 (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)
                                                                                                                  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)
                                                                                                                       gCTTGCATTAAAAAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGTACTTTATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTcc                               <  2:937540/150‑1 (MQ=255)
                                                                                                                       gCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTcc                               <  2:155625/150‑1 (MQ=255)
                                                                                                                           gCAGAACACAAAAATGAAACTCTCTGTTTTTTTCATTAATTTTGTGAAATATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGACGTGACTTTTATTGCTGTACAGATTATTTGGTTTTTCAGTGGCATTAAGAGCATATCTTcaagac                           >  2:837495/1‑145 (MQ=255)
                                                                                                                             aGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCcgccgc                         >  2:768283/1‑150 (MQ=255)
                                                                                                                              gAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCgccgcc                        <  2:39393/150‑1 (MQ=255)
                                                                                                                                   caAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGc                   <  1:965066/150‑1 (MQ=255)
                                                                                                                                       aaTGAAACAAACAGAGTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGc                                        <  2:67932/125‑1 (MQ=255)
                                                                                                                                       aaTGAAACAAACAGAGTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGc                                        >  1:67932/1‑125 (MQ=255)
                                                                                                                                        aTGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTcc              <  1:1258471/150‑1 (MQ=255)
                                                                                                                                        aTGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCGCTGCATTcc              <  1:1561939/150‑1 (MQ=255)
                                                                                                                                          gAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTg            <  1:379171/150‑1 (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‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGaaa     <  1:706562/150‑1 (MQ=255)
                                                                                                                                                 tctGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCAAGCCGCCTCTGCATTCCTGTAGGaaa     >  1:144464/1‑150 (MQ=255)
                                                                                                                                                  ctGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAaat    >  2:842920/1‑150 (MQ=255)
                                                                                                                                                    gtttgtttCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAaatta  >  2:1424678/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                        
CGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACA‑‑‑‑ATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTA  >  NC_000913/4275007‑4275304

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: