Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,844,152 A→G 6.6% T240A (ACA→GCA)  PP_5122 → GMC family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,844,1520AG 6.6% 339.2 / 8.6 117T240A (ACA→GCA) PP_5122GMC family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (52/56);  new base G (5/3);  total (58/59)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 119 total aligned reads displayed.

TGGGCTACTGCGGCATGGGCTGCCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGG  >  NC_002947/5844003‑5844300
                                                                                                                                                     |                                                                                                                                                      
tGGGCTACTGCGGCATGGGCTGCCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGa                                                                                                                                                        <  2:805114/150‑1 (MQ=255)
        tGCGGCATGGGCTGCCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAg                                                                                                                                              <  2:86390/150‑1 (MQ=255)
          cGGCATGGGCTGCCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAAATGCTGTATCTGGCGCGGGCCGGGCGCTTCGCTCACAACGGCGAATGCATTCAGGGCCTTGCA‑TG‑CCAGGc                                                                                                                                            >  1:938221/1‑150 (MQ=255)
                gggCTGCCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGa                                                                                                                                      <  1:1765507/150‑1 (MQ=255)
                gggCTGCCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGa                                                                                                                                      <  1:73747/150‑1 (MQ=255)
                gggCTGCCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGa                                                                                                                                      <  1:2451273/150‑1 (MQ=255)
                   cTGCCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGc                                                                                                                                   <  1:704635/150‑1 (MQ=255)
                       ccGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAg                                                                                                                               <  2:1216805/150‑1 (MQ=255)
                        cGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAgg                                                                                                                              <  1:660341/150‑1 (MQ=255)
                        cGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAgg                                                                                                                              <  1:180976/150‑1 (MQ=255)
                           tGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCa                                                                                                                           <  1:2295820/150‑1 (MQ=255)
                             aaTGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATa                                                                                                                         <  1:1260652/150‑1 (MQ=255)
                             aaTGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATa                                                                                                                         <  1:1260414/150‑1 (MQ=255)
                               tGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATaca                                                                                                                       <  1:2127088/150‑1 (MQ=255)
                               tGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACCAATACCGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGGCGTCCAGGGcctggc                                                                                                                       <  1:1637366/150‑6 (MQ=255)
                                  cAAACAATCGATGCTGGTGACCCGTATACCTGCCAACATAGCGCAAGGGGACGAACATCTGTACCTGGCGCTTGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACAccc                                                                                                                    >  1:647351/1‑150 (MQ=255)
                                   aaaCAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCt                                                                                                                   >  1:110256/1‑150 (MQ=255)
                                     aCAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTAc                                                                                                                 <  2:241116/150‑1 (MQ=255)
                                       aaTCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCg                                                                                                               >  1:365583/1‑150 (MQ=255)
                                        aTCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAAATGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCgg                                                                                                              >  2:2039430/1‑150 (MQ=255)
                                         tCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCAAGCAGCGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCggc                                                                                                             >  1:2253415/1‑150 (MQ=255)
                                         tCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGAGCTTCGCTCACAACGGCGACCACA‑AC‑GATGCGGTCATTGACCAGTCCCTGGACGCCCAGGGCATACACCCTACCggc                                                                                                             >  1:1740999/1‑150 (MQ=255)
                                          cGATGCTGGTGACCCGTATAGCGGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCggca                                                                                                            <  1:1608305/150‑1 (MQ=255)
                                           gATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCggcag                                                                                                           <  2:2393200/150‑1 (MQ=255)
                                             tGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGTCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCggcaggc                                                                                                         >  2:196198/1‑150 (MQ=255)
                                              gCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGAGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCggcaggca                                                                                                        <  1:2183496/150‑1 (MQ=255)
                                               cTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGgcaggcag                                                                                                       <  2:314606/150‑1 (MQ=255)
                                                tGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACg                                                                                                                                    <  2:1857950/120‑1 (MQ=255)
                                                tGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACg                                                                                                                                    >  1:1857950/1‑120 (MQ=255)
                                                  gTGACCCGTATACCTGCCACCCTAGAACTAGTCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTc                                                                                                    >  1:2275878/1‑150 (MQ=255)
                                                   tGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCa                                                                                                   <  1:2039430/150‑1 (MQ=255)
                                                       ccGTATACCTGCCACCCTCGAGCAAGGTGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGt                                                                                               >  2:467691/1‑150 (MQ=255)
                                                       ccGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGt                                                                                               >  1:1415527/1‑150 (MQ=255)
                                                       ccGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGt                                                                                               <  2:344700/150‑1 (MQ=255)
                                                         gTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCTGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGc                                                                                             >  1:377024/1‑150 (MQ=255)
                                                              ccTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAgg                                                                                                                              >  1:2001953/1‑112 (MQ=255)
                                                              ccTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAgg                                                                                                                              <  2:2001953/112‑1 (MQ=255)
                                                               cTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCAc                                                                                       >  1:964583/1‑150 (MQ=255)
                                                               cTGCCACCCTCGAGCAAGGCGGCGAACTGATGTATCTGGCGCGGACTAATATCTGCTCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCAc                                                                                       <  2:62695/150‑1 (MQ=255)
                                                                tGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑c                                                                                                                                                 <  1:1433632/91‑1 (MQ=255)
                                                                   cACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCggcagg                                                                                                          >  1:1943382/1‑127 (MQ=255)
                                                                   cACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCggcagg                                                                                                          <  2:1943382/127‑1 (MQ=255)
                                                                   cACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCCTGCACGctg                                                                                   <  2:365583/150‑3 (MQ=255)
                                                                   cACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCa                                                                                   >  1:447389/1‑150 (MQ=255)
                                                                      ccTCGAGCAAGGCGGCGACCTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTa                                                                                <  2:1946254/150‑1 (MQ=255)
                                                                      ccTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTa                                                                                <  2:1415527/150‑1 (MQ=255)
                                                                       cTCGAGCAATGCGGCGAACTGCTCTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCATGCCCTGGACGTCAAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAc                                                                               <  1:505430/150‑1 (MQ=255)
                                                                       cTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTTCAGGCACGCCATTAc                                                                               <  1:1147876/150‑1 (MQ=255)
                                                                       cTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAc                                                                               >  1:1794979/1‑150 (MQ=255)
                                                                        tCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAcc                                                                              >  1:824877/1‑150 (MQ=255)
                                                                                ggcggcGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTa                                                                                <  1:1174163/140‑1 (MQ=255)
                                                                                ggcggcGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTa                                                                                >  2:1174163/1‑140 (MQ=255)
                                                                                      gAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCggca                                                                                                            >  2:887393/1‑106 (MQ=255)
                                                                                      gAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCa                                                                                   >  1:123821/1‑131 (MQ=255)
                                                                                      gAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCa                                                                                   <  2:123821/131‑1 (MQ=255)
                                                                                      gAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGGCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGc                                                                >  1:1975182/1‑150 (MQ=255)
                                                                                      gAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCAGTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCggca                                                                                                            <  1:887393/106‑1 (MQ=255)
                                                                                       aaCTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCg                                                               <  2:110256/150‑1 (MQ=255)
                                                                                         ctgctgTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCATGCCCTGGGCGTCCAGGGcctgg                                                                                                                        <  1:1097494/91‑5 (MQ=38)
                                                                                         ctgctgTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCATGCCCTGGGCGTCCAGGGcctgg                                                                                                                        >  2:1097494/1‑87 (MQ=38)
                                                                                          tgctgTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGcc                                                                      >  2:940304/1‑140 (MQ=255)
                                                                                          tgctgTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGcc                                                                      <  1:940304/140‑1 (MQ=255)
                                                                                           gctgTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTAt                                                           <  1:986031/150‑1 (MQ=255)
                                                                                            ctgTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGAc                                                                                                                                                       >  1:1861364/1‑59 (MQ=255)
                                                                                            ctgTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATc                                                          >  1:692688/1‑150 (MQ=255)
                                                                                               tATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAAc                                                       >  2:2314360/1‑150 (MQ=255)
                                                                                                 tCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑CA‑AAATGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAg                                                     >  2:919083/1‑150 (MQ=255)
                                                                                                     gcgcTGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCACACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGCCCGGTGGCGGTATCAACAGCCCg                                                 <  2:1185257/150‑1 (MQ=255)
                                                                                                       gcgGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGc                                               >  1:499189/1‑150 (MQ=255)
                                                                                                                    ttCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGt                                                                                                     >  2:2092887/1‑83 (MQ=255)
                                                                                                                    ttCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTgcgc                                  >  2:140974/1‑150 (MQ=255)
                                                                                                                      cGCTCACAACGAAGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGctc                                <  1:1695126/150‑1 (MQ=255)
                                                                                                                         tCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACTGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGc                             >  1:244935/1‑150 (MQ=255)
                                                                                                                         tCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGc                             >  2:1612103/1‑150 (MQ=255)
                                                                                                                         tCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGc                             >  2:1544448/1‑150 (MQ=255)
                                                                                                                             aaCGGCGAACGCATACAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCCGTGTGCAGGCACGCCATTACCCCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCcgcg                         >  2:1145108/1‑150 (MQ=255)
                                                                                                                              aCGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCTGCCCTGCTGCTGCGCTCTGCCgcgc                        <  1:566325/150‑1 (MQ=255)
                                                                                                                                ggCGAGCGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGccc                      <  2:448543/150‑1 (MQ=255)
                                                                                                                                  cGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTAc                                                                                                                 >  2:2340846/1‑57 (MQ=255)
                                                                                                                                  cGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTAc                                                                                                                 <  1:2340846/57‑1 (MQ=255)
                                                                                                                                    aaCGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCa                                                                                            <  2:2370258/76‑1 (MQ=255)
                                                                                                                                    aaCGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCa                                                                                            >  1:2370258/1‑76 (MQ=255)
                                                                                                                                       gCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGc               >  1:794631/1‑150 (MQ=255)
                                                                                                                                         aTTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCAtg                                                                                                                         <  2:1652866/42‑2 (MQ=255)
                                                                                                                                         aTTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCAtg                                                                                                                         >  1:1652866/1‑41 (MQ=255)
                                                                                                                                          ttCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATg            >  2:1807493/1‑150 (MQ=255)
                                                                                                                                          tccagggCCTGGCA‑TG‑TCAGGCTCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGAc                   >  2:1363784/7‑143 (MQ=255)
                                                                                                                                          tccAGGGCCTGGCA‑TG‑TCAGGCTCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGAc                   <  1:1363784/141‑1 (MQ=255)
                                                                                                                                            cagggCCTGGCA‑TG‑TCAGGCTCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAcctc                                                                            <  1:1187403/80‑1 (MQ=38)
                                                                                                                                            cAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTTGCCGTTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGTGCCCTACCCTCATGGt          <  2:938221/150‑1 (MQ=255)
                                                                                                                                             agggCCTGGCA‑TG‑TCAGGCTCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTg                                                                  >  1:2061535/4‑93 (MQ=255)
                                                                                                                                                 ccTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACTGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCt     <  2:244935/150‑1 (MQ=255)
                                                                                                                                                 ccTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCt     >  1:1910005/1‑150 (MQ=255)
                                                                                                                                                 ccTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGAACGGTCGCCt     >  1:701277/1‑150 (MQ=255)
                                                                                                                                                 ccTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGAAACCTGTGGTCGCCt     >  1:2204815/1‑150 (MQ=255)
                                                                                                                                                  cTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTTCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTg    >  2:79745/1‑150 (MQ=255)
                                                                                                                                                  cTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTg    <  2:229219/150‑1 (MQ=255)
                                                                                                                                                   tGACA‑TG‑CGAGGCCATGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTAAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGCTATCAACCGCCCGGCCCTACTGCTGCGCTCTGCCGCGCCCGACCCGCCTGGGCGCCTgt   >  2:1112784/1‑149 (MQ=255)
                                                                                                                                                   tGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTgg   >  1:754600/1‑150 (MQ=255)
                                                                                                                                                    gACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTggg  >  1:845986/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TGGGCTACTGCGGCATGGGCTGCCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGG  >  NC_002947/5844003‑5844300

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