Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,026,541 T→C 10.0% intergenic (+16/+16) helD → / ← mgsA DNA helicase IV/methylglyoxal synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,026,5410TC10.0% 206.0 / 13.5 80intergenic (+16/+16)helD/mgsADNA helicase IV/methylglyoxal synthase
Reads supporting (aligned to +/- strand):  ref base T (35/37);  new base C (4/4);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

AAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCG  >  NC_000913/1026392‑1026690
                                                                                                                                                     |                                                                                                                                                     
aagaaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGttt                                                                                                                                                       >  2:234865/1‑150 (MQ=255)
 agaaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTg                                                                                                                                                      <  1:1235259/150‑1 (MQ=255)
 agaaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTg                                                                                                                                                      <  1:868963/150‑1 (MQ=255)
  gaaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGc                                                                                                                                                     <  2:389647/150‑1 (MQ=255)
  gaaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGc                                                                                                                                                     <  1:1194979/150‑1 (MQ=255)
  gaaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGc                                                                                                                                                     <  1:348180/150‑1 (MQ=255)
   aaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGcc                                                                                                                                                    <  1:2847/150‑1 (MQ=255)
      ggcggTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCTCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCAcc                                                                                                                                                 <  2:1194207/150‑1 (MQ=255)
      ggcggTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCAcc                                                                                                                                                 <  1:1022888/150‑1 (MQ=255)
      ggcggTTAATGGACGGGGCGCTGACCCGGGCACGCCATCGGTTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCAcc                                                                                                                                                 <  2:263824/150‑1 (MQ=255)
        cggTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtg                                                                                                                                               <  1:158679/150‑1 (MQ=255)
         ggTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAATAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCtgt                                                                                                                                              <  2:1157158/150‑1 (MQ=255)
         ggTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtgt                                                                                                                                              <  1:1560731/150‑1 (MQ=255)
          gTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtgtt                                                                                                                                             <  2:1206677/150‑2 (MQ=255)
          gTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtgtt                                                                                                                                             >  1:1206677/1‑149 (MQ=255)
          gTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTCAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCtgtg                                                                                                                                             <  2:829368/150‑1 (MQ=255)
            tAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtgtttc                                                                                                                                           <  2:1245911/150‑4 (MQ=255)
                gTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAatat                                                                                                                                       >  2:1650005/1‑150 (MQ=255)
                    gTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACt                                                                                                                                   >  2:814644/1‑150 (MQ=255)
                      ggCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTc                                                                                                                                 <  1:939129/150‑1 (MQ=255)
                          cTGACCCGGGCACGCCATCGGGACTGGGCACTGAATACCCAAGAGAATCCCTCTCCCTTTGTGGAAATACAGAAAAATCTGGATGTGCCGGCGGCGAGAAAACCGTAAGAAACAGGTGGCGCTTGCCACCTGTGCAATATTACTTCAGAc                                                                                                                             <  1:315800/150‑1 (MQ=255)
                             aCCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGt                                                                                                                          >  2:1384834/1‑150 (MQ=255)
                                 gggCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCcgc                                                                                                                      >  2:808213/1‑150 (MQ=255)
                                      cGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGAt                                                                                                                 >  2:356941/1‑150 (MQ=255)
                                       gCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATa                                                                                                                >  2:684529/1‑150 (MQ=255)
                                        ccATCGGGTATGGGCACTGTTTAACAATGATAATCCCTCTCCCTTTGTGGAATTACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGGAATAAACAGGTGGCGTTTGCCAGCTGTGCAATATTAATTAAGACGGTCCGCGATATaa                                                                                                               >  2:943446/1‑150 (MQ=255)
                                         cATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAAc                                                                                                              <  1:1485406/150‑1 (MQ=255)
                                              gggaTGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTCAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGAAATATTACTTCAGACGGTCCGCGAGATAACGCTGa                                                                                                         <  1:1073875/147‑1 (MQ=255)
                                              ggTATGTGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGa                                                                                                         <  1:469643/150‑1 (MQ=255)
                                               gTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGAt                                                                                                        <  2:1376473/150‑1 (MQ=255)
                                                   gggCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATc                                                                                                    <  1:573831/150‑1 (MQ=255)
                                                     gCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCAACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCgg                                                                                                  >  1:203636/1‑150 (MQ=255)
                                                              aaCAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGaa                                                                                         >  1:1538891/1‑150 (MQ=255)
                                                               aCAAATAGAAGCCCTCTCCCTTTGTGGAAATACTGCAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTTCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAat                                                                                        <  2:1077223/150‑1 (MQ=255)
                                                                cAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAata                                                                                       >  2:1564189/1‑150 (MQ=255)
                                                                 aaaGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAatat                                                                                      <  1:695239/150‑1 (MQ=255)
                                                                 aaaGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGAATCAGAatat                                                                                      >  1:78914/1‑150 (MQ=255)
                                                                      gaATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGccc                                                                                 >  2:1267754/1‑150 (MQ=255)
                                                                        aTCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACcgc                                                                               >  1:1611796/1‑150 (MQ=255)
                                                                        aTCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACcgc                                                                               >  2:943278/1‑150 (MQ=255)
                                                                          cccTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGcgt                                                                             <  1:547699/150‑1 (MQ=255)
                                                                           ccTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAACCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGcgtc                                                                            >  1:343286/1‑150 (MQ=255)
                                                                            ctctCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGcgtcg                                                                           <  1:808213/150‑1 (MQ=255)
                                                                               tcCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTg                                                                        >  1:803887/1‑150 (MQ=255)
                                                                                cccTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGa                                                                       <  2:1140713/150‑1 (MQ=255)
                                                                                cccTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGa                                                                       >  2:604413/1‑150 (MQ=255)
                                                                                       tgGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATgcgg                                                                <  1:234865/150‑1 (MQ=255)
                                                                                        ggAAATACTGAAAAATCTGGATGTGCAGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCTGGGATCAGAATATCGACCGCGTCGTTGAAATgcggc                                                               >  2:1253918/1‑150 (MQ=255)
                                                                                            aTACTGAAAACTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACt                                                           <  2:60624/150‑1 (MQ=255)
                                                                                               cTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGa                                                        >  2:842375/1‑150 (MQ=255)
                                                                                                 taaaaaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCTCCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATAGGGGAGCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAtt                                                      <  1:790156/149‑1 (MQ=255)
                                                                                                 gAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAtt                                                      >  2:1315732/1‑150 (MQ=255)
                                                                                                   aaaaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTAt                                                    >  2:257563/1‑150 (MQ=255)
                                                                                                     aaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGa                                                  <  2:504562/150‑1 (MQ=255)
                                                                                                     aaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGa                                                  <  1:547013/150‑1 (MQ=255)
                                                                                                         tGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGta                                              >  2:590141/1‑149 (MQ=255)
                                                                                                          ggATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCt                                             >  2:631357/1‑150 (MQ=255)
                                                                                                           gATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCATTGAAATGCGGCAACTGGATTATGAAGTCTg                                            <  2:771203/150‑1 (MQ=255)
                                                                                                            aTGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGc                                           >  1:1068557/1‑150 (MQ=255)
                                                                                                            aTGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGc                                           <  1:814644/150‑1 (MQ=255)
                                                                                                             tgtgCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCgggg                                                                                                >  1:710668/1‑96 (MQ=255)
                                                                                                             tgtgCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCgggg                                                                                                <  2:710668/96‑1 (MQ=255)
                                                                                                                gCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTACGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACAGCGTCGTTGAAATTCGGCGACTGGATTATGAATTCTTCCGtt                                       >  2:755993/1‑150 (MQ=255)
                                                                                                                 ccGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTg                                      >  2:1323672/1‑150 (MQ=255)
                                                                                                                   ggTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGcc                                    >  2:62925/1‑150 (MQ=255)
                                                                                                                     tGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCAc                                  >  2:5093/1‑150 (MQ=255)
                                                                                                                      ggCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACg                                 <  1:1576855/150‑1 (MQ=255)
                                                                                                                       gCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGt                                >  1:1700235/1‑150 (MQ=255)
                                                                                                                        ggagaAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGtt                               <  2:1068557/149‑1 (MQ=255)
                                                                                                                        cGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGtt                               >  1:893055/1‑150 (MQ=255)
                                                                                                                          agaAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtgg                             <  2:210878/150‑1 (MQ=255)
                                                                                                                              aaCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtggtgg                          >  2:1439120/1‑149 (MQ=255)
                                                                                                                              aaCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtggtgg                          <  1:1439120/149‑1 (MQ=255)
                                                                                                                                 cGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGAc                      >  2:1132774/1‑150 (MQ=255)
                                                                                                                                       gcaCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGCCCGCGTCGTTGAAATGCGGCGACTGGATTAGGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAAt                >  2:120516/3‑150 (MQ=255)
                                                                                                                                        caCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATg               >  1:507357/2‑150 (MQ=255)
                                                                                                                                        aaCAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATg               >  1:590616/1‑150 (MQ=255)
                                                                                                                                         aCAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGt              >  1:632592/1‑150 (MQ=255)
                                                                                                                                         aCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGt              >  2:401195/1‑150 (MQ=255)
                                                                                                                                            ggTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTcc           <  1:1313039/150‑1 (MQ=255)
                                                                                                                                                 cGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATAcc      <  1:1564189/150‑1 (MQ=255)
                                                                                                                                                     tGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCg  <  2:803887/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCG  >  NC_000913/1026392‑1026690

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: