Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,311,328 C→G 5.4% A71P (GCA→CCA)  PP_1144 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,311,3280CG 5.4% 271.9 / 4.5 85A71P (GCA→CCA) PP_1144membrane protein
Reads supporting (aligned to +/- strand):  ref base C (36/44);  new base G (2/3);  total (38/47)
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.99e-01

TTGCTGCGTTCGAGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGTTCATGAC  >  NC_002947/1311180‑1311471
                                                                                                                                                    |                                                                                                                                               
ttGCTGCGTTCGAGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGcc                                                                                                                                                >  1:2139470/1‑150 (MQ=255)
 tgctgcGTTCGAGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCa                                                                                                                                               <  1:2040792/150‑1 (MQ=255)
  gctgcGTTCGAGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAg                                                                                                                                              <  2:1173843/150‑1 (MQ=255)
    tgcGTTCGAGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGgc                                                                                                                                            <  1:945599/150‑1 (MQ=255)
        ttCGAGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCTGCCTGGCAGgcg                                                                                                                                           <  1:465564/147‑2 (MQ=255)
        ttCGAGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCTGCCTGGCAGgcg                                                                                                                                           >  2:465564/1‑146 (MQ=255)
        ttCGAGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGgcagca                                                                                                                                        >  1:1037334/1‑150 (MQ=255)
          cGAGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAgg                                                                                                                                      <  2:652533/150‑1 (MQ=255)
            aGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCTGCCTGGCAGgcggcatatg                                                                                                                                    <  1:962811/150‑5 (MQ=255)
            aGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCt                                                                                                                                    <  2:1435534/150‑1 (MQ=255)
            aGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCt                                                                                                                                    >  1:493115/1‑150 (MQ=255)
                tgtgCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCaa                                                                                                                                <  2:2139470/150‑1 (MQ=255)
                tgtgCATGGCCAACCAGGCAACGCCTAGTGCTACGAGCAGTGAAACAATGGTTAGTGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTTCag                                                                                                                                <  1:1289962/150‑2 (MQ=255)
                    cATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCAt                                                                                                                            >  2:647158/1‑150 (MQ=255)
                        gCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTc                                                                                                                        <  1:1349415/150‑1 (MQ=255)
                        gCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTc                                                                                                                        >  1:1157050/1‑150 (MQ=255)
                            aCCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCAtt                                                                                                                    >  1:856976/1‑150 (MQ=255)
                                    aCGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTcggtgg                                                                                                            <  1:411896/150‑1 (MQ=255)
                                       cccAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCCCGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTcggtgggcg                                                                                                         >  2:286916/1‑150 (MQ=255)
                                         cAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTcggtgggcggt                                                                                                       <  1:167426/150‑1 (MQ=255)
                                          aGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTcggtgggcggtg                                                                                                      <  2:493115/150‑1 (MQ=255)
                                           gTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTcggtgggcggtgg                                                                                                     <  2:1039233/150‑1 (MQ=255)
                                           gTGCGATAAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAATGCCAGGATGCTACCACACTGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTcggtgggcggtgg                                                                                                     <  1:1976357/150‑1 (MQ=255)
                                            tGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTcggtgggcggtggg                                                                                                    <  1:780918/150‑1 (MQ=255)
                                            tGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTcggtgggcggtggg                                                                                                    <  1:1662240/150‑1 (MQ=255)
                                                      aGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGgca                                                                                                                                           <  2:673362/101‑1 (MQ=255)
                                                      aGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGgca                                                                                                                                           >  1:673362/1‑101 (MQ=255)
                                                         gAAAGAATGGTTAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGACAATCTGCCGCCTGCCAGGCAGCAGGCTTCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTgg                                                                                       <  1:781646/150‑1 (MQ=255)
                                                          aaaGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGgtc                                                                                                 <  1:2001458/139‑1 (MQ=255)
                                                          aaaGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGgtc                                                                                                 >  2:2001458/1‑139 (MQ=255)
                                                          aaaGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTggg                                                                                      >  1:78743/1‑150 (MQ=255)
                                                           aaGAATGGTTAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGa                                                                                     >  2:2304542/1‑150 (MQ=255)
                                                           aaGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGa                                                                                     >  1:133791/1‑150 (MQ=255)
                                                                tGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGAc                                                                                <  1:1282762/150‑1 (MQ=255)
                                                                tGGTCAGGGATACGCCATAGTGAACGTCCCCCGGGTCCCGAAAGGCCAGCATCCCGACGCAGCGCGTACCTCATAGGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGAc                                                                                <  1:44456/150‑1 (MQ=255)
                                                                       ggATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGc                                                                         >  1:2128252/1‑150 (MQ=255)
                                                                       ggATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGc                                                                         >  1:869936/1‑150 (MQ=255)
                                                                         aTACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCa                                                                       >  2:84481/1‑150 (MQ=255)
                                                                             gtcgtAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCa                                                                                                                      <  2:423090/99‑1 (MQ=255)
                                                                             gtcgtAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCa                                                                                                                      >  1:423090/1‑99 (MQ=255)
                                                                             gtcgtAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTc                                                                                                                 >  2:2202824/1‑104 (MQ=255)
                                                                             gtcgtAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTc                                                                                                                 <  1:2202824/104‑1 (MQ=255)
                                                                             gtcgtAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGgtcgtc                                                                                              <  2:1277824/123‑1 (MQ=255)
                                                                             gtcgtAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGgtcgtc                                                                                              >  1:1277824/1‑123 (MQ=255)
                                                                             gtcgtAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTc                                                                   >  2:1552730/1‑150 (MQ=255)
                                                                             gtcgtAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTc                                                                   >  1:1913065/1‑150 (MQ=255)
                                                                                gtAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCCTGTCTGTTGGACTGACGCTCGGCCATGTCGAg                                                                <  1:674377/150‑1 (MQ=255)
                                                                                 tAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTgg                                                                                       <  1:156857/126‑1 (MQ=255)
                                                                                 tAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTgg                                                                                       >  2:156857/1‑126 (MQ=255)
                                                                                      cACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAgg                                                                                                                                      >  2:1885809/1‑74 (MQ=255)
                                                                                      cACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAgg                                                                                                                                      <  1:1885809/74‑1 (MQ=255)
                                                                                         gTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCgg                                                                          >  1:1558406/1‑131 (MQ=255)
                                                                                         gTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCgg                                                                          <  2:1558406/131‑1 (MQ=255)
                                                                                         gTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGaa                                                       >  2:1017600/1‑150 (MQ=255)
                                                                                            cACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGt                                                                    >  1:280208/1‑134 (MQ=255)
                                                                                            cACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGt                                                                    <  2:280208/134‑1 (MQ=255)
                                                                                             aCCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCCGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAg                                                   <  1:1029628/150‑1 (MQ=255)
                                                                                              ccGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAgg                                                  <  1:286916/150‑1 (MQ=255)
                                                                                                    cccGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTg                                                            >  2:562174/1‑134 (MQ=255)
                                                                                                    cccGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTg                                                            <  1:562174/134‑1 (MQ=255)
                                                                                                     ccGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTc                                                                                                                 <  2:142438/80‑1 (MQ=255)
                                                                                                     ccGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTc                                                                                                                 >  1:142438/1‑80 (MQ=255)
                                                                                                      cGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTgggg                                                                                                   >  2:1758447/1‑93 (MQ=255)
                                                                                                      cGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTgggg                                                                                                   <  1:1758447/93‑1 (MQ=255)
                                                                                                                 cATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATa                               <  2:856976/150‑1 (MQ=255)
                                                                                                                    gCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGc                            <  1:1017600/150‑1 (MQ=255)
                                                                                                                    gCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGc                            <  1:647158/150‑1 (MQ=255)
                                                                                                                      tgatgaAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAg                                                   <  1:404477/125‑1 (MQ=255)
                                                                                                                      tgatgaAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAg                                                   >  2:404477/1‑125 (MQ=255)
                                                                                                                      tgatgaAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGa                          >  2:1381070/1‑150 (MQ=255)
                                                                                                                         ggaAAGGCAAGGCCCACATGCCGCGTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCa                       <  2:133791/149‑1 (MQ=255)
                                                                                                                            aaTGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGaa                    >  2:1340368/1‑150 (MQ=255)
                                                                                                                                 aTGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGAATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAag               >  2:1392253/1‑150 (MQ=255)
                                                                                                                                  tGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAaga              >  1:2165064/1‑150 (MQ=255)
                                                                                                                                  gggCCCATATGCCGGCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTTGGACGGCCGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAaga              <  2:273017/149‑1 (MQ=255)
                                                                                                                                    gCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGc            >  2:2008657/1‑150 (MQ=255)
                                                                                                                                      ccATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCgg          <  1:1340368/150‑1 (MQ=255)
                                                                                                                                        gcctgctgccTGGCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAagag             >  2:355921/8‑145 (MQ=255)
                                                                                                                                        gcctgctgccTGGCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAagag             <  1:355921/138‑1 (MQ=255)
                                                                                                                                        atatGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGtt                                                                                                                         <  1:189001/37‑1 (MQ=255)
                                                                                                                                        atatGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGtt                                                                                                                         >  2:189001/1‑37 (MQ=255)
                                                                                                                                        atatGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAAc                                     <  1:1850610/121‑1 (MQ=255)
                                                                                                                                        atatGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAAc                                     >  2:1850610/1‑121 (MQ=255)
                                                                                                                                        atatGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGtt        >  2:1260176/1‑150 (MQ=255)
                                                                                                                                             ccgccTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGTTCATGa   <  1:1260176/150‑1 (MQ=255)
                                                                                                                                              cgccTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGTTCATGAc  <  2:1785592/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TTGCTGCGTTCGAGGCTGTGCATGGCCAACCAGGCCACGCCCAGTGCGATGAGCAGTGAAAGAATGGTCAGGGATACGTCGTAGTGCACGTCCACCGGGGCCCGAAAGGCCAGCATGCTGATGAAATGCATGGCCCATATGCCGCCTGCCAGGCAGCAGGCTCCAAGCATGTTCCATTGTCGGTGGGCGGTGGGGTCGTCGCTGTGGGACTGACGCTCGGCCATGTCGAGCGTGGCGAAACAGGCCGCGCAAGCAACGAAATAGGCGAGCAGAACCAAGAGCGGTTCATGAC  >  NC_002947/1311180‑1311471

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