Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,634,508 T→G 19.9% intergenic (+33/+350) putA → / ← PP_5718 proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,634,5080TG19.9% 161.6 / 45.7 94intergenic (+33/+350)putA/PP_5718proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (37/38);  new base G (10/9);  total (47/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-01

CACGGGCTGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGGTCCTT  >  NC_002947/5634359‑5634657
                                                                                                                                                     |                                                                                                                                                     
cACGGGCTGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCt                                                                                                                                                       <  2:679878/150‑1 (MQ=255)
cACGGGCTGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCt                                                                                                                                                       <  2:681985/150‑1 (MQ=255)
    ggCTGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGc                                                                                                                                                   >  1:2153617/1‑150 (MQ=255)
      cTGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGCCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGctc                                                                                                                                                 <  1:330495/150‑1 (MQ=255)
         tCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTTATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGCTCTcc                                                                                                                                              <  2:1788428/150‑1 (MQ=255)
           cAGCGGGGATCACCAGATTGCGCTGGAGCGCCCGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtt                                                                                                                                            >  1:975026/1‑150 (MQ=255)
            aGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCttt                                                                                                                                           <  1:6052/150‑1 (MQ=255)
            aGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGCTCTCCttc                                                                                                                                           <  1:1397331/150‑2 (MQ=255)
             gCGGGGATCACCCGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCACTGGCTAATGGCTGAGGGTGGTTGAAGAAGAAGGATAGCCTAGGCGCTGCtttt                                                                                                                                          >  1:296373/1‑150 (MQ=255)
              cGGGGATCACCAGCTTGCGATGGAGCGACTGGTCATCGATCGCGCCGTGAGCGTGAACAATGCAGCTCATGCATTGAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGCTCTCcttctt                                                                                                                                         <  1:414176/150‑1 (MQ=255)
              cGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAGAGGACGAAGGGATCGGAAGAGTGAGGAGAGCCGCGGCTCTCCttttt                                                                                                                                         >  2:1840749/1‑150 (MQ=255)
               ggggATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtttttt                                                                                                                                        <  2:2153617/150‑1 (MQ=255)
                     cACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGt                                                                                                                                  <  2:1178299/150‑1 (MQ=255)
                      aCCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTg                                                                                                                                 <  1:953005/150‑1 (MQ=255)
                        cAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCa                                                                                                                               >  2:1601056/1‑150 (MQ=255)
                        cAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCa                                                                                                                               >  1:419290/1‑150 (MQ=255)
                           aTTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGt                                                                                                                            <  1:465817/150‑1 (MQ=255)
                                gcTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGCTCTcc                                                                                                                                              <  1:1565338/127‑1 (MQ=255)
                                gcTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGCTCTcc                                                                                                                                              >  2:1565338/1‑127 (MQ=255)
                                      gcgcCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCCCCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAgg                                                                                                                 <  1:1840749/150‑1 (MQ=255)
                                         ccTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGc                                                                                                              <  1:782267/150‑1 (MQ=255)
                                             gTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGa                                                                                                          <  2:1609919/150‑1 (MQ=255)
                                             gTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCATCATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGa                                                                                                          <  1:626521/150‑1 (MQ=255)
                                             gTCATCGAGCGCGCGGCGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGa                                                                                                          <  2:1220085/150‑1 (MQ=255)
                                              tCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGAc                                                                                                         >  1:1063437/1‑150 (MQ=255)
                                                aTCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTCATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGACGCAGGGGCTGGACtc                                                                                                       >  1:826574/1‑149 (MQ=255)
                                                      cgcgcgGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGt                                                                                                 <  2:2303212/150‑1 (MQ=255)
                                                       gcgcgGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTg                                                                                                >  1:1705310/1‑150 (MQ=255)
                                                        cgcgGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTATCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGc                                                                                               <  1:52723/150‑1 (MQ=255)
                                                           ggTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGt                                                                                                                            >  2:1629273/1‑118 (MQ=255)
                                                           ggTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGt                                                                                                                            >  2:1629291/1‑118 (MQ=255)
                                                           ggTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGt                                                                                                                            <  1:1629291/118‑1 (MQ=255)
                                                           ggTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGt                                                                                                                            <  1:1629273/118‑1 (MQ=255)
                                                           ggTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCa                                                                                            <  2:2032076/150‑1 (MQ=255)
                                                            gTGAGCGTGAACACTGCAGCGGCGGGCGGTAACGCCAGCCTGATGACCATTGGCGAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTCTATGCGTGCATGTCCTGAAGCAGAGGCTGGACTAACAGATGCCCAc                                                                                           >  2:1399204/1‑150 (MQ=255)
                                                              gAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTAATGAGTGCATGTTCTGAAGCATGGGCTTGTCTAATAGGTGCTCAccc                                                                                         >  2:1204150/1‑149 (MQ=255)
                                                                 cGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGaa                                                                                      <  2:766053/150‑1 (MQ=255)
                                                                        acTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTc                                                                               >  2:1110018/1‑150 (MQ=255)
                                                                          tGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCAt                                                                             >  1:2393040/1‑150 (MQ=255)
                                                                          tGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCAt                                                                             <  2:1351900/150‑1 (MQ=255)
                                                                            cAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGATCAACAGAACGTAGTCATGc                                                                           <  2:524681/150‑1 (MQ=255)
                                                                             aGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGCTCTCCttc                                                                                                                                           <  1:1343428/85‑2 (MQ=39)
                                                                             aGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGCTCTCCttc                                                                                                                                           >  2:1343428/1‑84 (MQ=39)
                                                                                 gcggGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACg                                                                      <  2:1472449/150‑1 (MQ=255)
                                                                                        gtAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGATGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTACTCACCAGAACGTGGTCATGCCAACGGTGACtc                                                               >  2:785613/1‑149 (MQ=255)
                                                                                         tAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTg                                                              >  2:1487884/1‑150 (MQ=255)
                                                                                          aaCGCCAGCCGGATGGCCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTGATGCGTGCATGTCCTTAAGCAGGGGCTGGTCTAACAGGTGCGCACCAGAACGTGGTCATGCCAAAGGTGaagtgc                                                             <  2:975026/150‑6 (MQ=255)
                                                                                           aCGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAgg                                                                                                                 >  2:1630551/1‑97 (MQ=255)
                                                                                           aCGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAgg                                                                                                                 <  1:1630551/97‑1 (MQ=255)
                                                                                            ggCCAGCCTGATGACCATGGGGGAATGGCTGGGGGTTGTTGAAGAAGAAGGAGAGCCGCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgt                                                           <  2:826574/149‑1 (MQ=255)
                                                                                             gCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgtg                                                          >  2:2387570/1‑150 (MQ=255)
                                                                                             gCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgtg                                                          >  1:1173187/1‑150 (MQ=255)
                                                                                                 gCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCAt                                                                                                                              <  2:1145557/78‑1 (MQ=255)
                                                                                                 gCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCAt                                                                                                                              >  1:1145557/1‑78 (MQ=255)
                                                                                                 gCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGtctt                                                      <  1:817630/150‑1 (MQ=255)
                                                                                                  ccTGATGACCATTTTCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGtcttc                                                     <  2:1357674/150‑1 (MQ=255)
                                                                                                   cTGATGACCATTGGCTAATGGCTGAGGGTTGTTGTAAAACAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGtcttct                                                    <  2:2017769/150‑1 (MQ=255)
                                                                                                    tgatgaCCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCtg                                                   <  2:419290/150‑1 (MQ=255)
                                                                                                      atgaCCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGa                                                                  <  1:2398507/133‑1 (MQ=255)
                                                                                                      atgaCCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGa                                                                  >  2:2398507/1‑133 (MQ=255)
                                                                                                      atgaCCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCtgtt                                                 <  2:1049432/150‑1 (MQ=255)
                                                                                                      atgaCCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCCGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCtgtt                                                 >  2:999269/1‑150 (MQ=255)
                                                                                                          ccATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAg                                                                                                                  <  2:1203672/81‑1 (MQ=255)
                                                                                                          ccATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAg                                                                                                                  >  1:1203672/1‑81 (MQ=255)
                                                                                                          ccATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTgg                                                             >  1:520799/1‑134 (MQ=255)
                                                                                                          ccATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTgg                                                             <  2:520799/134‑1 (MQ=255)
                                                                                                           cATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTAc                                            >  1:1582941/1‑150 (MQ=255)
                                                                                                           cATTGGCTAATGGCTGAGGGGTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTAc                                            >  1:1734559/1‑150 (MQ=255)
                                                                                                             ttGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTAcca                                          >  1:2357167/1‑149 (MQ=255)
                                                                                                               ggCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGtt                                        <  2:1582941/150‑1 (MQ=255)
                                                                                                                    aTGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgt                                                           >  2:1247195/1‑126 (MQ=255)
                                                                                                                    aTGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgt                                                           <  1:1247195/126‑1 (MQ=255)
                                                                                                                                  ttgAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCg                     <  1:1677317/150‑1 (MQ=255)
                                                                                                                                     ataaaaaagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCt                  >  2:2441589/7‑150 (MQ=255)
                                                                                                                                     aagaagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCt                  >  1:2433381/1‑150 (MQ=255)
                                                                                                                                     aagaagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCt                  >  1:1506938/1‑150 (MQ=255)
                                                                                                                                       gaagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACg                                                                                    <  1:2314179/82‑1 (MQ=255)
                                                                                                                                       gaagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACg                                                                                    <  1:1050032/82‑1 (MQ=255)
                                                                                                                                       gaagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACg                                                                                    >  2:1050032/1‑82 (MQ=255)
                                                                                                                                       gaagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACg                                                                                    >  2:2314179/1‑82 (MQ=255)
                                                                                                                                       gaagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGc                <  1:1110018/150‑1 (MQ=255)
                                                                                                                                         aaaagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCt              >  2:49237/3‑150 (MQ=255)
                                                                                                                                         aaaagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCt              >  2:803617/3‑150 (MQ=255)
                                                                                                                                          gaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTg             >  2:528767/1‑150 (MQ=255)
                                                                                                                                          gaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTg             >  2:652524/1‑150 (MQ=255)
                                                                                                                                          aaagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgt                                                           >  2:1254395/2‑104 (MQ=255)
                                                                                                                                          aaagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgt                                                           <  1:1254395/103‑1 (MQ=255)
                                                                                                                                           aagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTg                                                                                                            >  1:440977/1‑54 (MQ=37)
                                                                                                                                           aagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTg                                                                                                            <  2:440977/54‑1 (MQ=37)
                                                                                                                                            agGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTggg           >  1:1358981/1‑150 (MQ=255)
                                                                                                                                            agGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTggg           >  2:614076/1‑150 (MQ=255)
                                                                                                                                             ggAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGc          >  2:933490/1‑150 (MQ=255)
                                                                                                                                              gagaGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGcc         >  1:1653665/1‑150 (MQ=255)
                                                                                                                                                gagCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCCGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCgg       <  1:999269/150‑1 (MQ=255)
                                                                                                                                                   ccTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGGTcc    >  1:2224919/1‑150 (MQ=255)
                                                                                                                                                    cTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGGTCCt   <  2:1653665/150‑1 (MQ=255)
                                                                                                                                                     tCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGGTCCtt  >  1:1053025/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACGGGCTGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGGTCCTT  >  NC_002947/5634359‑5634657

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