Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,818,526 A→G 6.3% Y209C (TAC→TGC)  PP_1622 → M23B subfamily metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,818,5260AG 6.3% 186.9 / 2.2 63Y209C (TAC→TGC) PP_1622M23B subfamily metallopeptidase
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base G (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

AGTTTGAATAAAGGCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGT  >  NC_002947/1818378‑1818672
                                                                                                                                                    |                                                                                                                                                  
aGTTTGAATAAAGGCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCCAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTAc                                                                                                                                                   <  1:1127961/150‑1 (MQ=255)
 gTTTGAATAAAGGCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACg                                                                                                                                                  <  1:780400/150‑1 (MQ=255)
  tttGAATAAAGGCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACgg                                                                                                                                                 >  1:881114/1‑150 (MQ=255)
   ttGAATAAAGGCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGt                                                                                                                                                >  1:86787/1‑150 (MQ=255)
       aTAAAGGCATTGACAACGCCGGTGAGCTGGGACAGCCTGCTTTTGCTGCGCCTGATGGCGTGGTGGCTTATGCCGGCTGCGGTGTCAGGGGCTACGGCGAGCTGATCATCATCAAGTACAGCCACACCAACGTCAGCGCCTACggctata                                                                                                                                            >  2:125734/1‑145 (MQ=255)
        tAAAGGCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATaa                                                                                                                                           <  2:1069839/150‑1 (MQ=255)
           aGGCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCg                                                                                                                                        >  2:1485702/1‑150 (MQ=255)
              cATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTGCGGTTATGACCGTAg                                                                                                                                     <  2:647665/150‑1 (MQ=255)
              cATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAg                                                                                                                                     >  2:1434230/1‑150 (MQ=255)
               aTTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTGCGGTTATGACCGTAgg                                                                                                                                    <  1:224349/150‑1 (MQ=255)
                   atatCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGt                                                                                                                                <  2:207592/150‑1 (MQ=255)
                   atatCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGt                                                                                                                                <  1:1568101/150‑1 (MQ=255)
                     atCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTg                                                                                                                              <  2:722835/150‑1 (MQ=255)
                         ccGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTc                                                                                                                          <  1:2384566/150‑1 (MQ=255)
                           ggTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCgg                                                                                                                        <  1:1485702/150‑1 (MQ=255)
                              gATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCGACGGTCATAACCGCAGGCTGTTGGTTCGGGAg                                                                                                                     <  1:1171243/150‑1 (MQ=255)
                                   gggACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGgca                                                                                                                >  2:2232241/1‑150 (MQ=255)
                                          ccTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTc                                                                                                         >  1:451266/1‑150 (MQ=255)
                                           cTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATaa                                                                                                                                           <  2:1533498/115‑1 (MQ=255)
                                           cTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATaa                                                                                                                                           >  1:1533498/1‑115 (MQ=255)
                                           cTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCa                                                                                                        >  2:402702/1‑150 (MQ=255)
                                                  tgctgcGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGTTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggcagg                                                                                                 <  2:881114/150‑1 (MQ=255)
                                                  tgctgcGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggcagg                                                                                                 <  1:2324939/150‑1 (MQ=255)
                                                  tgctgcGTCTGATGGTGCGGTGGTTTATGACGGGAGTTGTGTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggcagg                                                                                                 <  2:266175/150‑1 (MQ=255)
                                                    ctgcGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggcagggc                                                                                               >  1:370263/1‑150 (MQ=255)
                                                           tGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCaa                                                                                                       <  1:2333113/135‑1 (MQ=255)
                                                           tGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCaa                                                                                                       >  2:2333113/1‑135 (MQ=255)
                                                                 tGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGa                                                                                  >  2:1824208/1‑150 (MQ=255)
                                                                     gtggtTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTa                                                                                                                                                    <  1:831119/80‑1 (MQ=255)
                                                                     gtggtTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTa                                                                                                                                                    >  2:831119/1‑80 (MQ=255)
                                                                       ggtTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATggg                                                                            <  2:451266/150‑1 (MQ=255)
                                                                                    aGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGcaca                                                               <  1:1824208/150‑1 (MQ=255)
                                                                                        gTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATc                                                           >  2:1874767/1‑150 (MQ=255)
                                                                                         tttGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCg                                                          >  1:230673/1‑150 (MQ=255)
                                                                                          ttGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCgg                                                         <  1:1434230/150‑1 (MQ=255)
                                                                                              ggggCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCg                                                          <  1:1348530/145‑1 (MQ=255)
                                                                                              ggggCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCg                                                          >  2:1348530/1‑145 (MQ=255)
                                                                                               gggCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGTGTCTACGGGCACAGATCGGGTGaa                                                    <  1:1874767/150‑1 (MQ=255)
                                                                                                  cgaCGGAGAGTGGATCCTCCTGAAGTAAAGCGTGACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCt                                                 <  1:125734/148‑1 (MQ=255)
                                                                                                    aCGGCGAGCTGATCATCATCAAGCAAAGCGACACCTACGTCAGTGCCTACGGTCCTACCCGCAGGCTGATGGCGCGGGAGGGGCAGCAGGGCAGGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCAGAGATCGGGTGAAGCTGc                                               >  1:736673/1‑150 (MQ=255)
                                                                                                      ggCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGa                                                                                                                      >  1:1046950/1‑77 (MQ=255)
                                                                                                      ggCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGa                                                                                                                      <  2:1046950/77‑1 (MQ=255)
                                                                                                              gATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGAt                                                                                        >  2:1471400/1‑99 (MQ=255)
                                                                                                              gATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGAt                                                                                        <  1:1471400/99‑1 (MQ=255)
                                                                                                               atcatcatcAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAAt                                                                               <  2:978809/107‑1 (MQ=255)
                                                                                                               atcatcatcAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAAt                                                                               >  1:978809/1‑107 (MQ=255)
                                                                                                               atcatcatcAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTAc                                                                      >  2:726075/1‑116 (MQ=255)
                                                                                                               atcatcatcAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTAc                                                                      <  1:726075/116‑1 (MQ=255)
                                                                                                                  atcatcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTcgc                                 >  2:1115123/1‑150 (MQ=255)
                                                                                                                     atcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGt                                                                                                                                >  2:1818946/1‑52 (MQ=255)
                                                                                                                     atcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGt                                                                                                                                <  1:1818946/52‑1 (MQ=255)
                                                                                                                                 gACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCACGGCAGGGTAGTCGATTGCTGAAATGGGGTCTACGGCCACAGATCAGTTGAAGCCGCATTTCGAGATTCGCCGCCAAGGCAAAccc                  >  2:451159/1‑150 (MQ=255)
                                                                                                                                 gACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAAccc                  <  2:230673/150‑1 (MQ=255)
                                                                                                                                      cTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGa             >  2:438406/1‑150 (MQ=255)
                                                                                                                                      cTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGAAGCAGGTCAAGGCAGGGCAGTCGATTg                                                                                      >  2:2044684/1‑77 (MQ=255)
                                                                                                                                      cTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGAAGCAGGTCAAGGCAGGGCAGTCGATTg                                                                                      <  1:2044684/77‑1 (MQ=255)
                                                                                                                                       tACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGACATTCGCCGCCAGGGCAAACCCGTCGAt            >  2:533570/1‑150 (MQ=255)
                                                                                                                                           caacagCCTGCGGTTATGACCGTAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCAc        >  2:1682578/6‑150 (MQ=255)
                                                                                                                                               tGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCgg                                                         >  1:2283659/1‑97 (MQ=255)
                                                                                                                                               tGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCgg                                                         <  2:2283659/97‑1 (MQ=255)
                                                                                                                                               tGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCa    >  1:467821/1‑150 (MQ=255)
                                                                                                                                                gCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAg   >  1:1849864/1‑150 (MQ=255)
                                                                                                                                                 ccTGCGGTTATGACCGTAGGCTGTTGGTTCGGGAGGGGCANCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGt  >  1:1425759/1‑150 (MQ=255)
                                                                                                                                                   tACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGt                >  1:792158/1‑134 (MQ=255)
                                                                                                                                                   tACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGt                <  2:792158/134‑1 (MQ=255)
                                                                                                                                                    aCGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAAc                    >  2:1281180/1‑129 (MQ=255)
                                                                                                                                                    aCGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAAc                    <  1:1281180/129‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGTTTGAATAAAGGCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGT  >  NC_002947/1818378‑1818672

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