Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,634,508 T→G 11.1% 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,5080TG11.1% 239.0 / 21.1 99intergenic (+33/+350)putA/PP_5718proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (46/40);  new base G (7/4);  total (54/45)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

CTGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGC  >  NC_002947/5634365‑5634649
                                                                                                                                               |                                                                                                                                             
cTGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGctc                                                                                                                                         >  2:1729717/1‑150 (MQ=255)
   tCCAGCGGGGATCACCAGATTGCGCTGTAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAAGGGCTGAGTGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTcc                                                                                                                                      <  1:1522543/150‑1 (MQ=255)
   tCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTCGTAACCGTGCAATGCACGCGTACGAACACTGAACGGGTTTGTGGTAACGCCAGCCTGATGACCATTGGCTAAGGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTcc                                                                                                                                      >  1:36273/1‑150 (MQ=255)
     cAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtt                                                                                                                                    >  1:642363/1‑150 (MQ=255)
      aGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCttt                                                                                                                                   >  1:1783820/1‑150 (MQ=255)
       gCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtttt                                                                                                                                  <  1:1729717/150‑1 (MQ=255)
       gCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtttt                                                                                                                                  >  1:1387971/1‑150 (MQ=255)
       gCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGACTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtttt                                                                                                                                  >  1:1876675/1‑150 (MQ=255)
           ggATCACCAAATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTAt                                                                                                                              >  2:775234/1‑150 (MQ=255)
                aCCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTg                                                                                                                         >  2:1344919/1‑150 (MQ=255)
                aCCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTg                                                                                                                         <  2:409140/150‑1 (MQ=255)
                 ccAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGc                                                                                                                        >  2:345742/1‑150 (MQ=255)
                 ccAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGc                                                                                                                        >  1:1308838/1‑150 (MQ=255)
                 ccAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGACCCTCGGCTCTCCTTTTTTATGCGTGc                                                                                                                        <  1:345742/150‑1 (MQ=255)
                   aGATTGCGCTGGTGCGCCTGGTCATCGAGCGAGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCAt                                                                                                                      <  2:36273/150‑1 (MQ=255)
                    gATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGGGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGGTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATg                                                                                                                     >  2:1239636/1‑150 (MQ=255)
                     aTTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGt                                                                                                                    <  1:775234/150‑1 (MQ=255)
                         cgcTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTg                                                                                                                <  2:1532833/150‑1 (MQ=255)
                           cTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGaa                                                                                                              >  1:358342/1‑150 (MQ=255)
                                gcgcCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAgg                                                                                                         <  2:1387971/150‑1 (MQ=255)
                                     tGGTCATCGAGCGCGCGGTGAGCGTGAGCACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGAGTGAGGGTTGTTGAGGATGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTg                                                                                                    <  2:197344/150‑1 (MQ=255)
                                     tGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTg                                                                                                    <  2:479429/150‑1 (MQ=255)
                                       gTCATCGAGCGCGCGGTGAGCGTGAACACTGCCGCGGTGGGTAGACGCGCCTGTTTGATTATCATTTTCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGCGCCACGGCTCTCCTTCTTTATGAATGCATGTCCTGAAGCAGGGGCTGGa                                                                                                  <  2:1876675/150‑1 (MQ=255)
                                       gTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGa                                                                                                  <  1:1420405/150‑1 (MQ=255)
                                          aTCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTa                                                                                               >  2:178718/1‑150 (MQ=255)
                                           tCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGATGTTGAAGAAGAAGGAGAGCCCAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTaa                                                                                              >  1:128356/1‑150 (MQ=255)
                                             gAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGa                                                                                                  >  2:13980/1‑144 (MQ=255)
                                             gAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGa                                                                                                  <  1:13980/144‑1 (MQ=255)
                                              tgtgcgcgcTGAGCGTGAACCCTTCAGCGGCGGGTGGTAATGCCAGCCTGATGAGGATTGGCTAATGGCTGAGGGTTGCTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAg                                                                                           <  2:115064/147‑1 (MQ=255)
                                               gcgcgcGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAgg                                                                                          >  1:1915736/1‑150 (MQ=255)
                                                cgcgcgGCGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGt                                                                                         <  2:1929298/150‑1 (MQ=255)
                                                       tGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCAcc                                                                                  >  1:1170441/1‑150 (MQ=255)
                                                        gAGCGTGAACACTGCAGTGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCa                                                                                 <  1:1469446/150‑1 (MQ=255)
                                                         aGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtttt                                                                                                                                  >  2:693173/1‑100 (MQ=255)
                                                         aGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtttt                                                                                                                                  <  1:693173/100‑1 (MQ=255)
                                                          gCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGa                                                                               >  1:6970/1‑150 (MQ=255)
                                                              gAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGCGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGt                                                                           >  2:1719223/1‑150 (MQ=255)
                                                               aaCACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCGTGTCCTGAAGCAGGGGCTGGACTAAc                                                                                             >  1:604907/1‑131 (MQ=255)
                                                               aaCACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCGTGTCCTGAAGCAGGGGCTGGACTAAc                                                                                             <  2:604907/131‑1 (MQ=255)
                                                               aaCACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTg                                                                          <  2:1783820/150‑1 (MQ=255)
                                                                 cacTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGt                                                                        <  2:1805026/150‑1 (MQ=255)
                                                                    tGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAgg                                                                                                         <  1:447982/114‑1 (MQ=255)
                                                                    tGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAgg                                                                                                         >  2:447982/1‑114 (MQ=255)
                                                                    tGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCAt                                                                     <  2:444648/150‑1 (MQ=255)
                                                                      cAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTg                                                                                                    <  2:817284/117‑1 (MQ=255)
                                                                      cAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTg                                                                                                    >  1:817284/1‑117 (MQ=255)
                                                                           gcggGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGt                                                                                         >  2:1758560/1‑123 (MQ=255)
                                                                           gcggGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGt                                                                                         <  1:1758560/123‑1 (MQ=255)
                                                                           gcggGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACg                                                              <  1:1344919/150‑1 (MQ=255)
                                                                               gtggtAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGa                                                          <  2:1308838/150‑1 (MQ=255)
                                                                                  gtAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACtt                                                       <  1:1719223/150‑1 (MQ=255)
                                                                                    aaCGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTgg                                                     >  1:55022/1‑150 (MQ=255)
                                                                                       gCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGCTCATGCCAACGGTGACTTGGgtg                                                  >  2:195604/1‑150 (MQ=255)
                                                                                         cAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTg                                                                                        >  2:326809/1‑110 (MQ=255)
                                                                                         cAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTg                                                                                        <  1:326809/110‑1 (MQ=255)
                                                                                          aGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTTGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAAc                                                               <  1:280691/134‑1 (MQ=255)
                                                                                          aGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTTGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAAc                                                               >  2:280691/1‑134 (MQ=255)
                                                                                             cTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTgg                                                     <  1:874353/141‑1 (MQ=255)
                                                                                             cTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTgg                                                     >  2:874353/1‑141 (MQ=255)
                                                                                                  gaCCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCtgttgt                                       <  1:1239636/150‑1 (MQ=255)
                                                                                                   aCCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAgg                                                                                                         <  1:1672639/83‑1 (MQ=255)
                                                                                                   aCCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAgg                                                                                                         >  2:1672639/1‑83 (MQ=255)
                                                                                                        tGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGt                                 <  2:1915736/150‑1 (MQ=255)
                                                                                                         ggCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTg                                                      <  2:550134/128‑1 (MQ=255)
                                                                                                         ggCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTg                                                      >  1:550134/1‑128 (MQ=255)
                                                                                                            tAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCNTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTAt                             >  1:1787787/1‑150 (MQ=255)
                                                                                                            tAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTa                                     >  1:1158072/1‑142 (MQ=255)
                                                                                                            tAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTa                                     <  2:1158072/142‑1 (MQ=255)
                                                                                                            tAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTAt                             >  2:1153047/1‑150 (MQ=255)
                                                                                                             aaTGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATc                            <  1:1153047/150‑1 (MQ=255)
                                                                                                                 gCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATcc                           >  2:76800/1‑147 (MQ=255)
                                                                                                                 gCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATcc                           <  1:76800/147‑1 (MQ=255)
                                                                                                                   tGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTCTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTgg                      >  1:1272124/1‑150 (MQ=255)
                                                                                                                       ggTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTCTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAAt                  >  1:1567929/1‑150 (MQ=255)
                                                                                                                           gttgAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCAt                                                                     >  2:1224457/1‑95 (MQ=255)
                                                                                                                           gttgAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCAt                                                                     <  1:1224457/95‑1 (MQ=255)
                                                                                                                            ttgAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCAt                                                                     <  2:424253/94‑1 (MQ=255)
                                                                                                                            ttgAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCAt                                                                     >  1:424253/1‑94 (MQ=255)
                                                                                                                            ttgAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCg             >  2:1541645/1‑150 (MQ=255)
                                                                                                                              cataaaaaagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGc           >  2:737283/8‑150 (MQ=255)
                                                                                                                               aagaagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgtg                                                  <  1:1519371/110‑1 (MQ=255)
                                                                                                                               aagaagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgtg                                                  >  2:1519371/1‑110 (MQ=255)
                                                                                                                                 gaagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGc        >  2:1699063/1‑150 (MQ=255)
                                                                                                                                  aagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGcc       >  1:1142163/1‑150 (MQ=255)
                                                                                                                                   aaaagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGtc                                                >  2:938826/3‑108 (MQ=255)
                                                                                                                                   aaaagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGtc                                                <  1:938826/106‑1 (MQ=255)
                                                                                                                                   aaaagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCt      >  1:22644/3‑150 (MQ=255)
                                                                                                                                    gaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTg     >  1:404475/1‑150 (MQ=255)
                                                                                                                                    aaagGAGAGCCGAGGCTCTCCTTCTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCAAACGGTATCTTGGGAGTCTTCTGTTGTTACCGTTTATCCCATGGCAATTCTCGGCTGCCTg     >  2:1022000/2‑150 (MQ=255)
                                                                                                                                      agGAGATCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTggg   >  2:1781568/1‑150 (MQ=255)
                                                                                                                                      agGAGAGCCTCGGCTCTCCTTTTTTATGTGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTggg   >  1:1643159/1‑150 (MQ=255)
                                                                                                                                      agGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCa                                                                                 <  2:1810034/72‑1 (MQ=39)
                                                                                                                                      agGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCa                                                                                 >  1:1810034/1‑72 (MQ=39)
                                                                                                                                       ggcgagCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCtgtt                                         <  2:1326646/108‑1 (MQ=255)
                                                                                                                                       ggAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCtgtt                                         >  1:1326646/1‑111 (MQ=255)
                                                                                                                                       ggAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGc  <  2:534223/150‑1 (MQ=255)
                                                                                                                                       ggAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGc  >  1:810290/1‑150 (MQ=255)
                                                                                                                                        gagaGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTg     >  2:1314794/1‑146 (MQ=255)
                                                                                                                                        gagaGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTg     <  1:1314794/146‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
CTGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGC  >  NC_002947/5634365‑5634649

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