Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,062,881 G→C 6.8% A585P (GCT→CCT)  PP_1838 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,062,8810GC 6.8% 215.5 / 4.8 72A585P (GCT→CCT) PP_1838membrane protein
Reads supporting (aligned to +/- strand):  ref base G (27/39);  new base C (2/3);  total (29/43)
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.66e-01

AGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCC  >  NC_002947/2062732‑2063026
                                                                                                                                                      |                                                                                                                                                 
agcagGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAgg                                                                                                                                                   <  2:1157754/150‑1 (MQ=255)
                ggACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTg                                                                                                                                   <  1:491450/150‑1 (MQ=255)
                   ccTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTg                                                                                                                                      >  2:1425445/1‑144 (MQ=255)
                   ccTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTg                                                                                                                                      <  1:1425445/144‑1 (MQ=255)
                         ccGCTTCAACGTGCCGCTGCTGCTGCTTGCCCCGGGCATTCAGGAGAATGTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGa                                                                                                                          <  2:867505/150‑1 (MQ=255)
                          cGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGAc                                                                                                                         <  2:964826/150‑1 (MQ=255)
                             ttCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTg                                                                                                                      >  1:1030757/1‑150 (MQ=255)
                             ttCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTg                                                                                                                      >  1:1238972/1‑150 (MQ=255)
                               cAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGCAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCt                                                                                                                    >  1:1358458/1‑150 (MQ=255)
                                   gTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAAc                                                                                                                <  2:173245/150‑1 (MQ=255)
                                    tGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACAGCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAAcc                                                                                                               >  1:455585/1‑150 (MQ=255)
                                      ccGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCACCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTg                                                                                                             >  2:116445/1‑150 (MQ=255)
                                       cGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGc                                                                                                                                                  <  1:256407/112‑1 (MQ=255)
                                       cGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGc                                                                                                                                                  >  2:256407/1‑112 (MQ=255)
                                       cGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGc                                                                                                            <  1:862571/150‑1 (MQ=255)
                                        gctgctgctGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGAc                                                                                                                         <  2:930312/136‑1 (MQ=255)
                                        gctgctgctGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGAc                                                                                                                         >  1:930312/1‑136 (MQ=255)
                                        gctgctgctGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGcc                                                                                                           <  2:1358458/150‑1 (MQ=255)
                                             tgctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAgg                                                                                                      <  1:1200955/150‑1 (MQ=255)
                                              gctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTACTCAACCTGCACGAggg                                                                                                     >  2:364465/1‑150 (MQ=255)
                                                     gCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGACCCCCATGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTCGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAg                                                                                              <  1:1155559/150‑1 (MQ=255)
                                                       cccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAggg                                                                                            >  1:684911/1‑150 (MQ=255)
                                                        ccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGc                                                                                           >  1:415001/1‑150 (MQ=255)
                                                           ggCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGTCTCGCCACCAGTGCTggg                                                                                                                                 <  2:837771/109‑1 (MQ=255)
                                                           ggCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTggg                                                                                                                                 >  1:837771/1‑109 (MQ=255)
                                                            gCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTCAACCTGTTGGCATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAg                                                                                       <  2:37010/150‑1 (MQ=255)
                                                                  aGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGa                                                                                 <  2:455585/150‑1 (MQ=255)
                                                                   ggAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTa                                                                                        <  1:458801/142‑1 (MQ=255)
                                                                   ggAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTa                                                                                        >  2:458801/1‑142 (MQ=255)
                                                                      gaAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGt                                                                                                    <  1:642841/127‑1 (MQ=255)
                                                                      gaAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGt                                                                                                    >  2:642841/1‑127 (MQ=255)
                                                                      gaAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCaa                                                                             <  2:739286/150‑1 (MQ=255)
                                                                        aGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGc                                                                           <  2:301275/150‑1 (MQ=255)
                                                                               tgcggtgAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCCCCAGCTGCGGCGGCCGCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCggg                                                                    <  1:738402/150‑1 (MQ=255)
                                                                                        tCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGTGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCa                                                           >  1:155045/1‑150 (MQ=255)
                                                                                         cacaCCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAg                                                          <  1:729971/150‑1 (MQ=255)
                                                                                           cacCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGAc                                                                                                 >  2:283414/1‑109 (MQ=255)
                                                                                           cacCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGAc                                                                                                 <  1:283414/109‑1 (MQ=255)
                                                                                            acCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGAtt                                                       >  1:1047114/1‑150 (MQ=255)
                                                                                                    tACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCt                                               <  2:1047114/150‑1 (MQ=255)
                                                                                                    tACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCt                                               >  2:52657/1‑150 (MQ=255)
                                                                                                       ccAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGt                                            >  1:1388636/1‑150 (MQ=255)
                                                                                                       ccAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGt                                            >  2:140268/1‑150 (MQ=255)
                                                                                                            tCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAgg                                       <  2:1238972/150‑1 (MQ=255)
                                                                                                                 aTCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCGCCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGAc                                  <  2:1030757/150‑1 (MQ=255)
                                                                                                                  tCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGAcc                                 <  2:1293700/150‑1 (MQ=255)
                                                                                                                  tCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGAcc                                 <  2:1388636/150‑1 (MQ=255)
                                                                                                                  tCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGAcc                                 <  1:52657/150‑1 (MQ=255)
                                                                                                                  tCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATCCTCAAGCCGTCGGGCAGCGAGCCGATTGTCTGGCTGGTGCAGGGTCAcc                                 <  2:415001/150‑1 (MQ=255)
                                                                                                                         gACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAg                                                                            >  2:1061857/1‑100 (MQ=255)
                                                                                                                         gACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAg                                                                            <  1:1061857/100‑1 (MQ=255)
                                                                                                                          aCCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCGCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGAcc                                 <  2:474929/142‑1 (MQ=255)
                                                                                                                          aCCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGAcc                                 >  1:474929/1‑142 (MQ=255)
                                                                                                                            catcatGGGGCGGCTGTTTGGCTAGGCGGGCCACGAGTGCGGGGGGCGCGACCTGCTCCACGTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATGCTGAt                       <  2:582222/150‑1 (MQ=255)
                                                                                                                               catGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGa                    <  2:684911/150‑1 (MQ=255)
                                                                                                                                atGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGc                                                                   <  1:578493/102‑1 (MQ=255)
                                                                                                                                atGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGc                                                                   >  2:578493/1‑102 (MQ=255)
                                                                                                                                     ccgccTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGGCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCaa              <  1:116445/150‑1 (MQ=255)
                                                                                                                                      cgccTGGGTGGCAAGGCTCGCCCCCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGGGTCCAAg             <  1:638347/150‑1 (MQ=255)
                                                                                                                                      cgccTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAg             <  1:369844/150‑1 (MQ=255)
                                                                                                                                         cactggTGGCGAGCCTTGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCaa              <  1:1311294/143‑1 (MQ=255)
                                                                                                                                         cactggTGGCGAGCCTTGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCaa              >  2:1311294/5‑147 (MQ=255)
                                                                                                                                            tggTGGCGAGCCTTGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGt                                                     >  1:946526/2‑105 (MQ=255)
                                                                                                                                            tggTGGCGAGCCTTGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGt                                                     <  2:946526/104‑1 (MQ=255)
                                                                                                                                              gTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGa      >  1:5988/1‑150 (MQ=255)
                                                                                                                                                ggCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGc    <  1:910706/150‑1 (MQ=255)
                                                                                                                                                 gCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCNAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGAt                       >  1:165720/1‑130 (MQ=255)
                                                                                                                                                 gCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTa                                                                                        <  2:105455/65‑1 (MQ=255)
                                                                                                                                                 gCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTa                                                                                        >  1:105455/1‑65 (MQ=255)
                                                                                                                                                 gCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGTCTGGTGCAGGGTGACCGTATCCTGAt                       <  2:165720/130‑1 (MQ=255)
                                                                                                                                                  cAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGccc  <  1:1363616/150‑1 (MQ=255)
                                                                                                                                                   aaGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAggg                                      >  1:831752/1‑113 (MQ=255)
                                                                                                                                                   aaGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAggg                                      <  2:831752/113‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
AGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCC  >  NC_002947/2062732‑2063026

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: