Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,102,762 C→A 5.3% R67R (CGG→AGG)  ttg2E → toluene‑tolerance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,102,7620CA 5.3% 242.4 / 2.4 76R67R (CGG→AGG) ttg2Etoluene‑tolerance protein
Reads supporting (aligned to +/- strand):  ref base C (32/39);  new base A (2/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

TGGCGTGCTGGACTACCGCAGCGGCCCGGCCTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGCAAGCTAGTGGCCAGGCACGGATGATGGCCCCTCGGTCAGCGAAGCCTTC  >  NC_002947/1102614‑1102908
                                                                                                                                                    |                                                                                                                                                  
tGGCGTGCTGGACTACCGCAGCGGCCCGGCCTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGATATTTTCGTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCg                                                                                                                                                   >  1:404200/1‑150 (MQ=255)
 ggCGTGCTGGACTACCGCAGCGGCCCGGCCTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCgg                                                                                                                                                  <  1:1000765/150‑1 (MQ=255)
     tGCTGGACTACCGCAGCGGCCCGGCCTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTTGCGTACATGCGGGATg                                                                                                                                              >  1:687428/1‑150 (MQ=255)
        tGGACTACCGCAGCGGCCCGGCCTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTc                                                                                                                                           >  1:985060/1‑150 (MQ=255)
        tGGACTACCGCAGCGGCCCGGCCTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTc                                                                                                                                           >  1:39494/1‑150 (MQ=255)
         ggACTACCGCAGCGGCCCGGCCTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATTCGGTATGCTCa                                                                                                                                          <  2:316194/150‑1 (MQ=255)
         ggACTACCGCAGCGGCCCGGCCTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCa                                                                                                                                          <  1:747741/150‑1 (MQ=255)
               ccGCAGCGGCCCGGCCTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTTGCGTACATGCGGGATGCTCAGGCCAc                                                                                                                                    <  2:687428/150‑1 (MQ=255)
               ccGCAGCGGCCCGGCCTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCAc                                                                                                                                    >  1:575148/1‑150 (MQ=255)
               ccGCAGCGGCCCGGCCTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCAc                                                                                                                                    >  1:629349/1‑150 (MQ=255)
                 gCAGCGGCCCGGCCTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTg                                                                                                                                  <  2:647969/150‑1 (MQ=255)
                       gcccggccTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCCgg                                                                                                                                        <  1:1210709/138‑1 (MQ=255)
                       gcccggccTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAgg                                                                                                                                        >  2:1210709/1‑138 (MQ=255)
                       gcccggccTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAgg                                                                                                                            <  2:629349/150‑1 (MQ=255)
                        cccggccTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGc                                                                                                                           >  1:772142/1‑150 (MQ=255)
                         ccggccTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCa                                                                                                                                     >  2:620223/1‑139 (MQ=255)
                         ccggccTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCa                                                                                                                                     <  1:620223/139‑1 (MQ=255)
                         ccggccTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCg                                                                                                                          >  1:722851/1‑150 (MQ=255)
                          cggccTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGt                                                                                                                         <  2:39494/150‑1 (MQ=255)
                              ctggcgcAAGCAAGGCAAGGTGCTGGTCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTGGGCCTGTCGCGACTGCTGGAGTACATGCGGGATGCTCAGGTCACTGGCAAGGCGTGCGa                                                                                                                     <  2:434890/147‑1 (MQ=255)
                                   gcaagcaaGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGc                                                                                                                <  2:183016/150‑1 (MQ=255)
                                     aagcaagGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGt                                                                                                              >  2:819069/1‑150 (MQ=255)
                                     aagcaagGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGt                                                                                                              <  1:317705/150‑1 (MQ=255)
                                     aagcaagGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGt                                                                                                              >  2:168265/1‑150 (MQ=255)
                                      agcaagGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTg                                                                                                             <  2:200440/150‑1 (MQ=255)
                                       gcaaggcaagGTGCTGAACGATGGAACGCGCGAGCCGCGAATGGTGATTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCATGTCGCTACTGCTGGCATACATGCGGGATGCCCAGACCACTGGCAAGGCGTGCGAGGTGCGTgg                                                                                                            >  2:1227042/1‑150 (MQ=255)
                                          aggcaaggTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCAt                                                                                                         <  2:480189/150‑1 (MQ=255)
                                             caaggTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGcc                                                                                                      <  1:992146/150‑1 (MQ=255)
                                              aaggTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCt                                                                                                     >  2:407852/1‑150 (MQ=255)
                                                 gTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCtgat                                                                                                  <  1:588953/150‑1 (MQ=255)
                                                    cTGATCGCGCGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTgatgat                                                                                               <  1:695662/150‑1 (MQ=255)
                                                     tGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGata                                                                                              >  2:257791/1‑150 (MQ=255)
                                                     tGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGata                                                                                              >  1:200667/1‑150 (MQ=255)
                                                      gATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCa                                                                                                                                          >  1:758963/1‑105 (MQ=255)
                                                      gATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCa                                                                                                                                          <  2:758963/105‑1 (MQ=255)
                                                       aTCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGa                                                                                                                     >  1:1193537/1‑125 (MQ=255)
                                                       aTCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGa                                                                                                                     <  2:1193537/125‑1 (MQ=255)
                                                        tCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGc                                                                                           <  2:722851/150‑1 (MQ=255)
                                                        tCGATGGCTCCCGAGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGc                                                                                           <  1:1227042/150‑1 (MQ=255)
                                                         cGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGGAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCg                                                                                          <  1:716167/150‑1 (MQ=255)
                                                               cTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAAt                                                                                    <  1:712674/150‑1 (MQ=255)
                                                                tCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGGCGAGTGTTGGCCTGTCGCTACGCCTGGCGAACATCTGTTCTGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATACGCGGGAAAtt                                                                                   >  1:91052/1‑150 (MQ=255)
                                                                   cgcgAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGcc                                                                                <  1:819069/150‑1 (MQ=255)
                                                                      gAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAg                                                                             >  1:180410/1‑150 (MQ=255)
                                                                           gcgACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATTATATGCGGGAAATTGCCGAGGTCTa                                                                        <  1:60021/150‑1 (MQ=255)
                                                                                      ttGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATg                                                             >  2:41425/1‑150 (MQ=255)
                                                                                        ggttGCTCGGCGGTGGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGaa                                                           <  2:1036073/148‑1 (MQ=255)
                                                                                             cTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCt                                                      <  2:985060/150‑1 (MQ=255)
                                                                                             cTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCt                                                      <  2:575148/150‑1 (MQ=255)
                                                                                              tCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTg                                                     <  1:407852/150‑1 (MQ=255)
                                                                                                  cTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGCaa                                                 <  1:168265/150‑1 (MQ=255)
                                                                                                  cTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGCaa                                                 <  1:257791/150‑1 (MQ=255)
                                                                                                        agaAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCAACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGCAAGCTAGt                                           >  2:80002/1‑150 (MQ=255)
                                                                                                               gtcgAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGCaa                                                 <  2:399489/137‑1 (MQ=255)
                                                                                                               gtcgAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGCaa                                                 >  1:399489/1‑137 (MQ=255)
                                                                                                                tcgAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGc                                                   <  2:1281362/134‑1 (MQ=255)
                                                                                                                tcgAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGc                                                   >  1:1281362/1‑134 (MQ=255)
                                                                                                                tcgAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGCAAGCTAGTGGCCAGGc                                   >  2:904861/1‑150 (MQ=255)
                                                                                                                       ttGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTa                                                                        >  2:133473/1‑106 (MQ=255)
                                                                                                                       ttGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTa                                                                        <  1:133473/106‑1 (MQ=255)
                                                                                                                         ggCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGCAAGCTAGTGGCCAGGCACGgatgat                          >  1:524940/1‑150 (MQ=255)
                                                                                                                          gCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGc                                                                                                                                >  2:1185656/1‑47 (MQ=255)
                                                                                                                          gCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGc                                                                                                                                <  1:1185656/47‑1 (MQ=255)
                                                                                                                          gCCTGTCGCTACTGCTGGCGTACATGAGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGAt                                                                    <  1:705910/107‑1 (MQ=255)
                                                                                                                          gCCTGTCGCTACTGCTGGCGTACATGAGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGAt                                                                    >  2:705910/1‑107 (MQ=255)
                                                                                                                          gCCTGTCGCTACTGCTGGCGTACATGAGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGAt                                                                    >  2:706047/1‑107 (MQ=255)
                                                                                                                          gCCTGTCGCTACTGCTGGCGTACATGAGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGAt                                                                    <  1:706047/107‑1 (MQ=255)
                                                                                                                              gTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGCAAGCTAGTGGCCAGGCACGGATGATGGccc                     <  2:404200/150‑1 (MQ=255)
                                                                                                                                 gCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGCAAGCTAGTGGCCAGGCACGGATGATGGCCCCTc                  >  2:288280/1‑150 (MQ=255)
                                                                                                                                  cTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGCAAGCTAGTGGCCAGGCACGGATGATGGCCCCTCg                 <  2:772142/150‑1 (MQ=255)
                                                                                                                                      tgctgGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGAATGATGATATGAGGGCCATCAAGGCCAGCTAGCAGTTCGATGAAGTGCTGGCAAGCTAGTGGCCAGGCACGGATGATGTCCCCTCGGTCa             >  1:1254498/1‑150 (MQ=255)
                                                                                                                                         tgGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGCAAGCTAGTGGCCAGGCACGGATGATGGCCCCTCGGTCAGCg          >  1:687083/1‑150 (MQ=255)
                                                                                                                                           gCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGCAAGCTAGTGGCCAGGCACGGATGATGGCCCCTCGGTCAGCGaa        <  1:904861/150‑1 (MQ=255)
                                                                                                                                                 aTGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGCAAGCTAGTGGCCAGGCACGGATGATGGCCCCTCGGTCAGCGAAGCCTTc  >  2:673764/1‑150 (MQ=255)
                                                                                                                                                  tGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGAtctc                                                                 <  2:556596/86‑1 (MQ=255)
                                                                                                                                                  tGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGAtctc                                                                 >  1:556596/1‑86 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGGCGTGCTGGACTACCGCAGCGGCCCGGCCTTGCGCAAGCAAGGCAAGGTGCTGATCGATGGCTCGCGCGAGCCGCGACTGGTGCTTGATTGCTCGGCTGTCGAGAAGTCGTCGAGTGTTGGCCTGTCGCTACTGCTGGCGTACATGCGGGATGCTCAGGCCACTGGCAAGGCGTGCGAGGTGCGTGGCATGCCTGATGATATGCGGGAAATTGCCGAGGTCTATGATCTCGATGAAGTGCTGGCAAGCTAGTGGCCAGGCACGGATGATGGCCCCTCGGTCAGCGAAGCCTTC  >  NC_002947/1102614‑1102908

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: