Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,844,157 C→T 6.1% C241C (TGC→TGT)  PP_5122 → GMC family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,844,1570CT 6.1% 330.9 / 8.0 114C241C (TGC→TGT) PP_5122GMC family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (51/55);  new base T (4/3);  total (56/58)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 118 total aligned reads displayed.

TGCGGCATGGGCTGCCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGC  >  NC_002947/5844011‑5844305
                                                                                                                                                    |                                                                                                                                                    
tGCGGCATGGGCTGCCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAg                                                                                                                                                   <  2:86390/150‑1 (MQ=255)
  cGGCATGGGCTGCCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAAATGCTGTATCTGGCGCGGGCCGGGCGCTTCGCTCACAACGGCGAATGCATTCAGGGCCTTGCA‑TG‑CCAGGc                                                                                                                                                 >  1:938221/1‑150 (MQ=255)
      aTGGGCTGCTAGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTg                                                                                                                                             >  1:1760574/1‑150 (MQ=255)
        gggCTGCCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGa                                                                                                                                           <  1:2451273/150‑1 (MQ=255)
          gCTGCCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCa                                                                                                                                                    <  1:994203/139‑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:194073/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:1260414/150‑1 (MQ=255)
                     aaTGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATa                                                                                                                              <  1:1260652/150‑1 (MQ=255)
                       tGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACCAATACCGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGGCGTCCAGGGcctggc                                                                                                                            <  1:1637366/150‑6 (MQ=255)
                         ccAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACAcc                                                                                                                          >  1:229219/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)
                                   gATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCggcag                                                                                                                <  2:2393200/150‑1 (MQ=255)
                                   gATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCggcag                                                                                                                <  2:2393458/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                                                                                                                                         >  1:1857950/1‑120 (MQ=255)
                                        tGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACg                                                                                                                                         <  2:1857950/120‑1 (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)
                                               ccGTATACCTGCCACCCTCGAGCAAGATCAACAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTAACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGt                                                                                                    >  1:62695/1‑150 (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)
                                                       cTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCAc                                                                                            <  1:320811/150‑1 (MQ=255)
                                                       cTGCCACCCTCGAGCAAGGCGGCGAACTGATGTATCTGGCGCGGACTAATATCTGCTCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCAc                                                                                            <  2:62695/150‑1 (MQ=255)
                                                        tGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑c                                                                                                                                                      >  2:1433632/1‑91 (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)
                                                            aCCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCTGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCCTGCACGctga                                                                                       <  2:377024/150‑4 (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                                                                                     >  2:1174163/1‑140 (MQ=255)
                                                                        ggcggcGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTa                                                                                     <  1:1174163/140‑1 (MQ=255)
                                                                              gAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCggca                                                                                                                 >  2:887393/1‑106 (MQ=255)
                                                                              gAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCa                                                                                        <  2:123821/131‑1 (MQ=255)
                                                                              gAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCa                                                                                        >  1:123821/1‑131 (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)
                                                                                   gctgTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTAt                                                                <  1:986031/150‑1 (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)
                                                                                                    ccGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCctgc                                               <  1:708479/150‑1 (MQ=255)
                                                                                                            ttCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGt                                                                                                          <  1:2092887/83‑1 (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)
                                                                                                                      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)
                                                                                                                            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‑CCAGGCCCTGGACGCCCAGGGCATACACCCGACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATg                 <  1:1568507/150‑1 (MQ=255)
                                                                                                                                  tccAGGGCCTGGCA‑TG‑TCAGGCTCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGAc                        <  1:1363784/141‑1 (MQ=255)
                                                                                                                                    cagggCCTGGCA‑TG‑TCAGGCTCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAcctc                                                                                 <  1:1187403/80‑1 (MQ=38)
                                                                                                                                    cAGGGCCTGGCA‑TG‑TCAGGCTCTAGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAcctc                                                                                 >  2:1187403/1‑84 (MQ=37)
                                                                                                                                    cAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTTGCCGTTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGTGCCCTACCCTCATGGt               <  2:938221/150‑1 (MQ=255)
                                                                                                                                     agggCCTGGCA‑TG‑TCAGGCTCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTg                                                                       >  1:2061535/4‑93 (MQ=255)
                                                                                                                                     agggCCTGGCA‑TG‑TCAGGCTCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTg                                                                       <  2:2061535/90‑1 (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)
                                                                                                                                                    ccAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAAc                                                            <  1:2025149/91‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCGGCATGGGCTGCCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACA‑TG‑CCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGC  >  NC_002947/5844011‑5844305

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