Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,195,808 A→T 5.7% L101Q (CTG→CAG)  dhaT ← 1,3‑propanediol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,195,8080AT 5.7% 271.1 / 6.2 87L101Q (CTG→CAG) dhaT1,3‑propanediol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (43/39);  new base T (2/3);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

GCCGGCCGTGGTCGGGATGGCGACGATGGGCGGGATGGCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTGGCGTAGACGCTGTCCAGCAGGCCTGCT  >  NC_002947/3195663‑3195954
                                                                                                                                                  |                                                                                                                                                  
gccggccgTGGTCGGGATGGCGACGATGGGCGGGATGGCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCcggcag                                                                                                                                                <  1:7035/150‑1 (MQ=255)
 ccggccgTGGTCGGGATGGCGACGATGGGCGGGATGGCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGNTGGCCTNGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCNgg                                                                                                                                               <  1:2064457/150‑1 (MQ=255)
 ccggccgTGGTCGGGATGGCGACGATGGGCGGGATGGCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAgg                                                                                                                                               >  2:579528/1‑150 (MQ=255)
 ccggccgTGGTCGGGATGGCGACGATGGGCGGGATGGCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAgg                                                                                                                                               >  2:2291390/1‑150 (MQ=255)
  cggccgTGGTCGGGATGGCGACGATGGGCGGGATGGCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGc                                                                                                                                              >  2:2203109/1‑150 (MQ=255)
     ccgTGGTCGGGATGGCGACGATGGGCGGGATGGCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGt                                                                                                                                           <  1:764155/150‑1 (MQ=255)
     ccgTGGTCGGGATGGCGACGATGGGCGGGATGGCATGGGTAAGCGAACAACACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGt                                                                                                                                           <  1:145809/150‑1 (MQ=255)
       gTGGTCGGGATGGCGACGATGGGCGGGATGGCATGGGTAAGCGTGAACGACACTGCGGAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCTGCTGGCCGtgg                                                                                                                                         <  2:566167/150‑3 (MQ=255)
               gATGGCGACGATGGGCGGGATGGCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGt                                                                                                                                 >  2:1778369/1‑150 (MQ=255)
               gATGGCGACGATGGGCGGGATGGCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGt                                                                                                                                 <  2:2204135/150‑1 (MQ=255)
                     gacgaTGGGCGGGATGGCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCg                                                                                                                           >  1:922927/1‑150 (MQ=255)
                            ggCGGGATGGCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAgcgc                                                                                                                    >  2:1667923/1‑150 (MQ=255)
                             gCGGGATGGCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCt                                                                                                                   >  2:2174874/1‑150 (MQ=255)
                                 gATGGCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGc                                                                                                               >  2:2279545/1‑150 (MQ=255)
                                    ggCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGCTAGCGCTGGGCGgt                                                                                                            <  2:1851992/150‑1 (MQ=255)
                                    ggCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGgt                                                                                                            >  2:349854/1‑150 (MQ=255)
                                       aTGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTTTGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGtgtg                                                                                                         >  2:2301362/1‑150 (MQ=255)
                                        tGGGTAAGAGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCGGCAGGCTACTGGCCGTGGCCGCATCTAATGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTgg                                                                                                        <  1:168755/150‑1 (MQ=255)
                                         gggTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGt                                                                                                       <  2:1258556/150‑1 (MQ=255)
                                             aaGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGAt                                                                                                   <  2:922927/150‑1 (MQ=255)
                                                cGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATgtc                                                                                                >  1:170080/1‑150 (MQ=255)
                                                  tGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTgg                                                                                                        >  1:2124163/1‑140 (MQ=255)
                                                  tGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTgg                                                                                                        <  2:2124163/140‑1 (MQ=255)
                                                       gACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTg                                                                                         >  1:1809116/1‑150 (MQ=255)
                                                          ccTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTgcgc                                                                                      <  1:160883/150‑1 (MQ=255)
                                                          ccTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTgcgc                                                                                      >  2:2196935/1‑150 (MQ=255)
                                                          ccTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTgcgc                                                                                      >  1:2404356/1‑150 (MQ=255)
                                                          ccTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTATAGACATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTgcgc                                                                                      >  2:1401794/1‑150 (MQ=255)
                                                            ttCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATa                                                                                                                                      <  1:422631/100‑1 (MQ=255)
                                                            ttCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATa                                                                                                                                      >  2:422631/1‑100 (MQ=255)
                                                            ttCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGAt                                                                                                   <  2:934356/135‑1 (MQ=255)
                                                            ttCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGAt                                                                                                   >  1:934356/1‑135 (MQ=255)
                                                                   tCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGt                                                                                                       <  2:2308662/124‑1 (MQ=255)
                                                                   tCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGt                                                                                                       >  1:2308662/1‑124 (MQ=255)
                                                                   tCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGcc                                                                             >  2:905289/1‑150 (MQ=255)
                                                                     cTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCtt                                                                           <  2:1809116/150‑1 (MQ=255)
                                                                     cTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCtt                                                                           >  2:932482/1‑150 (MQ=255)
                                                                      tCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGa                                                                                                    >  1:851756/1‑124 (MQ=255)
                                                                      tCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGa                                                                                                    <  2:851756/124‑1 (MQ=255)
                                                                        gATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGcccc                                                                                                                            <  1:584779/98‑1 (MQ=255)
                                                                        gATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGcccc                                                                                                                            >  2:584779/1‑98 (MQ=255)
                                                                            cGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTcc                                                                    <  1:2291390/150‑1 (MQ=255)
                                                                             ggCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCg                                                                   >  1:16551/1‑150 (MQ=255)
                                                                             ggCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCg                                                                   >  1:826502/1‑150 (MQ=255)
                                                                             ggCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGACGCGACGATCTCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCg                                                                   >  2:1902802/1‑150 (MQ=255)
                                                                              gCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGc                                                                  <  2:2077535/150‑1 (MQ=255)
                                                                              gCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGc                                                                  <  1:905289/150‑1 (MQ=255)
                                                                              gCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGc                                                                  <  1:579528/150‑1 (MQ=255)
                                                                                cgtcgtGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCa                                                                <  1:2083132/150‑1 (MQ=255)
                                                                                 gtcgtGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCAc                                                               >  2:1569947/1‑150 (MQ=255)
                                                                                 gtcgtGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCAc                                                               <  2:452206/150‑1 (MQ=255)
                                                                                    gtGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTc                                                            <  2:962079/150‑1 (MQ=255)
                                                                                    gtGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATAGAGCTGAAGAGCACACGTCTGAGCAGGCGTTCCGGACCCTGGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTc                                                            >  1:1487631/1‑150 (MQ=255)
                                                                                     tGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTcc                                                           <  1:2196935/150‑1 (MQ=255)
                                                                                     tGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTcc                                                           <  1:1989993/150‑1 (MQ=255)
                                                                                      gggTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGt                                                                                                       >  1:1281179/1‑105 (MQ=255)
                                                                                      gggTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGt                                                                                                       <  2:1281179/105‑1 (MQ=255)
                                                                                      gggTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCt                                                          <  1:1667923/150‑1 (MQ=255)
                                                                                       ggTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTc                                                         >  2:1982407/1‑150 (MQ=255)
                                                                                        gTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCTGCTGGCCGtg                                                                                                                                          >  1:1074701/1‑67 (MQ=255)
                                                                                        gTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCTGCTGGCCGtg                                                                                                                                          <  2:1074701/68‑2 (MQ=255)
                                                                                         tGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGt                                                       <  1:349854/150‑1 (MQ=255)
                                                                                           agcagcAGGCTGATGGCCTTTGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAg                                                     <  1:2135929/150‑1 (MQ=255)
                                                                                            gcagcaGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAgg                                                    <  2:1648610/150‑1 (MQ=255)
                                                                                                   gCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCCCCGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATAcc                                             >  2:46350/1‑150 (MQ=255)
                                                                                                    cTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGACAGCAGGCCGTCTATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTAGTTACTACGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCg                                            <  1:436236/150‑1 (MQ=255)
                                                                                                     tGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCAACGGCCAGCAGGCCGTCGATACCCGTGCCCAGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCAACTTCCTCGTAGGCGATACCgg                                           >  2:38436/1‑150 (MQ=255)
                                                                                                          gCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATgtc                                                                                                <  2:284854/92‑1 (MQ=255)
                                                                                                          gCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATgtc                                                                                                >  1:284854/1‑92 (MQ=255)
                                                                                                                gCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTGGATGTCCGCCACTTCCTCGTCGGCGATACCGGCCTTGTCGacg                                >  2:123195/1‑148 (MQ=255)
                                                                                                                   gCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCTAGGCTg                             <  1:1401794/150‑1 (MQ=255)
                                                                                                                   gCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTg                             >  1:31368/1‑150 (MQ=255)
                                                                                                                    catccatcGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTgg                            >  2:2028823/1‑150 (MQ=255)
                                                                                                                      tccatcGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTGGCg                          <  1:932482/150‑1 (MQ=255)
                                                                                                                         atcGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTGGCGTAg                       <  2:2046766/150‑1 (MQ=255)
                                                                                                                         atcGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTGGCGTAg                       >  1:312234/1‑150 (MQ=255)
                                                                                                                            gcgcTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTGGCGTacaca                    >  2:603777/1‑149 (MQ=255)
                                                                                                                            gcgcTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTGGCGTAGACg                    >  2:603497/1‑150 (MQ=255)
                                                                                                                                tGCCGCCACCCACGGCCAGCAGGCTGTCGATACCCGTGCCCCGATAGCGCTGGGAGGCATGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTGGCGTAGACGCTGt                <  1:983259/150‑1 (MQ=255)
                                                                                                                                tGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTGGCGTAGACGCTGt                >  1:195923/1‑150 (MQ=255)
                                                                                                                                tGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTGGCGTAGACGCTGt                >  1:632174/1‑150 (MQ=255)
                                                                                                                                      cacccacGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGGCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTGGCGTAGACGCTGTCcagca          <  1:2174874/150‑1 (MQ=255)
                                                                                                                                       acccacGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTGGCGTAGACGCTGTCcagcag         <  2:826502/150‑1 (MQ=255)
                                                                                                                                          gacGGCCTGCTGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCtt                                                                           >  1:2099277/2‑81 (MQ=255)
                                                                                                                                          gacGGCCTGCTGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCtt                                                                           <  2:2099277/80‑1 (MQ=255)
                                                                                                                                          cacGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTGGCGTAGACGCTGTc               <  1:238558/141‑1 (MQ=255)
                                                                                                                                          cacGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTGGCGTAGACGCTGTc               >  2:238558/1‑141 (MQ=255)
                                                                                                                                              gCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTGGCGTAGACGCTGTCCAGCAGGCctgct  >  1:2227877/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCGGCCGTGGTCGGGATGGCGACGATGGGCGGGATGGCATGGGTAAGCGTGAACGACCCTTCGTAGTCCTCGATACGGCCGTCGTGGGTGAGCAGCAGGCTGATGGCCTTGGCCGCATCCATCGCGCTGCCGCCACCCACGGCCAGCAGGCCGTCGATACCCGTGCCCCGATAGCGCTGGGCGGTGTGGTTGATGTCGTCACTGCGCGGGTTGGCCTTGATGTCCGCCACTTCCTCGTAGGCGATACCGGCCTTGTCGAGGCTGGCGTAGACGCTGTCCAGCAGGCCTGCT  >  NC_002947/3195663‑3195954

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