Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,732,472 G→C 5.7% T28T (ACC→ACG hutI ← imidazolonepropionase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,732,4720GC 5.7% 219.4 / 3.7 70T28T (ACC→ACGhutIimidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base G (38/28);  new base C (2/2);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

CTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTCGGCCGGCTCCGCTACCGGCTGGCCCTGCCCTGCG  >  NC_002947/5732329‑5732619
                                                                                                                                               |                                                                                                                                                    
cTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGAc                                                                                                                                                 <  2:344411/149‑1 (MQ=255)
      ttGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATcgc                                                                                                                                           >  6:12823/1‑149 (MQ=255)
                  acCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGAGCGGGCCTATCCACACGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGAt                                                                                                                               >  2:303570/1‑149 (MQ=255)
                         gggTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAg                                                                                                                        >  4:118478/1‑149 (MQ=255)
                          ggTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAg                                                                                                                        <  3:118478/148‑1 (MQ=255)
                             gtgGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGc                                                                                                                    <  4:141649/149‑1 (MQ=255)
                             gtgGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCAATCGATAAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGc                                                                                                                    >  1:280252/1‑149 (MQ=255)
                              tgGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGcc                                                                                                                           <  6:118105/141‑1 (MQ=255)
                              tgGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGcc                                                                                                                           >  5:118105/1‑141 (MQ=255)
                                   gTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGt                                                                                                                                                    >  6:158782/1‑111 (MQ=255)
                                   gTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGt                                                                                                                                                    <  5:158782/111‑1 (MQ=255)
                                           gCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATg                                                                                                      >  6:6456/1‑149 (MQ=255)
                                               gggAATGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTg                                                                                                  >  1:291604/1‑149 (MQ=255)
                                                    tGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACa                                                                                            <  5:6456/149‑1 (MQ=255)
                                                    tGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACa                                                                                            <  5:44854/149‑1 (MQ=255)
                                                      aCCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATg                                                                                          <  1:113128/149‑1 (MQ=255)
                                                        ccAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATgg                                                                                         >  7:176642/1‑148 (MQ=255)
                                                         cAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGc                                                                                        >  6:189417/1‑148 (MQ=255)
                                                              accaccCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAAtgct                                                                                  <  8:176642/149‑1 (MQ=255)
                                                                  cccAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAAtgctgc                                                                                >  2:94594/1‑147 (MQ=255)
                                                                   ccAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAAtgctgc                                                                                >  1:285013/1‑146 (MQ=255)
                                                                   ccAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAAtgctgc                                                                                >  8:226771/1‑146 (MQ=255)
                                                                    cAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCaga                                                                            <  1:94594/149‑1 (MQ=255)
                                                                    cAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCaga                                                                            <  5:12823/149‑1 (MQ=255)
                                                                          aaCCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTtc                                                                      >  4:29775/1‑149 (MQ=255)
                                                                          aaCCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGTTGCAACCATGGccgatg                                                                                   <  3:282121/136‑6 (MQ=255)
                                                                            ccGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTtctc                                                                    <  1:303570/149‑1 (MQ=255)
                                                                            ccGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTtctc                                                                    <  2:285013/149‑1 (MQ=255)
                                                                              gTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCAt                                                                  <  2:291604/149‑1 (MQ=255)
                                                                               tGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTg                                                                                                                                                      >  1:57911/1‑65 (MQ=255)
                                                                               tGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTg                                                                                                                                                      <  2:57911/65‑1 (MQ=255)
                                                                                  ggTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAAtgctgc                                                                                <  6:188351/131‑1 (MQ=255)
                                                                                  ggTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAAtgctgc                                                                                >  5:188351/1‑131 (MQ=255)
                                                                                      aGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAAt                                                           <  5:189417/148‑1 (MQ=255)
                                                                                         ttCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTcc                                                       >  7:321928/1‑149 (MQ=255)
                                                                                           cGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCTGCCATCGGCCATGGTTGCaa                                                                                               <  8:194525/108‑1 (MQ=255)
                                                                                           cGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCtt                                                     >  8:108801/1‑149 (MQ=255)
                                                                                           cGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCAGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCaa                                                                                               >  7:194525/1‑108 (MQ=255)
                                                                                            gACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCtt                                                     <  7:108801/148‑1 (MQ=255)
                                                                                              cGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGc                                                  >  3:100898/1‑149 (MQ=255)
                                                                                              cGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGc                                                  >  1:330911/1‑149 (MQ=255)
                                                                                               gggCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGc                                                  >  2:91602/1‑148 (MQ=255)
                                                                                                 gcgcCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAg                                               >  3:50091/1‑149 (MQ=255)
                                                                                                    tcAGCTGGGCGCGGGGGCCTATCCACTCGATCAGCCCGGCGCGGGTGTCGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCgg                                            <  8:321928/148‑1 (MQ=255)
                                                                                                       gCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGc                                         <  3:29775/149‑1 (MQ=255)
                                                                                                          cGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCttt                                                    >  6:217550/1‑135 (MQ=255)
                                                                                                          cGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCttt                                                    <  5:217550/135‑1 (MQ=255)
                                                                                                           ggCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCt                                     >  2:291885/1‑149 (MQ=255)
                                                                                                            gcgcgcGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGCTGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGAGGTGCAGCTc                                    >  5:246890/1‑149 (MQ=255)
                                                                                                                      cTATCCACTCGATCAGCCCGGCTCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCTTGCACAATCTCCTTTTCTAGCGGAGCAGCAGGGACGGCTcc                          >  8:310416/1‑149 (MQ=255)
                                                                                                                      cTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTcc                                                       <  6:102713/120‑1 (MQ=255)
                                                                                                                      cTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTcc                                                       >  5:102713/1‑120 (MQ=255)
                                                                                                                       tATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTCGGCCGGCTCCg                         >  3:46867/1‑149 (MQ=255)
                                                                                                                         tCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTcggcc                                <  5:82542/140‑1 (MQ=255)
                                                                                                                         tCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTcggcc                                >  6:82542/1‑140 (MQ=255)
                                                                                                                         tCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTCGGCCGGCTCCg                         <  1:91602/147‑1 (MQ=255)
                                                                                                                         tCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTCGGCCGGCTCCGc                        >  2:165498/1‑148 (MQ=255)
                                                                                                                            aCTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTCGGCCGGCTCCGCTAcc                    >  6:362142/1‑149 (MQ=255)
                                                                                                                            aCTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTCGGCAGGCTCCGCTa                      >  8:194018/1‑147 (MQ=255)
                                                                                                                             cTCGATCAGCCCGGCGCTGGTGACGATCGCAGCATCCTCGATGGCCCAGTAGCCGCCATCGGCCATGGTTGCAACATGGCAATGCTGGCATAGGGTTCTAATGCACAATCTACTTTGCTAGCGGAGCAGCTCGGCACGCTCAGCTACCg                   >  1:204076/1‑149 (MQ=255)
                                                                                                                              tCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGt                                           <  8:93422/124‑1 (MQ=255)
                                                                                                                              tCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGt                                           >  7:93422/1‑124 (MQ=255)
                                                                                                                                 aTCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTCCAGCTCGGCCGGCTCCGCTACCGGCTg               <  7:226771/149‑1 (MQ=255)
                                                                                                                                   ggatgcGGCGATCGTCACCAGCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTCGGCCGGCTCCGCTACCgg                  <  7:176405/139‑1 (MQ=255)
                                                                                                                                   ggatgcGGCGATCGTCACCAGCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTCGGCCGGCTCCGCTACCgg                  >  8:176405/6‑144 (MQ=255)
                                                                                                                                   caacccGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTCGGCCGGCTCCGCTACCGGCTGgc             >  2:107911/4‑149 (MQ=255)
                                                                                                                                      cccGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGtt                                                                                                   >  2:193560/1‑61 (MQ=255)
                                                                                                                                      cccGGCGCTGGTGACGATCGCCGCATCCTCGAAGGCCGAGTAGCGGCCATCGGCCATGGtt                                                                                                   <  1:193560/61‑1 (MQ=255)
                                                                                                                                       gcGGCGATCGTCACCAGCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTCGGCCGGCTCCGCTACCGGCTgg              <  4:180762/143‑1 (MQ=255)
                                                                                                                                       gcGGCGATCGTCACCAGCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTCGGCCGGCTCCGCTACCGGCTgg              >  3:180762/2‑144 (MQ=255)
                                                                                                                                        cGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTCGGCCGGCTCCGCTACCGGCTGgccctgc        >  5:274460/1‑149 (MQ=255)
                                                                                                                                              tGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTCGGCCGGCTCCGCTACCGGCTGGCCCTGCCCTgcg  >  3:121378/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTCGACGGGCGCCAGCTCGGCGCGCGGGCCTATCCACTCGATCAGCCCGGCGCTGGTGACGATCGCCGCATCCTCGATGGCCGAGTAGCGGCCATCGGCCATGGTTGCAACATGGCAATGCTGCCAGAGGGTTCTCATGCACAATCTCCTTTGCTAGCGGTGCAGCTCGGCCGGCTCCGCTACCGGCTGGCCCTGCCCTGCG  >  NC_002947/5732329‑5732619

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