Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 87,493 T→A 7.5% N98I (AAT→ATT)  PP_0076 ← choline betaine‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294787,4930TA 7.5% 247.4 / 10.5 80N98I (AAT→ATT) PP_0076choline betaine‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (40/34);  new base A (3/3);  total (43/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 = 6.34e-01

GATTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTAACAGG  >  NC_002947/87348‑87642
                                                                                                                                                 |                                                                                                                                                     
gATTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGc                                                                                                                                                   <  1:1331748/150‑1 (MQ=255)
   tCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGAt                                                                                                                                                <  1:620792/150‑1 (MQ=255)
   tCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGAt                                                                                                                                                <  2:726011/150‑1 (MQ=255)
     cGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGTGCACTGCGAGGTTGAACTCGGTTCCCTCCAAGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGa                                                                                                                                              >  1:778366/1‑150 (MQ=255)
          aGTTTTTTGTCAAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTg                                                                                                                                         >  2:1042730/1‑150 (MQ=255)
             tttttGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCa                                                                                                                                      <  2:678395/150‑1 (MQ=255)
                    cGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGc                                                                                                                               <  2:1524263/150‑1 (MQ=255)
                     gAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGcc                                                                                                                              >  1:1528981/1‑150 (MQ=255)
                          gCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTtgct                                                                                                                         <  2:1139536/150‑1 (MQ=255)
                              ggCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTgcgc                                                                                                                     >  1:1586535/1‑150 (MQ=255)
                                     ttttGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATc                                                                                                              >  1:1924496/1‑150 (MQ=255)
                                         gCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGt                                                                                                          <  1:1042730/150‑1 (MQ=255)
                                          cAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGaa                                                                                                                                 <  1:1397416/126‑1 (MQ=255)
                                          cAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGaa                                                                                                                                 >  2:1397416/1‑126 (MQ=255)
                                           aGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGttt                                                                                                        >  1:601527/1‑150 (MQ=255)
                                             gTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTc                                                                                                                                       <  2:67963/117‑1 (MQ=255)
                                             gTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTc                                                                                                                                       >  1:67963/1‑117 (MQ=255)
                                                aGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAg                                                                                                   >  2:1015224/1‑150 (MQ=255)
                                                aGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAg                                                                                                   >  2:946394/1‑150 (MQ=255)
                                                aGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAg                                                                                                   >  1:1801880/1‑150 (MQ=255)
                                                   gAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGaa                                                                                                <  2:1801880/150‑1 (MQ=255)
                                                      gTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGc                                                                                             >  2:1347695/1‑150 (MQ=255)
                                                      gTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGc                                                                                             <  2:1503036/150‑1 (MQ=255)
                                                       tCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGa                                                                                                    >  2:705664/1‑142 (MQ=255)
                                                       tCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGa                                                                                                    <  1:705664/142‑1 (MQ=255)
                                                        ccTTCACCCTGGTATTCCACACATAGTCGCGCACTGCGAGGGTCCACCCGGTTCCCTGCAGTTTGCGCGACAGCTGCTGCACATCGCCATTGGCGGCTAACTTGCCATGGAAGCCTTGCTGCGCCCCGATCCAGTTTCCGAGGAATGCGt                                                                                           <  2:778366/150‑1 (MQ=255)
                                                        ccTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGt                                                                                           >  2:119199/1‑150 (MQ=255)
                                                         cTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCAATGGCGATgtg                                                                                                                                             <  1:1495232/99‑3 (MQ=255)
                                                         cTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCAATGGCGATgtg                                                                                                                                             <  1:1806768/99‑3 (MQ=255)
                                                         cTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCAATGGCGATgtg                                                                                                                                             >  2:1806768/1‑97 (MQ=255)
                                                         cTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCAATGGCGATgcg                                                                                                                                             >  2:1495232/1‑97 (MQ=255)
                                                            cACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGAc                                                                                       <  2:1586535/150‑1 (MQ=255)
                                                             aCCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGATCTCGGTTCCCTCCATGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGAcc                                                                                      >  1:1847665/1‑150 (MQ=255)
                                                               cccGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTg                                                                                    >  2:16203/1‑150 (MQ=255)
                                                                     aTTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTc                                                                              >  1:70855/1‑150 (MQ=255)
                                                                          acacaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGa                                                                         <  2:70855/150‑1 (MQ=255)
                                                                            acaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAg                                                                        <  2:1810081/149‑1 (MQ=255)
                                                                            acaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAg                                                                        >  1:1810081/1‑149 (MQ=255)
                                                                              aTAGTCGGGCACTGCGCGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGgcc                                                                     >  1:1701147/1‑150 (MQ=255)
                                                                              aTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGgcc                                                                     >  1:1709584/1‑150 (MQ=255)
                                                                                   cGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGt                                                                <  2:1064860/150‑1 (MQ=255)
                                                                                      gCACTGCGAGTGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAACCCTTGCTGCGCCGGCATCCAGTTTCCGCGGAAGCCGTCGCCCTGGCCGTCCTCGAGGCCACCCTAga                                                             >  1:120704/1‑150 (MQ=255)
                                                                                          tGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAgatgat                                                         >  1:1018144/1‑150 (MQ=255)
                                                                                             gAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTgg                                                                                   <  2:1126187/121‑1 (MQ=255)
                                                                                             gAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTgg                                                                                   >  1:1126187/1‑121 (MQ=255)
                                                                                                 gTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGccgccg                                                                 >  2:857050/1‑135 (MQ=255)
                                                                                                 gTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGccgccg                                                                 <  1:857050/135‑1 (MQ=255)
                                                                                                  tGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTa                                                               <  2:1212890/136‑1 (MQ=255)
                                                                                                  tGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTa                                                               >  1:1212890/1‑136 (MQ=255)
                                                                                                   gAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGc                                                >  1:326386/1‑150 (MQ=255)
                                                                                                     aCTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCa                                              >  2:353886/1‑150 (MQ=255)
                                                                                                      cTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAg                                             >  1:705709/1‑150 (MQ=255)
                                                                                                       tCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAgg                                            >  1:927370/1‑150 (MQ=255)
                                                                                                       tCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAgg                                            >  1:765788/1‑150 (MQ=255)
                                                                                                               cTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGa                                                                         >  2:1776547/1‑113 (MQ=255)
                                                                                                               cTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGa                                                                         <  1:1776547/113‑1 (MQ=255)
                                                                                                                tCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATc                                   <  1:1347695/150‑1 (MQ=255)
                                                                                                                tCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATc                                   <  1:946394/150‑1 (MQ=255)
                                                                                                                   aGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGt                                                                                           >  2:1666433/1‑91 (MQ=255)
                                                                                                                   aGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGt                                                                                           <  1:1666433/91‑1 (MQ=255)
                                                                                                                     gTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGATGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCAc                              <  2:1727077/150‑1 (MQ=255)
                                                                                                                     gTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCAc                              <  2:1701147/150‑1 (MQ=255)
                                                                                                                     gTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTCCCGAg                                                                                                   >  2:1701363/1‑81 (MQ=255)
                                                                                                                     gTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTCCCGAg                                                                                                   <  1:1701363/81‑1 (MQ=255)
                                                                                                                     gTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTCGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGCCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCAc                              <  1:658249/150‑1 (MQ=255)
                                                                                                                          gcgACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAga                                                             >  1:393772/1‑114 (MQ=255)
                                                                                                                          gcgACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAga                                                             <  2:393772/114‑1 (MQ=255)
                                                                                                                          gcgACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGt                         >  2:753995/1‑150 (MQ=255)
                                                                                                                          gcgACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGt                         >  2:342642/1‑150 (MQ=255)
                                                                                                                              cAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAAcc                     <  2:705709/150‑1 (MQ=255)
                                                                                                                                 ctgctgCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAg                  >  1:649937/1‑150 (MQ=255)
                                                                                                                                      gttcaTCGCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCa                                              >  1:1130382/4‑117 (MQ=255)
                                                                                                                                      gttcaTCGCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCa                                              <  2:1130382/114‑1 (MQ=255)
                                                                                                                                       cacaTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTc            <  2:1924496/150‑1 (MQ=255)
                                                                                                                                         caTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCa          <  2:1847665/150‑1 (MQ=255)
                                                                                                                                          aTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTACCCCAGGCCTTCCAg         >  2:1696496/1‑150 (MQ=255)
                                                                                                                                            cGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTa       >  2:1062400/1‑150 (MQ=255)
                                                                                                                                            cGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTa       >  1:1932819/1‑150 (MQ=255)
                                                                                                                                             gCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTaa      <  2:326386/150‑1 (MQ=255)
                                                                                                                                              ccATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTAAc     <  1:753995/150‑1 (MQ=255)
                                                                                                                                                 ttGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTAACAgg  >  1:713038/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GATTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTAACAGG  >  NC_002947/87348‑87642

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