Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,069,546 A→G 7.6% V21A (GTG→GCG) ‡ 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,5460AG 7.6% 186.3 / 5.0 67V21A (GTG→GCG) ‡PP_5323M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base A (27/34);  new base G (2/3);  total (29/38)
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.85e-01

GCCACCGTCGAAGGTTTTCGCGTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAAC  >  NC_002947/6069398‑6069692
                                                                                                                                                    |                                                                                                                                                    
gCCACCGTCGAAGGTTTTCGCGTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTAGGTGATCTTGTAGAGCGTCa                                                                                                                                                      >  4:273157/1‑149 (MQ=255)
      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                                                                                                                               <  3:338670/137‑7 (MQ=255)
                                   ggCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCGAGCCAGGGCGTGACGCtctacaa                                                                                                                               >  4:338670/1‑131 (MQ=255)
                                    gCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCGTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGcac                                                                                                                <  6:322120/149‑1 (MQ=255)
                                      cTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAg                                                                                                              <  7:155735/149‑1 (MQ=255)
                                            gCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGccc                                                                                                                                                 >  2:283885/1‑110 (MQ=255)
                                            gCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGccc                                                                                                                                                 <  1:283885/110‑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                                                                                 <  3:54876/144‑1 (MQ=255)
                                                                        cAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGt                                                                                 >  4:54876/1‑144 (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                                                     >  2:172471/1‑141 (MQ=255)
                                                                                                       gAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGaa                                                     <  1:172471/141‑1 (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)
                                                                                                                                                    |                                                                                                                                                    
GCCACCGTCGAAGGTTTTCGCGTTCAGCAGTACCTGGCCTTCGAGCGTTTCGACCATGGGCTCGAGCAACACCAAGGGCTTGGCGCCGGGGCCAGGGCGATCGGAATAGGAGATGACACCGGAATCGTCGGTGATCTTGTAGAGCGTCACGCCCTGGCTCGACGCTGAGGCCATGAGCAATGCACAGAACAGCAGCAGGCGCGCGGGCATGATAGTGGCTTTTTGACAGTTATGATCTGTCGAAGCCTAGCAGCGGGTTTTACGGTAGGTATTGCAGTTGGTCGGATGGGGCAAC  >  NC_002947/6069398‑6069692

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: