Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,546,869 G→A 8.2% intergenic (+25/+15) PP_2238 → / ← PP_5504 metallopeptidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,546,8690GA 8.2% 246.8 / 10.0 87intergenic (+25/+15)PP_2238/PP_5504metallopeptidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/45);  new base A (3/4);  total (38/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

AGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCAGGCGCTTGC  >  NC_002947/2546720‑2547018
                                                                                                                                                      |                                                                                                                                                     
aGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTCGCGCCCCTGCTCCAAGGGGATGCGCTACCCTGGCTGGGGATCCGTACTGCACCCCTGTAACGCAGCCGAAAGTGGGCAGTTTGCg                                                                                                                                                       <  2:1355203/150‑1 (MQ=255)
     cTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTg                                                                                                                                                  >  1:2047727/1‑150 (MQ=255)
       gACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGcc                                                                                                                                                <  1:2179675/150‑1 (MQ=255)
        aCCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGccc                                                                                                                                               >  1:2168841/1‑150 (MQ=255)
         ccAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCt                                                                                                                                              <  2:805791/150‑1 (MQ=255)
         ccAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCt                                                                                                                                              <  1:2048314/150‑1 (MQ=255)
          cAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTCTCTACCGTGTCCTTGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCtt                                                                                                                                             <  2:1871033/150‑1 (MQ=255)
           aaGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCttt                                                                                                                                            <  1:669980/150‑1 (MQ=255)
           aaGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGCCCttt                                                                                                                                            >  1:1888826/1‑150 (MQ=255)
            aggaggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCtttt                                                                                                                                           <  1:248197/150‑1 (MQ=255)
            aggaggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGCCCtttt                                                                                                                                           <  2:1888826/150‑1 (MQ=255)
              gaggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGCCCtttttc                                                                                                                                         <  1:1867475/150‑3 (MQ=255)
                ggaggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTctct                                                                                                                                       >  2:174026/1‑150 (MQ=255)
                   ggTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTAtt                                                                                                                                    >  1:210351/1‑150 (MQ=255)
                    gTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTtc                                                                                                                                          <  1:748801/143‑1 (MQ=255)
                    gTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTtc                                                                                                                                          >  2:748801/1‑143 (MQ=255)
                      ggAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTAc                                                                                                                                 <  2:2388375/150‑1 (MQ=255)
                      ggAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCAACGCTGTAAAGCAGCGGAAAAAGAGCAGTTTGCGAACTGCCCTTTTCTCTATTTAc                                                                                                                                 >  1:1364884/1‑150 (MQ=255)
                       gAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACa                                                                                                                                <  1:1194964/150‑1 (MQ=255)
                       gAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACa                                                                                                                                <  1:1558245/150‑1 (MQ=255)
                         gTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAga                                                                                                                              <  1:1909907/150‑1 (MQ=255)
                         gTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCCTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAga                                                                                                                              <  2:2100319/150‑1 (MQ=255)
                             cTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgat                                                                                                                          >  1:1847976/1‑150 (MQ=255)
                             cTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgat                                                                                                                          <  1:1626330/150‑1 (MQ=255)
                              ttaacttcTACCACGCCTGCGTACATGAGCGCCTGGCGACCTTGCTACAATGGGATGCGCTAGCCTGGCTGGTGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgatg                                                                                                                         <  2:1000027/143‑1 (MQ=255)
                              tGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgatg                                                                                                                         >  1:593158/1‑150 (MQ=255)
                               gAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgatga                                                                                                                        <  2:1847976/150‑1 (MQ=255)
                                aaCGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgatgat                                                                                                                       <  2:210351/150‑1 (MQ=255)
                                          aCGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATc                                                                                                                      >  1:1919175/1‑141 (MQ=255)
                                          aCGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATc                                                                                                                      <  2:1919175/141‑1 (MQ=255)
                                          aCGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCgg                                                                                                             >  2:1446090/1‑150 (MQ=255)
                                              cTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTgg                                                                                                         >  1:1590469/1‑150 (MQ=255)
                                               tGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGgt                                                                                                        >  2:1903646/1‑150 (MQ=255)
                                               tGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGgt                                                                                                        <  1:1554113/150‑1 (MQ=255)
                                                   tACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAg                                                                                                    <  1:93690/150‑1 (MQ=255)
                                                    aCATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGc                                                                                                   <  1:2210634/150‑1 (MQ=255)
                                                      aTGAGCGCCTGGCGCCCTTGCTCCAAGTGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGgt                                                                                                        <  1:2146599/143‑1 (MQ=255)
                                                      aTGAGCGCCTGGCGCCCTTGCTCCAAGTGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGgt                                                                                                        >  2:2146599/1‑143 (MQ=255)
                                                              cTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGTCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAga                                                                                                                              <  1:1213078/113‑1 (MQ=255)
                                                              cTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGTCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAga                                                                                                                              >  2:1213078/1‑113 (MQ=255)
                                                              cTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTa                                                                                                                                      <  1:1820480/105‑1 (MQ=255)
                                                              cTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTa                                                                                                                                      >  2:1820480/1‑105 (MQ=255)
                                                                 ggCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGa                                                                                       <  1:1903646/150‑1 (MQ=255)
                                                                  gcgcCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGAt                                                                                      >  1:1881618/1‑150 (MQ=255)
                                                                    gcCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGTATTTACGAACGTCATCGTTGTTAAGCAGCGGAAAAAGGGCCGCCAGCCATTACTCGGTGGCGCTGTTTACAGATGATGATCATGCTCCGGCTGGTGTAGCACGCCGGGTTGATGc                                                                                    >  1:2350911/1‑150 (MQ=255)
                                                                    gcCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGc                                                                                    <  1:845583/150‑1 (MQ=255)
                                                                    gcCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGATGGTGTAGCCCGCCGGGTTGATGc                                                                                    <  2:2021478/150‑1 (MQ=255)
                                                                     cccTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgatga                                                                                                                        >  1:72896/1‑113 (MQ=255)
                                                                     cccTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAgatgatga                                                                                                                        <  2:72896/113‑1 (MQ=255)
                                                                      ccTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCg                                                                                  >  2:133773/1‑150 (MQ=255)
                                                                            tCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCTGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCACGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACaaa                                                                            <  1:2296532/150‑1 (MQ=255)
                                                                            tCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATg                                                                                                                   <  1:1723103/111‑1 (MQ=255)
                                                                            tCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATg                                                                                                                   >  2:1723103/1‑111 (MQ=255)
                                                                            tCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCgg                                                                                                             >  1:1406810/1‑117 (MQ=255)
                                                                            tCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCgg                                                                                                             <  2:1406810/117‑1 (MQ=255)
                                                                            tCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACaaa                                                                            >  2:2296532/1‑150 (MQ=255)
                                                                              cAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATa                                                                          >  2:855375/1‑150 (MQ=255)
                                                                              cAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATa                                                                          <  2:2047727/150‑1 (MQ=255)
                                                                               aaGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAg                                                                         >  1:141227/1‑150 (MQ=255)
                                                                                aGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGt                                                                        <  1:36639/150‑1 (MQ=255)
                                                                                 ggggATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTc                                                                       <  1:133773/150‑1 (MQ=255)
                                                                                       gcgcTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGc                                                                 <  2:141227/150‑1 (MQ=255)
                                                                                           tAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCCGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTGACCGGCctct                                                             >  1:973163/1‑150 (MQ=255)
                                                                                           tAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCctct                                                             >  2:1066155/1‑150 (MQ=255)
                                                                                             gCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTc                                                           >  1:2141544/1‑150 (MQ=255)
                                                                                              ccTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATc                                                                                                                      >  2:165053/1‑90 (MQ=255)
                                                                                              ccTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATc                                                                                                                      <  1:165053/90‑1 (MQ=255)
                                                                                              ccTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGCCCtt                                                                                                                                             >  2:663531/1‑67 (MQ=38)
                                                                                              ccTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTCGCAAACTGCCCtt                                                                                                                                             <  1:663531/67‑1 (MQ=38)
                                                                                               ctggctggAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCAt                                                         >  2:2234714/1‑150 (MQ=255)
                                                                                                 ggctggAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCg                                                       <  1:1066155/150‑1 (MQ=255)
                                                                                                      gagaTGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATc                                                  <  1:2234714/150‑1 (MQ=255)
                                                                                                            cGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGa                                            <  1:1649866/150‑1 (MQ=255)
                                                                                                             gTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCATTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAg                                           >  1:1474421/1‑150 (MQ=255)
                                                                                                               aCTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTcc                                                                                                               <  1:1400141/80‑1 (MQ=255)
                                                                                                               aCTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTcc                                                                                                               >  2:1400141/1‑80 (MQ=255)
                                                                                                                      cGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAgcgcg                                       >  2:367808/1‑145 (MQ=255)
                                                                                                                      cGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAgcgcg                                       <  1:367808/145‑1 (MQ=255)
                                                                                                                       gCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCATTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGa                                 <  2:1474421/150‑1 (MQ=255)
                                                                                                                         tGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGAcc                               >  2:2064268/1‑150 (MQ=255)
                                                                                                                                 aGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGcc                       >  1:1880639/1‑150 (MQ=255)
                                                                                                                                     gAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCgg                            >  1:310554/1‑141 (MQ=255)
                                                                                                                                     gAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCgg                            <  2:310554/141‑1 (MQ=255)
                                                                                                                                         aaGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGtgc               <  1:370313/150‑1 (MQ=255)
                                                                                                                                         aaGGGCAGTTCGCAAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTa                                                                                                     >  1:770554/1‑64 (MQ=37)
                                                                                                                                         aaGGGCAGTTCGCAAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTa                                                                                                     <  2:770554/64‑1 (MQ=37)
                                                                                                                                               aGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCACGCGCTCGCAATGACCCGGTAGCCGCGCGTGCTGCAgg         <  1:572789/150‑1 (MQ=255)
                                                                                                                                                 tttGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCAGgcg       <  2:1880639/150‑1 (MQ=255)
                                                                                                                                                      gAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCAGGCGCTTGc  >  2:2158237/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCGGAAAAAGGGCAGTTTGCGAACTGCCCTTTTCTCTATTTACAGATGATGATCATGCTCCGGCTGGTGTAGCCCGCCGGGTTGATGCCGAACAAATAGTCACCGGCCTCTTCATCGCTATCGCCCGAGCGCGCAATGACCCGGTAGCCGCGCGTGCTGCAGGCGCTTGC  >  NC_002947/2546720‑2547018

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