Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,706,336 T→G 8.7% M134R (ATG→AGG)  PP_1500 → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,706,3360TG 8.7% 208.1 / 10.0 69M134R (ATG→AGG) PP_1500lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base G (2/4);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

AACGGCGGTAACCTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATGGTGTGC  >  NC_002947/1706188‑1706482
                                                                                                                                                    |                                                                                                                                                  
aaCGGCGGTAACCTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCAt                                                                                                                                                    >  5:278590/1‑149 (MQ=255)
aaCGGCGGTAACCTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGAGCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCAt                                                                                                                                                    <  6:415970/149‑1 (MQ=255)
          aCCTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCAGGCCATGGCTg                                                                                                                                          <  6:189785/149‑1 (MQ=255)
            cTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTgca                                                                                                                                        >  4:281606/1‑149 (MQ=255)
            cTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTg                                                                                                                                          >  5:301823/1‑147 (MQ=255)
            cTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCAGGCCATGGCTgca                                                                                                                                        <  1:239236/149‑1 (MQ=255)
             tGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTgcag                                                                                                                                       >  6:7255/1‑149 (MQ=255)
             tGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTgcag                                                                                                                                       >  6:33732/1‑149 (MQ=255)
             tGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACAGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGAAATTGCTGCAGCCATGGCCTGGCTgctg                                                                                                                                       >  3:246955/1‑147 (MQ=255)
              gCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTgcag                                                                                                                                       >  1:396654/1‑148 (MQ=255)
               ccAGGGCACCTGGGGGGTAGCTACTGGGTGCCGGAGCAGGGCTACCAGGTGCTGATGATGGATTACCGCGGCTATTGCTTGTCTCAGGGCCAGCCGAGCCTGACGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTgcagca                                                                                                                                     <  1:231701/149‑1 (MQ=255)
               ccAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCAGGCCATGTCTgcag                                                                                                                                       <  8:256883/147‑1 (MQ=255)
                           gggggTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCCGCAGGCGACCGAgg                                                                                                                          >  4:48077/1‑148 (MQ=255)
                           gggggTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAgg                                                                                                                          <  3:129381/148‑1 (MQ=255)
                           gggggTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAg                                                                                                                           <  4:335654/147‑1 (MQ=255)
                           gggggTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGt                                                                                                                         <  6:278590/149‑1 (MQ=255)
                               gTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAg                                                                                                                     <  5:234452/149‑1 (MQ=255)
                               gTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAg                                                                                                                     <  3:281606/149‑1 (MQ=255)
                                  gCTACTGGCTTCCGGAGCAGTGCTACCATGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGAAATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGc                                                                                                                  <  4:246955/149‑1 (MQ=255)
                                     aCTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAg                                                                                                               >  3:144515/1‑149 (MQ=255)
                                             ccGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTgg                                                                                                        >  7:147594/1‑148 (MQ=255)
                                              cGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTa                                                                                                      <  5:33732/149‑1 (MQ=255)
                                               ggAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGGTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGc                                                                                                              <  8:22259/140‑1 (MQ=255)
                                               ggAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGc                                                                                                              >  7:22259/1‑140 (MQ=255)
                                                gAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTAct                                                                                                    >  6:9535/1‑149 (MQ=255)
                                                 aGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTActg                                                                                                   <  5:7255/149‑1 (MQ=255)
                                                       gCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTgggg                                                                                             <  5:135487/149‑1 (MQ=255)
                                                              ggTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGcc                                                                                      >  8:468/1‑149 (MQ=255)
                                                              ggTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAACCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGcc                                                                                      >  1:286943/1‑149 (MQ=255)
                                                               gTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCt                                                                                     >  8:262744/1‑149 (MQ=255)
                                                                 gCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTgg                                                                                   >  3:136897/1‑149 (MQ=255)
                                                                  cTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTggg                                                                                  <  5:9535/149‑1 (MQ=255)
                                                                  cTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTgg                                                                                   <  4:144515/148‑1 (MQ=255)
                                                                     atgatCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCgg                                                                               >  5:177649/1‑149 (MQ=255)
                                                                       gatCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTg                                                                             >  3:114197/1‑149 (MQ=255)
                                                                            aTTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATg                                                                        >  5:137937/1‑149 (MQ=255)
                                                                               aCCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATg                                                                        >  5:129118/1‑146 (MQ=255)
                                                                               aCCGCGGATATGGCTGGTCTCAGGTCCAGGCGAGCCTGCCGGAGGTGTACCTGGACATTGCTGCAGCCATGGGCTGGCTGCAGCAGGCGGCCGAGGTCAAGGGCAAGGCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATg                                                                        <  6:129118/146‑1 (MQ=255)
                                                                                ccGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCca                                                                    >  2:376518/1‑149 (MQ=255)
                                                                                  gcgGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCcatt                                                                  <  4:39891/149‑1 (MQ=255)
                                                                                     gCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCcattcat                                                               <  7:468/149‑1 (MQ=255)
                                                                                               gTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTggcg                                                      <  2:286943/148‑1 (MQ=255)
                                                                                                tctcAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTggcggc                                                    >  2:390656/1‑149 (MQ=255)
                                                                                                tctcAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTggcggc                                                    >  5:371911/1‑149 (MQ=255)
                                                                                                  tcAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGccc                                                  >  2:261578/1‑149 (MQ=255)
                                                                                                  tcAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGcc                                                   >  4:71440/1‑148 (MQ=255)
                                                                                                       gCCAGCCGAGCCTGCCGGAGGTGTACCAGGACAGTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCAcccc                                             <  3:48077/149‑1 (MQ=255)
                                                                                                         cAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGa                                           >  5:384994/1‑149 (MQ=255)
                                                                                                         cAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGa                                           <  4:357538/149‑1 (MQ=255)
                                                                                                          aGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCTag                                          >  4:28505/1‑149 (MQ=255)
                                                                                                          aGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGag                                          >  1:57836/1‑149 (MQ=255)
                                                                                                           gCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGagc                                         <  2:57836/149‑1 (MQ=255)
                                                                                                            ccGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGagca                                        <  1:390656/149‑1 (MQ=255)
                                                                                                            ccGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTACTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGag                                          >  6:354831/1‑147 (MQ=255)
                                                                                                                 ccTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTc                                   >  7:358172/1‑149 (MQ=255)
                                                                                                                     ccGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCcattcatt                                                              >  4:62708/1‑118 (MQ=255)
                                                                                                                     ccGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCcattcatt                                                              <  3:62708/118‑1 (MQ=255)
                                                                                                                     ccGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAgcg                               <  3:71440/149‑1 (MQ=255)
                                                                                                                     ccGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAg                                 <  1:376518/147‑1 (MQ=255)
                                                                                                                            tgtACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAgc                        <  6:371911/149‑1 (MQ=255)
                                                                                                                              tACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAgcgc                      <  8:229283/149‑1 (MQ=255)
                                                                                                                                ccAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGcc                     <  8:147594/148‑1 (MQ=255)
                                                                                                                                ccAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGGCAAGGGCAAGCCTCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGcc                     <  4:206485/148‑1 (MQ=255)
                                                                                                                                  atgACATTGATGCAGCCATGGCCTGGCTTCAGAAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGTGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGTCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTg                  <  2:396654/147‑1 (MQ=255)
                                                                                                                                       cctgctgcAGCCAGGCCATGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGgtgt              >  8:350836/3‑148 (MQ=255)
                                                                                                                                       cctgctgcAGCCAGGCCATGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGtt             <  7:350836/147‑1 (MQ=255)
                                                                                                                                           ctgcAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGAt         <  6:384994/149‑1 (MQ=255)
                                                                                                                                           ctgcAGCCAGGCCATGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGAt         >  6:210733/1‑149 (MQ=255)
                                                                                                                                            tgcAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATg        >  3:399882/1‑149 (MQ=255)
                                                                                                                                                 ccATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATGGtgtg   >  3:23526/1‑149 (MQ=255)
                                                                                                                                                  cATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATGGTGTGc  <  6:137937/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AACGGCGGTAACCTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCCTGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATGGTGTGC  >  NC_002947/1706188‑1706482

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: