Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,353,932 C→T 8.8% H338H (CAC→CAT yrpB → 2‑nitropropane dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,353,9320CT 8.8% 189.4 / 8.0 68H338H (CAC→CATyrpB2‑nitropropane dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (27/35);  new base T (4/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

CGCGAGCTGGGGCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCTCAAGGAGCT  >  NC_002947/4353786‑4354080
                                                                                                                                                  |                                                                                                                                                     
cgcgAGCTGGGGCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATg                                                                                                                                                    >  2:194967/1‑150 (MQ=255)
     gCTGGGGCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCATGTGCcggc                                                                                                                                               <  2:289241/150‑1 (MQ=255)
         gggCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGt                                                                                                                                           <  1:194967/150‑1 (MQ=255)
         gggCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGt                                                                                                                                           <  2:299145/150‑1 (MQ=255)
               aTGCGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTg                                                                                                                                          <  2:610936/145‑1 (MQ=255)
               aTGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTg                                                                                                                                          >  1:610936/1‑145 (MQ=255)
                    cgagcTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGTAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCg                                                                                                                                <  1:633963/150‑1 (MQ=255)
                     gagcTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGTTGATACATGCCGGCGGGTGAGCTGACCCgg                                                                                                                               <  1:308682/150‑1 (MQ=255)
                      agcTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCggg                                                                                                                              <  1:723935/150‑1 (MQ=255)
                       gcTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGga                                                                                                                             <  2:822609/150‑1 (MQ=255)
                        cTTGCACCGCGCTTCCGGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGCTCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGCGCTGACCCGGgag                                                                                                                            >  2:775699/1‑150 (MQ=255)
                           gCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGAtt                                                                                                                         >  1:685626/1‑150 (MQ=255)
                            cACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTg                                                                                                                        >  2:1063001/1‑150 (MQ=255)
                             aCCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCATTTACGAGCGATCCCCAATCCGCAGGAAATGGATGTTTAATGGAATCAATCACAGTGGCGGTCTTGGCTGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGc                                                                                                                       >  2:679319/1‑150 (MQ=255)
                                   cTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTaa                                                                                                                 <  2:675834/150‑1 (MQ=255)
                                    ttCCCGCTGGCGGGCGGGGCATTGATGCCTTTCCGGCAGACAACAGATCGAAGGGGTAAGAGTGATTTCAGCGGTCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGGGAACTGACCCGGGAGATTGCCGGTAAg                                                                                                                >  1:466981/1‑150 (MQ=255)
                                    ttCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAg                                                                                                                <  2:685626/150‑1 (MQ=255)
                                    ttCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAg                                                                                                                >  1:355291/1‑150 (MQ=255)
                                        cGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGcct                                                                                                            <  2:715677/150‑3 (MQ=255)
                                           tGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTgg                                                                                                         <  1:433420/150‑1 (MQ=255)
                                            ggcgggcgGGGCATTGATGCCTTTAACGACGGTCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGc                                                                                                        >  1:308302/1‑150 (MQ=255)
                                             gcgggcggGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCa                                                                                                       <  1:173457/150‑1 (MQ=255)
                                                  cgggGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGAt                                                                                                  <  2:503440/150‑1 (MQ=255)
                                                   ggggCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGTCAGTGATc                                                                                                 <  1:705648/150‑1 (MQ=255)
                                                    gggCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCg                                                                                                <  1:606115/150‑1 (MQ=255)
                                                      gCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGt                                                                                              >  2:396908/1‑150 (MQ=255)
                                                        aTTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCgg                                                                                                                               >  1:506451/1‑115 (MQ=255)
                                                        aTTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCgg                                                                                                                               <  2:506451/115‑1 (MQ=255)
                                                           gATGCCTATACGGGCAAGCACCGATCCGCAGGGGAAGAGCGATATCAGCAATCTGTAGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGAGAGCTGACCCGGGAGATTCCCACTAAGGCATTGACAGTGATCgagcagca                                                                                         >  2:496916/1‑143 (MQ=255)
                                                                  ttACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCTGTATGGGCGCTTCGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTc                                                                                  <  2:308302/150‑1 (MQ=255)
                                                                  ttACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGTCCGGCTCATGCCTGTGGGTGCGCAAAGCCAGGAAATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTc                                                                                  <  1:570698/150‑1 (MQ=255)
                                                                  ttACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTc                                                                                  <  1:325210/150‑1 (MQ=255)
                                                                      gggCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGAc                                                                              >  2:452733/1‑150 (MQ=255)
                                                                       ggCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGAca                                                                             <  2:788187/150‑1 (MQ=255)
                                                                              aCCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACaggcaggc                                                                     >  1:396399/1‑150 (MQ=255)
                                                                                 gATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATg                                                                  <  2:548919/150‑1 (MQ=255)
                                                                                 gATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATg                                                                  <  2:355291/150‑1 (MQ=255)
                                                                                  aTCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATgg                                                                 <  2:1111550/150‑1 (MQ=255)
                                                                                            tAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACg                                                       <  1:938107/150‑1 (MQ=255)
                                                                                              agagTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGc                                                     <  2:1103469/150‑1 (MQ=255)
                                                                                                agTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCa                                                   >  1:174017/1‑150 (MQ=255)
                                                                                                agTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCa                                                   >  1:49302/1‑150 (MQ=255)
                                                                                                 gTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCaa                                                  <  1:1063001/150‑1 (MQ=255)
                                                                                                  tGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGAcaggca                                                                         <  2:221547/127‑1 (MQ=255)
                                                                                                  tGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGAcaggca                                                                         >  1:221547/1‑127 (MQ=255)
                                                                                                  tGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAAt                                                 <  1:396908/150‑1 (MQ=255)
                                                                                                  tGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAAt                                                 >  1:210637/1‑150 (MQ=255)
                                                                                                  tGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTACGCCGGCACATGCCGGCGGGTGAGCTGACCCg                                                                                                                                <  2:226781/72‑1 (MQ=255)
                                                                                                  tGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTACGCCGGCACATGCCGGCGGGTGAGCTGACCCg                                                                                                                                >  1:226781/1‑72 (MQ=255)
                                                                                                          gCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGc                                                                                                                       <  2:356667/73‑1 (MQ=255)
                                                                                                          gCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGc                                                                                                                       >  1:356667/1‑73 (MQ=255)
                                                                                                                    gTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGt                               >  2:476599/1‑150 (MQ=255)
                                                                                                                     tCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCa                                                                                                       <  2:799475/78‑1 (MQ=255)
                                                                                                                     tCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCa                                                                                                       >  1:799475/1‑78 (MQ=255)
                                                                                                                     tCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGtt                              >  1:792158/1‑150 (MQ=255)
                                                                                                                        gggCAGGCGCTACGGCTAGGCCGGCACATGCCGCCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCAGTGGAAGTGATCGGTGACCAGGCATTCTCACACGCAGGCATGGGATCTTCTCGGCCAATGGCTTATCGCTCTATAGTTAAg                           <  2:299558/150‑1 (MQ=255)
                                                                                                                         ggCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAgg                          >  1:812975/1‑150 (MQ=255)
                                                                                                                         ggCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAgg                          <  1:452733/150‑1 (MQ=255)
                                                                                                                          gCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGc                         >  1:35549/1‑150 (MQ=255)
                                                                                                                                 ttACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACg                  >  1:895045/1‑150 (MQ=255)
                                                                                                                                   aCGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGAt                <  2:562889/150‑1 (MQ=255)
                                                                                                                                         ccgCCGGCATGTGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGc                         <  2:282308/133‑1 (MQ=255)
                                                                                                                                         ccgCCGGCATGTGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGc                         >  1:282308/3‑135 (MQ=255)
                                                                                                                                         ccgCCGGCATGTGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCta          >  1:224733/3‑149 (MQ=255)
                                                                                                                                         ccgCCGGCATGTGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCTc          >  1:599130/3‑150 (MQ=255)
                                                                                                                                         ccgCCGGCATGTGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGAATAAAACTAGCTGTCCATGGCCTCTGGACGGCAAATGGACTCAAGCTTGCTAGTGAAGGCTATAACGATAACCTc          >  1:570209/3‑150 (MQ=255)
                                                                                                                                              ggCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCTCAAGGa     >  2:230514/1‑150 (MQ=255)
                                                                                                                                              ggCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCTCAAGGa     <  2:895045/150‑1 (MQ=255)
                                                                                                                                                 acaTGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCTCAAGGAGCt  <  2:812975/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGCGAGCTGGGGCCGATGAGCGAGCTTGCACCGCGCTTCCCGCTGGCGGGCGGGGCATTGATGCCTTTACGGGCGATCACCGATCCGCAGGGTAAGAGTGATTTCAGCAATCTGTGGTCGGGGCAGGCGTTACGGCTAGGCCGGCACATGCCGGCGGGTGAGCTGACCCGGGAGATTGCCGGTAAGGCATTGGCAGTGATCGGTCACCAGGCATTCTGACAGGCAGGCATGGCCTCTTCACGGCCAATGGCTTATCGCTATATAGTTAAGGCTATAACGATAACCTCAAGGAGCT  >  NC_002947/4353786‑4354080

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