Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,871,832 (G)6→7 coding (1448/1476 nt) cdgI → putative c‑di‑GMP binding protein CdgI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,871,8261.G99.5% 857.5 / ‑3.4 208coding (1442/1476 nt)cdgIputative c‑di‑GMP binding protein CdgI
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (94/113);  total (94/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 228 total aligned reads displayed.

ATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTTATTACTATTA  >  NC_000913/1871726‑1871921
                                                                                                     |                                                                                               
aTGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGa                                                                                                  <  1:2440714/101‑1 (MQ=255)
aTGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGa                                                                                                  >  1:2441158/1‑101 (MQ=255)
  gACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAgg                                                                                               >  2:2441160/1‑101 (MQ=255)
     cATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAggggg                                                                                            >  1:2441163/1‑101 (MQ=255)
      aTCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAgggggg                                                                                           >  1:2441165/1‑101 (MQ=255)
        cAGTATGGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAgggggggc                                                                                         <  2:2440797/101‑3 (MQ=255)
          gTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAgggggggcaa                                                                                       >  2:2441168/1‑97 (MQ=255)
            aTTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGcaaca                                                                                      <  1:2440826/101‑1 (MQ=255)
            aTTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGcaaca                                                                                      <  1:2440832/101‑1 (MQ=255)
             ttGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGcaacaa                                                                                     >  1:2441173/1‑101 (MQ=255)
              tGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACaaa                                                                                    <  1:2440931/101‑1 (MQ=255)
              tGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACaaa                                                                                    >  2:2441174/1‑101 (MQ=255)
               ggCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAg                                                                                   <  1:2440812/101‑1 (MQ=255)
                gCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGt                                                                                  >  1:2441176/1‑101 (MQ=255)
                gCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGt                                                                                  <  1:2440688/101‑1 (MQ=255)
                   cGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGat                                                                               >  2:2441178/1‑101 (MQ=255)
                    ggTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGatt                                                                              <  2:2440767/101‑1 (MQ=255)
                    ggTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGatt                                                                              >  2:2441180/1‑101 (MQ=255)
                    ggTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGatt                                                                              >  1:2441181/1‑101 (MQ=255)
                     gTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGatta                                                                             >  1:2441185/1‑101 (MQ=255)
                      ttttttCAACGGGTGACACGCGTAATATCTCGCTTGCCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGattat                                                                            <  1:2440640/101‑1 (MQ=255)
                       tttttCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGattatt                                                                           >  2:2441189/1‑101 (MQ=255)
                       tttttCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGattatt                                                                           <  1:2440602/101‑1 (MQ=255)
                        ttttCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATtc                                                                          <  1:2440926/101‑1 (MQ=255)
                        ttttCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATtc                                                                          <  2:2441072/101‑1 (MQ=255)
                        ttttCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATtc                                                                          >  1:2441191/1‑101 (MQ=255)
                         tttCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATtca                                                                         <  2:2441016/101‑1 (MQ=255)
                          ttCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATtcat                                                                        <  2:2440754/101‑1 (MQ=255)
                            cAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATtcatca                                                                      >  2:2441192/1‑101 (MQ=255)
                             aaCGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATtcat                                                                        >  2:2441195/1‑98 (MQ=255)
                              aCGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGcaaca                                                                                      >  2:2441199/1‑83 (MQ=255)
                              aCGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAg                                                                                   >  1:2441198/1‑86 (MQ=255)
                              aCGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCata                                                                    <  1:2441071/101‑1 (MQ=255)
                               cGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCatat                                                                   >  1:2441202/1‑101 (MQ=255)
                                gggTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATAtt                                                                  >  2:2441203/1‑101 (MQ=255)
                                      cacGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTg                                                            >  2:2441207/1‑101 (MQ=255)
                                        cgcgTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCa                                                          >  2:2441209/1‑101 (MQ=255)
                                         gcgTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCaa                                                         <  2:2440853/101‑1 (MQ=255)
                                           gTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCaaaa                                                       >  1:2441211/1‑101 (MQ=255)
                                            tAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCaaaaa                                                      <  1:2440760/101‑1 (MQ=255)
                                            tAATATCTCGCTTGTCATGACGGAAACAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCaaaaa                                                      <  2:2440902/101‑1 (MQ=255)
                                              atatCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAata                                                    <  2:2440947/101‑1 (MQ=255)
                                              atatCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAata                                                    >  2:2441213/1‑101 (MQ=255)
                                                atCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATAtt                                                  >  1:2441215/1‑101 (MQ=255)
                                                  ctcGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCa                                                <  2:2440941/101‑1 (MQ=255)
                                                   tcGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCag                                               <  2:2440891/101‑1 (MQ=255)
                                                    ctcTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATtcat                                                                        <  2:2440780/73‑1 (MQ=255)
                                                     gCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGt                                                             <  2:2440653/85‑1 (MQ=255)
                                                     gCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCagag                                             >  1:2441218/1‑101 (MQ=255)
                                                     gCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCagag                                             >  1:2441222/1‑101 (MQ=255)
                                                     gCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCagag                                             >  2:2441220/1‑101 (MQ=255)
                                                      cTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGc                                            <  1:2440986/101‑1 (MQ=255)
                                                       ttGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGcc                                           <  2:2440918/101‑1 (MQ=255)
                                                       ttGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGcc                                           >  2:2441226/1‑101 (MQ=255)
                                                       ttGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGcc                                           >  1:2441225/1‑101 (MQ=255)
                                                          tCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATg                                        >  1:2441229/1‑101 (MQ=255)
                                                          tCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATg                                        <  1:2440808/101‑1 (MQ=255)
                                                            aTGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCt                                      >  2:2441231/1‑101 (MQ=255)
                                                            aTGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCt                                      <  2:2440971/101‑1 (MQ=255)
                                                             tGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCGAAAATATTCAGAGCCATGCtt                                     <  1:2440788/101‑1 (MQ=255)
                                                             tGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCtt                                     <  2:2441131/101‑1 (MQ=255)
                                                             tGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCtt                                     <  1:2440844/101‑1 (MQ=255)
                                                              gACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCttt                                    >  2:2441235/1‑101 (MQ=255)
                                                              gACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCttt                                    >  1:2441234/1‑101 (MQ=255)
                                                                cGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCttttt                                  >  1:2441238/1‑101 (MQ=255)
                                                                cGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAggggggg                                                                                          >  1:2441241/1‑43 (MQ=255)
                                                                 ggAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTg                                 >  1:2441243/1‑101 (MQ=255)
                                                                   aaGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCg                               >  2:2441245/1‑101 (MQ=255)
                                                                    aGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGt                              <  1:2440852/101‑1 (MQ=255)
                                                                     gCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTg                             >  1:2441246/1‑101 (MQ=255)
                                                                      cAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTgg                            <  2:2440838/101‑1 (MQ=255)
                                                                        gATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCt                          <  2:2440906/101‑1 (MQ=255)
                                                                         aTAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCtt                         <  2:2440879/101‑1 (MQ=255)
                                                                         aTAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCtt                         >  1:2441250/1‑101 (MQ=255)
                                                                          tAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCttt                        <  2:2440706/101‑1 (MQ=255)
                                                                          tAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCttt                        <  2:2440721/101‑1 (MQ=255)
                                                                             aGCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATAATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGc                     >  2:2441252/1‑101 (MQ=255)
                                                                              gCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCa                    <  2:2440980/101‑1 (MQ=255)
                                                                              gCCTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCa                    >  1:2441254/1‑101 (MQ=255)
                                                                                cTTACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATa                  >  2:2441257/1‑101 (MQ=255)
                                                                                 ttACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATAc                 <  2:2440777/101‑1 (MQ=255)
                                                                                 ttACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATAc                 >  2:2441259/1‑101 (MQ=255)
                                                                                  tACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCa                    <  2:2440824/97‑1 (MQ=255)
                                                                                  tACGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACa                >  1:2441260/1‑101 (MQ=255)
                                                                                   aCGCGAAGCGAAAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATg                                        <  2:2440806/76‑1 (MQ=255)
                                                                                     gcgaagcgaaAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCaa                                                         <  1:2440988/57‑1 (MQ=255)
                                                                                      cgaagcgaaAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAAttt            <  1:2440771/101‑1 (MQ=255)
                                                                                        aagcgaaAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATttat          <  1:2441005/101‑1 (MQ=255)
                                                                                        aagcgaaAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATttat          <  2:2441012/101‑1 (MQ=255)
                                                                                         agcgaaAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATttatt         >  2:2441264/1‑101 (MQ=255)
                                                                                          gcgaaAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATttatta        <  1:2440920/101‑1 (MQ=255)
                                                                                          gcgaaAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATttatta        <  1:2440872/101‑1 (MQ=255)
                                                                                            gaaAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCtttt                       <  1:2440847/84‑1 (MQ=255)
                                                                                            gaaAAGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTtattact      >  1:2441265/1‑101 (MQ=255)
                                                                                             aaaaGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAatat                                                   <  1:2440746/55‑1 (MQ=255)
                                                                                             aaaaGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCt                                      <  1:2440859/68‑1 (MQ=255)
                                                                                             aaaaGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATttatt         <  2:2441113/97‑1 (MQ=255)
                                                                                              aaaGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTtattactat    <  2:2441045/101‑1 (MQ=255)
                                                                                               aaGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTtattactatt   <  1:2440911/101‑1 (MQ=255)
                                                                                               aaGCGAGGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTtattactatt   <  1:2440893/101‑1 (MQ=255)
                                                                                                     |                                                                                               
ATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTTATTACTATTA  >  NC_000913/1871726‑1871921

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