Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,069,552 T→G 7.6% Q19P (CAG→CCG) ‡ PP_5323 ← M23/M37 family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,069,5520TG 7.6% 193.8 / 7.4 66Q19P (CAG→CCG) ‡PP_5323M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base T (27/34);  new base G (2/3);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

GTCGAAGGTTTTCGCGTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAACTCGGGCTT  >  NC_002947/6069404‑6069700
                                                                                                                                                    |                                                                                                                                                      
gTCGAAGGTTTTCGCGTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCt                                                                                                                                                        >  8:155735/1‑149 (MQ=255)
           tCGCGGTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCTAGGGCTTGGCGCCGGGGCCAGGGCGATCTGAATAGGAGATGACACCGGAATCGTCGGTTATCTTGTAGAGCGTCACGCCCTGGCTCGACGc                                                                                                                                              <  7:123839/148‑1 (MQ=255)
             gcgTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCGAGCCAGGGCGTGACGCtct                                                                                                                                           <  4:175855/149‑3 (MQ=255)
             gcgTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGa                                                                                                                                           <  1:109087/149‑1 (MQ=255)
              cgTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAg                                                                                                                                         <  1:156591/149‑1 (MQ=255)
                      agTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAg                                                                                                                                 <  4:297567/149‑1 (MQ=255)
                       gTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCCTCACGCCCTGGCTCGACGCTGAGGCCATGAg                                                                                                                                 >  7:300046/1‑148 (MQ=255)
                          ccTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAAt                                                                                                                             >  4:36221/1‑149 (MQ=255)
                           cTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATg                                                                                                                            <  7:244491/149‑1 (MQ=255)
                           cTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATg                                                                                                                            <  3:273157/149‑1 (MQ=255)
                           cTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGAACTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATg                                                                                                                            >  6:288345/1‑149 (MQ=255)
                             ggCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGTGCCAGGGAGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTAGTAGAGCGTCACGACCTGGCTCGACGCTGAGTCAATGAGAAAtta                                                                                                                           >  8:146393/1‑146 (MQ=255)
                             ggCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCGAGCCAGGGCGTGACGCtctacaa                                                                                                                                       >  4:338670/1‑131 (MQ=255)
                             ggCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCGAGCCAGGGCGTGACGCtctacaa                                                                                                                                       <  3:338670/137‑7 (MQ=255)
                              gCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCGTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGcac                                                                                                                        <  6:322120/149‑1 (MQ=255)
                                cTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAg                                                                                                                      <  7:155735/149‑1 (MQ=255)
                                       cGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAAcagca                                                                                                               <  4:214513/149‑1 (MQ=255)
                                       cGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAAcagca                                                                                                               <  3:36221/149‑1 (MQ=255)
                                             gACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGgc                                                                                                         >  2:219450/1‑149 (MQ=255)
                                             gACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGgc                                                                                                         >  1:294686/1‑149 (MQ=255)
                                              aCCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGgcg                                                                                                        >  4:349323/1‑149 (MQ=255)
                                               ccATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGCGCAATGCACAGAACCGCAGCAGgcgc                                                                                                       >  8:139352/1‑149 (MQ=255)
                                                  tGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGcgcgcg                                                                                                    <  3:327675/149‑1 (MQ=255)
                                                     gCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCggg                                                                                                  >  7:340210/1‑148 (MQ=255)
                                                      cTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAAcagcagcag                                                                                                           >  4:153701/1‑138 (MQ=255)
                                                      cTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAAcagcagcag                                                                                                           <  3:153701/138‑1 (MQ=255)
                                                       tCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGc                                                                                                                                                     <  8:207037/97‑1 (MQ=255)
                                                          aGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCGAGCCAGGGCGTGACGCtctac                                                                                                                                         <  8:227146/106‑5 (MQ=38)
                                                          aGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATATTGTAGAGCGTCGAGCCAGGGCGTGACGCtctac                                                                                                                                         >  7:227146/1‑102 (MQ=37)
                                                            cAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCCACAGAACAGCAGCAGGCGCGCGGGCATGATa                                                                                           <  1:37612/149‑1 (MQ=255)
                                                              acacCAAGGGCTTGGCGCCGGGGCCAGGGCGATCTGAATATGAGATGACACCGCAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCTATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTg                                                                                        <  2:318236/149‑1 (MQ=255)
                                                              acacCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTg                                                                                        <  1:219450/149‑1 (MQ=255)
                                                                  cAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGt                                                                                         >  4:54876/1‑144 (MQ=255)
                                                                  cAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGt                                                                                         <  3:54876/144‑1 (MQ=255)
                                                                  cAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCtt                                                                                    >  5:186300/1‑149 (MQ=255)
                                                                   aaGGGCTTGGCGCTGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCttt                                                                                   <  7:134250/149‑1 (MQ=255)
                                                                    aGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGAc                                                                                                                                                >  3:259858/1‑89 (MQ=255)
                                                                    aGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGAc                                                                                                                                                <  4:259858/89‑1 (MQ=255)
                                                                     gggCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCttttt                                                                                 <  3:349323/149‑1 (MQ=255)
                                                                         ttGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGa                                                                               >  8:41227/1‑147 (MQ=255)
                                                                            gcgcCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTGTTTGACAGt                                                                           >  5:258502/1‑148 (MQ=255)
                                                                                     ccAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGt                                                                 <  8:93153/149‑1 (MQ=255)
                                                                                       aGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCg                                                               <  8:340210/149‑1 (MQ=255)
                                                                                         ggCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGGCg                                                               >  8:334722/1‑147 (MQ=255)
                                                                                          gCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATcaaatcgg                                                             >  3:237554/1‑141 (MQ=255)
                                                                                                 gAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGaa                                                             <  1:172471/141‑1 (MQ=255)
                                                                                                 gAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGaa                                                             >  2:172471/1‑141 (MQ=255)
                                                                                                  aaTAGTAGATGACACCGGAATCTTCGGTGATCTTGTAGAGCGTCTCGCCCTGGCTCGACGCTGATGCCATGAGCAATGCACAGAACAGCACCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTa                                                        <  7:146393/145‑1 (MQ=255)
                                                                                                   aTAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTagcagc                                                   >  1:10094/1‑149 (MQ=255)
                                                                                                        agaTGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGtt                                              <  5:288345/149‑1 (MQ=255)
                                                                                                         gaTGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGttt                                             <  6:186300/149‑1 (MQ=255)
                                                                                                            gACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAg                                                                                                                                         >  6:64565/1‑55 (MQ=255)
                                                                                                            gACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAg                                                                                                                                         <  5:64565/55‑1 (MQ=255)
                                                                                                            gACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTa                                           <  6:305492/148‑1 (MQ=255)
                                                                                                             acacCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACg                                         <  8:29004/149‑1 (MQ=255)
                                                                                                                ccGGAATCGTCGGTGATCTTGTAGAGCGTCACTCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACagta                                      <  6:258502/149‑1 (MQ=255)
                                                                                                                 cGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGAAGCGGGTTTTACggtag                                     >  4:197271/1‑149 (MQ=255)
                                                                                                                 aggAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACggtag                                     >  2:323500/2‑149 (MQ=255)
                                                                                                                               aTCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTgg                        <  2:10094/148‑1 (MQ=255)
                                                                                                                                   tGTAGAGCGTCACGCCCTGGCTCGACGCTGATGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTAtt                                >  1:38032/1‑137 (MQ=255)
                                                                                                                                   tGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTAtt                                <  2:38032/136‑1 (MQ=255)
                                                                                                                                    gTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGAt                  >  6:200259/1‑149 (MQ=255)
                                                                                                                                      agagCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATg                 <  7:242627/148‑1 (MQ=255)
                                                                                                                                         gCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATgggg              <  5:200259/148‑1 (MQ=255)
                                                                                                                                            tCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAAc          >  6:331396/1‑149 (MQ=255)
                                                                                                                                                gCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAACTCgg      >  2:12583/1‑149 (MQ=255)
                                                                                                                                                  ccTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAACTCGGGc    >  4:327150/1‑149 (MQ=255)
                                                                                                                                                    tGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCACCTCGGGCtt  >  5:36058/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GTCGAAGGTTTTCGCGTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAACTCGGGCTT  >  NC_002947/6069404‑6069700

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