Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,127,093 G→A 5.5% R807R (CGG→CGA)  PP_1887 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,127,0930GA 5.5% 220.1 / 2.4 73R807R (CGG→CGA) PP_1887hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/36);  new base A (2/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

GCCTGGCGTGAGGTGGAGCACTCACGGGTCAGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCAATACGGTCTATCGCACCTATGCCGGGCACAGTCGCAATGACCTG  >  NC_002947/2126944‑2127240
                                                                                                                                                     |                                                                                                                                                   
gcCTGGCGTGAGGTGGAGCACTCACGGGTCAGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCgg                                                                                                                                                     <  1:32886/150‑1 (MQ=255)
  cTGGCGTGAGGTGGAGCACTCACGGGTCAGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCt                                                                                                                                                   >  1:409474/1‑150 (MQ=255)
  cTGGCGTGAGGTGGAGCACTCACGGGTCAGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCt                                                                                                                                                   <  2:1537192/150‑1 (MQ=255)
        ggAGGTGGAGCACTCACGGGTCAGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAg                                                                                                                                             <  2:1596599/149‑1 (MQ=255)
          aGGTGGAGCACTCACGGGTCAGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAggg                                                                                                                                           <  2:222871/150‑1 (MQ=255)
              ggAGCACTCACGGGTCAGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCAAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGc                                                                                                                                       <  1:815181/150‑1 (MQ=255)
                 gCACTCACGGGTCAGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAg                                                                                                                                             >  1:1027155/1‑141 (MQ=255)
                 gCACTCACGGGTCAGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAg                                                                                                                                             <  2:1027155/141‑1 (MQ=255)
                 gCACTCACGGGTCAGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCa                                                                                                                                    <  2:327569/150‑1 (MQ=255)
                     tCACGGGTCAGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCgg                                                                                                                                                     >  1:204989/1‑129 (MQ=255)
                     tCACGGGTCAGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCgg                                                                                                                                                     <  2:204989/129‑1 (MQ=255)
                     tCACGGGTCAGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGa                                                                                                                                >  1:1362330/1‑150 (MQ=255)
                        cGGGTCAGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCGTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGc                                                                                                                             <  1:1090387/150‑1 (MQ=255)
                            tCAGCGGCATCACCGAGACGTGGTACAACCTGTTTGAGAAGCCGGTGCGGATCGGGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTGAATGACTACGATGGAGTTGGCCGGGTGGGTAGGGAAGCTCAAGGAAGCGGTg                                                                                                                         >  2:253207/1‑150 (MQ=255)
                              aGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCa                                                                                                                       <  1:206040/150‑1 (MQ=255)
                               gCGGCATCACCGAGACGTGGTACAACCTGTTGGAGCAGCCGGTGCGGATCGAGCGATTTGCCCTGGACGGAAGTAACAGCATCAGCTTGCGGGGCAAGGACTACGATGGGCTTGGCCGGCGGAGTAGGGAAGCTCAGGGAGGCGGTGCAg                                                                                                                      >  2:704717/1‑150 (MQ=255)
                                     tCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGc                                                                                                                >  1:757506/1‑150 (MQ=255)
                                             aCGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGtt                                                                                                        <  2:336083/150‑1 (MQ=255)
                                             aCGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGtt                                                                                                        <  2:254243/150‑1 (MQ=255)
                                               gTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGa                                                                                                      <  1:203090/150‑1 (MQ=255)
                                               gTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGa                                                                                                      >  1:910516/1‑150 (MQ=255)
                                                 ggTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGt                                                                                                    <  2:1097789/150‑1 (MQ=255)
                                                    aCAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATc                                                                                                 >  1:1472174/1‑150 (MQ=255)
                                                        ccTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAgg                                                                                                                                 <  2:504771/114‑1 (MQ=255)
                                                        ccTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAgg                                                                                                                                 >  1:504771/1‑114 (MQ=255)
                                                         cTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGa                                                                                                      <  1:857855/140‑1 (MQ=255)
                                                         cTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGa                                                                                                      >  2:857855/1‑140 (MQ=255)
                                                               gagaAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCAAACGAACGGGCGTGGAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGGCAGCGTGAGGTTGAGTATCGTTACGATGCt                                                                                      >  2:272368/1‑150 (MQ=255)
                                                                            ggATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTcc                                                                         >  1:189708/1‑150 (MQ=255)
                                                                            ggATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTcc                                                                         >  1:1626707/1‑150 (MQ=255)
                                                                              aTCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGt                                                                       >  1:287524/1‑150 (MQ=255)
                                                                                 gAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAg                                                                    <  2:834482/150‑1 (MQ=255)
                                                                                   gCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGa                                                                                                      >  2:1705542/1‑114 (MQ=255)
                                                                                   gCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGa                                                                                                      <  1:1705542/114‑1 (MQ=255)
                                                                                     gATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGt                                                                                                                          >  1:899591/1‑92 (MQ=255)
                                                                                     gATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGt                                                                                                                          <  2:899591/92‑1 (MQ=255)
                                                                                      aTTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAAc                                                                <  1:682828/149‑1 (MQ=255)
                                                                                      aTTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAAc                                                                >  2:682828/1‑149 (MQ=255)
                                                                                           cccTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTgg                                                          <  2:409474/150‑1 (MQ=255)
                                                                                            ccTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGc                                                         <  1:253207/150‑1 (MQ=255)
                                                                                             cTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGcc                                                        >  1:603655/1‑150 (MQ=255)
                                                                                                gatgggAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAgg                                                                                                                                 <  1:1727134/74‑1 (MQ=255)
                                                                                                gatggaAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAgg                                                                                                                                 >  2:1727134/1‑74 (MQ=255)
                                                                                                   ggaAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGa                                                  <  2:1472174/150‑1 (MQ=255)
                                                                                                     aaGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGt                                                <  2:603655/150‑1 (MQ=255)
                                                                                                     aaGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGt                                                <  2:189708/150‑1 (MQ=255)
                                                                                                       gCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCa                                              >  1:1126733/1‑150 (MQ=255)
                                                                                                          ccAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCAATa                                           <  2:363127/150‑1 (MQ=255)
                                                                                                           cAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCAATAc                                          >  2:1024185/1‑150 (MQ=255)
                                                                                                           cAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGACCGAGTCAATAc                                          >  1:1644381/1‑150 (MQ=255)
                                                                                                            aGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCAATACg                                         <  2:757506/150‑1 (MQ=255)
                                                                                                            aGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCAATACg                                         >  1:1179622/1‑150 (MQ=255)
                                                                                                               aTCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAg                                                                                                                      <  1:1652215/70‑1 (MQ=255)
                                                                                                               aTCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAg                                                                                                                      >  2:1652215/1‑70 (MQ=255)
                                                                                                                tCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCAATACGGTCt                                     >  1:1160162/1‑150 (MQ=255)
                                                                                                                  aGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGACTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTACATCCGGTcgag                                   >  2:138110/1‑147 (MQ=255)
                                                                                                                  aGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGACTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTAAATACGGCCTAt                                   >  2:1107854/1‑150 (MQ=255)
                                                                                                                   gCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGAc                                                                                >  2:308384/1‑104 (MQ=255)
                                                                                                                   gCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGAc                                                                                <  1:308384/104‑1 (MQ=255)
                                                                                                                         aGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCAATAAGGTCTATcaaaact                            >  2:1078024/1‑144 (MQ=255)
                                                                                                                         aGGTCAATGACTACGATGGGCTTGGCCGACTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGTCCCTCTCCGTGAGGAACGGCTGGCCGATCGAGTCAATACGGTCTATCGCACCt                            <  1:1107854/150‑1 (MQ=255)
                                                                                                                         aGGTCAATGACTACGATGGGCTTGGCCGACTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTCCAAACGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCAATACGGTCTATCGCACCt                            <  1:138110/150‑1 (MQ=255)
                                                                                                                                tGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTCGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGCGGAACGGCTGGCCGATCGAGTCAATACGGTCTATCGCACCTATGCCgg                     <  1:679746/150‑1 (MQ=255)
                                                                                                                                 gACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCAATACGGTCTATCGCACCTATGCCggg                    <  1:1024185/150‑1 (MQ=255)
                                                                                                                                  aCTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCAATACGGTCTATCGCACCTATGCCGGGc                   >  2:536397/1‑150 (MQ=255)
                                                                                                                                   cTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCAATACGGTCTATCGCACCTATGCCGGGca                  <  2:1635894/150‑1 (MQ=255)
                                                                                                                                     aCGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCAATACGGTCTATCGCACCTATGCCGGGcaca                <  1:536397/150‑1 (MQ=255)
                                                                                                                                       gATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATc                                                                                                 <  1:1756828/67‑1 (MQ=255)
                                                                                                                                       gATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATc                                                                                                 >  2:1756828/1‑67 (MQ=255)
                                                                                                                                       gATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGt                                                <  2:792872/116‑1 (MQ=255)
                                                                                                                                       gATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGt                                                >  1:792872/1‑116 (MQ=255)
                                                                                                                                        aTGGGCTTGGCCGGCTGATTATGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCAATACGGTCTATCGCACCTATGCCGGGCACAGTc             >  2:1483555/1‑150 (MQ=255)
                                                                                                                                        aTGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCgtgt                                                                           <  1:1508940/88‑1 (MQ=255)
                                                                                                                                        aTGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCgtgt                                                                           >  2:1508940/1‑88 (MQ=255)
                                                                                                                                            gCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATACTCTTCAGCGTTGCCGTGAGTAACGGCTGGCCGATCGAGTCAATACGGTCTATCGCACCTATGCCGGGCACAGTCGCaa         >  1:1204511/1‑150 (MQ=255)
                                                                                                                                                   cGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCAATACGGTCTATCGCACCTATGCCGGGCACAGTCGCAATGACCTg  <  2:1160162/150‑1 (MQ=255)
                                                                                                                                                    ggCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGa                                                  <  2:32511/101‑1 (MQ=255)
                                                                                                                                                    ggCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGa                                                  >  1:32511/1‑101 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCCTGGCGTGAGGTGGAGCACTCACGGGTCAGCGGCATCACCGAGACGTGGTACAACCTGTTCGAGAAGCCGGTGCGGATCGAGCGATTTGCCCTGGATGGAAGCACCAGTATCAGCTTGCAGGTCAATGACTACGATGGGCTTGGCCGGCTGAGTAGGGAAGCTCAAGGAAGCGGTGCAGGCCAGCGTGTGGTTGAGTATCGTTACGATGCTTTCGACCGTGTCCGTGAGGAACGGCTGGCCGATCGAGTCAATACGGTCTATCGCACCTATGCCGGGCACAGTCGCAATGACCTG  >  NC_002947/2126944‑2127240

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: ‑