Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,893,249 G→C 5.1% intergenic (‑14/‑118) yoaA ← / → yoaB putative 5' to 3' DNA helicase implicated in DNA repair/RutC family protein YoaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,893,2490GC 5.1% 309.3 / 4.3 90intergenic (‑14/‑118)yoaA/yoaBputative 5' to 3' DNA helicase implicated in DNA repair/RutC family protein YoaB
Reads supporting (aligned to +/- strand):  ref base G (49/36);  new base C (3/2);  total (52/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTTCCACCACCAGCGGCTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTATCGTTCGTATCGATGCTGAAGCCCG  >  NC_000913/1893101‑1893398
                                                                                                                                                    |                                                                                                                                                     
tGCTTCCACCACCAGCGGCTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCATACTACCCCCTGTTGa                                                                                                                                                      <  2:300855/150‑1 (MQ=255)
tGCTTCCACCACCAGCGGCTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGa                                                                                                                                                      >  2:155650/1‑150 (MQ=255)
  cTTCCACCACCAGCGGCTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAACGCCTGGTCTCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCCCCCTACCCCCGGTTGagt                                                                                                                                                    >  1:929322/1‑148 (MQ=255)
  cTTCCACCACCAGCGGCTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGAtt                                                                                                                                                    <  2:886852/150‑1 (MQ=255)
    tCCACCACCAGCGGCTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTg                                                                                                                                                  <  2:275257/150‑1 (MQ=255)
     ccaccaccaGCGGCTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGa                                                                                                                                                 >  2:276610/1‑150 (MQ=255)
     ccaccaccaGCGGCTGGCCTTTTTCTATCGCCTGGGTGACGACTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGa                                                                                                                                                 >  2:533315/1‑150 (MQ=255)
      caccaccaGCGGCTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGaa                                                                                                                                                >  2:1307733/1‑150 (MQ=255)
      caccaccaGCGGCTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGaa                                                                                                                                                >  2:1046416/1‑150 (MQ=255)
         caccaGCGGCTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAg                                                                                                                                             <  2:368396/150‑1 (MQ=255)
           ccaGCGGCTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAggg                                                                                                                                           <  2:733930/150‑1 (MQ=255)
            caGCGGCTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTCAAATCAACAgggg                                                                                                                                          <  1:855927/150‑2 (MQ=255)
             aGCGGCTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGAt                                                                                                                                         >  2:275872/1‑150 (MQ=255)
              gCGGCTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGAtt                                                                                                                                        <  2:284647/150‑1 (MQ=255)
                ggCTGGCCATTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTAGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTAt                                                                                                                                      >  1:1169857/1‑150 (MQ=255)
                  cTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTCAAATCAACAgggggtagtg                                                                                                                                    <  1:351584/150‑8 (MQ=255)
                      ccTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAgg                                                                                                                                >  2:983787/1‑150 (MQ=255)
                       cTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGa                                                                                                                               >  1:431678/1‑150 (MQ=255)
                        tttttCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGAt                                                                                                                              >  2:1367736/1‑150 (MQ=255)
                           ttCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTAACCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAg                                                                                                                           >  2:966970/1‑150 (MQ=255)
                                 cGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCt                                                                                                                     <  2:868030/150‑1 (MQ=255)
                                  gCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCTGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCtt                                                                                                                    >  2:1264448/1‑150 (MQ=255)
                                  gCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCtt                                                                                                                    >  1:1138267/1‑150 (MQ=255)
                                  gCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCtt                                                                                                                    >  2:4252/1‑150 (MQ=255)
                                    cTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTc                                                                                                                  <  1:276610/150‑1 (MQ=255)
                                     tGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCg                                                                                                                 >  2:15806/1‑150 (MQ=255)
                                     tGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCg                                                                                                                 >  2:1392227/1‑150 (MQ=255)
                                      gggTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGc                                                                                                                >  2:647981/1‑150 (MQ=255)
                                       ggTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGcc                                                                                                               <  1:1307733/150‑1 (MQ=255)
                                       ggTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGcc                                                                                                               <  1:745229/150‑1 (MQ=255)
                                       ggTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGcc                                                                                                               <  1:750196/150‑1 (MQ=255)
                                       ggTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGACTGGTATCGCTTTCGCCAGCTGACAGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGcc                                                                                                               >  2:680979/1‑150 (MQ=255)
                                            cGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTg                                                                                                          <  1:1367736/150‑1 (MQ=255)
                                            cGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATCATGTCAGGATGAGGGCGCTTTCGCCAGTTg                                                                                                          >  2:395784/1‑150 (MQ=255)
                                            cGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCAGTCACACTACCCCCTGTTGATTTGAACAGGGATCATGTCAGGATGAGGGCGCTTTCGCCAGTTg                                                                                                          >  2:400648/1‑150 (MQ=255)
                                             ggCTACCGCCATCTGTCGCTGTGGTTCTCTCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGa                                                                                                         >  1:213217/1‑150 (MQ=255)
                                              gCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGaa                                                                                                        >  1:1243162/1‑150 (MQ=255)
                                               cTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAg                                                                                                       >  2:17626/1‑150 (MQ=255)
                                                        tCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTTGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACgg                                                                                              >  2:466111/1‑150 (MQ=255)
                                                          tGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCg                                                                                            >  1:1199361/1‑150 (MQ=255)
                                                          tGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCg                                                                                            <  1:1046416/150‑1 (MQ=255)
                                                               cTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTc                                                                                       <  2:1172122/150‑1 (MQ=255)
                                                               cTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTc                                                                                       <  2:729330/150‑1 (MQ=255)
                                                               cTGTGGTTATCGAGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACAGTCTGGTGAAAAATCGTCCGTCACACTACCCACTGTTCATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTc                                                                                       >  2:1134662/1‑150 (MQ=255)
                                                                   ggTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACAt                                                                                   >  1:879610/1‑150 (MQ=255)
                                                                    gTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACAtt                                                                                  >  1:932433/1‑150 (MQ=255)
                                                                       ctcGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATtgtg                                                                               >  2:888246/1‑150 (MQ=255)
                                                                        tcGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCTGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTgg                                                                              <  1:1264448/150‑1 (MQ=255)
                                                                         cgcgGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGc                                                                             >  2:21936/1‑150 (MQ=255)
                                                                         cgcgGCTTAAAGCCGGGTATCGCTTTCGCCAGCTGCCCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGc                                                                             <  2:929322/150‑1 (MQ=255)
                                                                           cgGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTTGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAg                                                                           <  1:466111/150‑1 (MQ=255)
                                                                             gCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTc                                                                         >  2:790666/1‑150 (MQ=255)
                                                                                   aGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTgca                                                                   <  2:314078/150‑1 (MQ=255)
                                                                                      cTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTgcagca                                                                >  1:692872/1‑150 (MQ=255)
                                                                                          tATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTAGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGCTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGaa                                                            >  2:723433/1‑150 (MQ=255)
                                                                                             cGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAAt                                                         >  2:371027/1‑150 (MQ=255)
                                                                                             cGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAAt                                                         >  1:861902/1‑150 (MQ=255)
                                                                                             cGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAAt                                                         >  1:341090/1‑150 (MQ=255)
                                                                                              gCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACGCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATg                                                        <  1:680979/150‑1 (MQ=255)
                                                                                               cTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATgg                                                       <  2:341090/150‑1 (MQ=255)
                                                                                                  tCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAg                                                    <  1:647981/150‑1 (MQ=255)
                                                                                                   cGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGt                                                   >  2:100521/1‑150 (MQ=255)
                                                                                                        gCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAAt                                              <  1:275872/150‑1 (MQ=255)
                                                                                                           gACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAg                                           <  1:983787/150‑1 (MQ=255)
                                                                                                              cGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGAc                                        <  2:1138267/150‑1 (MQ=255)
                                                                                                                  tGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACaaaa                                    >  2:52306/1‑150 (MQ=255)
                                                                                                                   ggTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGCCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGCACAAAat                                   >  2:277344/1‑150 (MQ=255)
                                                                                                                    gTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAatg                                  >  2:278270/1‑150 (MQ=255)
                                                                                                                    gTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAatg                                  >  1:1287556/1‑150 (MQ=255)
                                                                                                                      gCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAatgat                                >  2:767355/1‑150 (MQ=255)
                                                                                                                        aaaaTCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAAtgatga                              <  1:533315/150‑1 (MQ=255)
                                                                                                                          aaTCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACt                            <  2:861902/150‑1 (MQ=255)
                                                                                                                              gtccgtcACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTAtcg                        <  2:431678/150‑1 (MQ=255)
                                                                                                                                 cgtcACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTAtcgttc                     >  2:354675/1‑150 (MQ=255)
                                                                                                                                  gtcACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTAtcgttcg                    <  2:932433/150‑1 (MQ=255)
                                                                                                                                     acacTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTAtcgt                       >  2:322200/1‑144 (MQ=255)
                                                                                                                                     acacTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTAtcgt                       <  1:322200/144‑1 (MQ=255)
                                                                                                                                      cacTACCCCCTGTTGATTTGAACAGTGATTATGTCAGGATGAGTGCGCTTTCTCCAGTTGAAGTGGTCACGGCGACCTCACATTGTGGCAGTATTTGAAGCACGAAAATGGAAGTATAATTAGTAAAAAATGATGACTATCGTTCGTAta                >  2:165318/1‑149 (MQ=255)
                                                                                                                                      cacTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTATCGTTCGTATc                >  2:1369598/1‑150 (MQ=255)
                                                                                                                                      cacTACCCCATGGTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAACTGATGACTATCGTTCGTATc                >  2:915008/1‑150 (MQ=255)
                                                                                                                                       cataatcccTGTTCAAATCAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACATCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTATCGTTGGTATCg               >  2:906647/7‑150 (MQ=255)
                                                                                                                                       acTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTATCGTTCGTATCg               <  1:52306/150‑1 (MQ=255)
                                                                                                                                         tACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTATCGTTCGTATCGAt             <  1:15806/150‑1 (MQ=255)
                                                                                                                                          aCCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCt                                                                        <  2:1197760/90‑1 (MQ=255)
                                                                                                                                          aCCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCt                                                                        >  1:1197760/1‑90 (MQ=255)
                                                                                                                                           cccccTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTATCGTTCGTATCGATGc           <  1:4252/150‑1 (MQ=255)
                                                                                                                                           atcccTGTTCAAATCAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTATCGTTCGTATCGATGc           >  2:956102/3‑150 (MQ=255)
                                                                                                                                              ccTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTATCGTTCGTATCGATGCTGa        <  1:278270/150‑1 (MQ=255)
                                                                                                                                                tGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTATCGTTCGTATCGATGCTGAAg      <  1:155650/150‑1 (MQ=255)
                                                                                                                                                 gTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTATCGTTCGTATCGATGCTGAAGc     >  1:496662/1‑150 (MQ=255)
                                                                                                                                                    gATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAg                                                    >  1:729043/1‑100 (MQ=255)
                                                                                                                                                    gATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAg                                                    <  2:729043/100‑1 (MQ=255)
                                                                                                                                                    gATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTATCGTTCGTATCGATGCTGAAGCCCg  >  1:449317/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCTTCCACCACCAGCGGCTGGCCTTTTTCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCGCTTTCGCCAGCTGACCGTCTGGTGCAAAATCGTCCGTCACACTACCCCCTGTTGATTTGAACAGGGATTATGTCAGGATGAGGGCGCTTTCGCCAGTTGAAGTGGTGACGGCGACCTCACATTGTGGCAGTCTTTGCAGCACGAAAATGGAAGTATAATGAGGACAAAATGATGACTATCGTTCGTATCGATGCTGAAGCCCG  >  NC_000913/1893101‑1893398

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