Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,538,610 C→G 5.6% A668P (GCG→CCG)  PP_1350 ← two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,538,6100CG 5.6% 232.6 / 3.9 72A668P (GCG→CCG) PP_1350two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base C (34/34);  new base G (2/2);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

ACCTCGCCAATCGCTTCGATCCGCGGGTGGTCCTTGGGCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAGGCCCAGCCCGGTACCCTGGCCGATTTCCTTGGT  >  NC_002947/1538464‑1538757
                                                                                                                                                  |                                                                                                                                                   
aCCTCGCCAATCGCTTCGATCCGCGGGTGGTCCTTGGGCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGc                                                                                                                                                  >  2:737535/1‑150 (MQ=255)
   tCGCCAATCGCTTCGATCCGCGGGTGGTCCTTGGGCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCgg                                                                                                                                               >  2:425500/1‑150 (MQ=255)
      ccAATCGCTTCGATCCGCGGGTGGTCCTTGGGCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCAGCGGCGGTGCCgg                                                                                                                                               >  2:954189/1‑147 (MQ=255)
      ccAATCGCTTCGATCCGCGGGTGGTCCTTGGGCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCAGCGGCGGTGCCgg                                                                                                                                               <  1:954189/147‑1 (MQ=255)
            gCTTCGATCCGCGGGTGGTCCTTGGGCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAg                                                                                                                                      >  1:1410550/1‑150 (MQ=255)
             cTTCGATCCGCGGGTGGTCCTTGGGCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAgg                                                                                                                                     <  1:29360/150‑1 (MQ=255)
             cTTCGATCCGCGGGTGGTCCTTGGGCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAgg                                                                                                                                     <  2:90183/150‑1 (MQ=255)
                         ggtggtCCTTGGGCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAgg                                                                                                                         <  1:994225/150‑1 (MQ=255)
                         ggtggtCCTTGGGCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAgg                                                                                                                         <  1:985537/150‑1 (MQ=255)
                             gtCCTTGGGCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTaa                                                                                                                     >  2:986680/1‑150 (MQ=255)
                               ccTTGGGCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCACATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAg                                                                                                                   >  2:1250021/1‑150 (MQ=255)
                                 ttGGGCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCa                                                                                                                 <  1:737535/150‑1 (MQ=255)
                                    ggCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGt                                                                                                              <  2:390099/150‑1 (MQ=255)
                                    ggCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGt                                                                                                              <  1:700466/150‑1 (MQ=255)
                                         tCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCtt                                                                                                         >  2:1169333/1‑150 (MQ=255)
                                          cGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTg                                                                                                        <  2:1410550/150‑1 (MQ=255)
                                              ggCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCt                                                                                                    >  1:713588/1‑150 (MQ=255)
                                                ccGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTgg                                                                                                  >  1:1376066/1‑150 (MQ=255)
                                                      cccAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGa                                                                                            <  1:986680/150‑1 (MQ=255)
                                                      cccAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGACCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGa                                                                                            <  1:92683/150‑1 (MQ=255)
                                                        cAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGAtt                                                                                          >  2:1147804/1‑150 (MQ=255)
                                                              ggTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACtt                                                                                    >  2:772547/1‑150 (MQ=255)
                                                               gTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTc                                                                                   <  1:539185/150‑1 (MQ=255)
                                                                   cGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATc                                                                               <  1:1126375/150‑1 (MQ=255)
                                                                    gCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCt                                                                              <  1:706422/150‑1 (MQ=255)
                                                                    gCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCt                                                                              <  2:280497/150‑1 (MQ=255)
                                                                     ccTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTg                                                                             >  1:230189/1‑150 (MQ=255)
                                                                       tttGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAgg                                                                                                               >  2:1184382/1‑114 (MQ=255)
                                                                       tttGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAgg                                                                                                               <  1:1184382/114‑1 (MQ=255)
                                                                       tttGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGcc                                                                           >  2:1202937/1‑150 (MQ=255)
                                                                        ttGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGccc                                                                          <  2:713588/150‑1 (MQ=255)
                                                                             gCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGACTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTtgt                                                                     <  1:963241/150‑1 (MQ=255)
                                                                                  tGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGCTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGgttg                                                                <  2:919945/150‑1 (MQ=255)
                                                                                    cGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGgttgtt                                                              <  1:669063/150‑1 (MQ=255)
                                                                                    cGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGgttgtt                                                              <  2:1376066/150‑1 (MQ=255)
                                                                                      ggTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTgg                                                            <  2:1260802/150‑1 (MQ=255)
                                                                                       gTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGt                                                           >  1:1296757/1‑150 (MQ=255)
                                                                                           aGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGtgatg                                                       >  2:1004697/1‑150 (MQ=255)
                                                                                           aGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGtgatg                                                       >  2:426528/1‑150 (MQ=255)
                                                                                            gATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTggtc                                                                                                     <  1:1214539/103‑3 (MQ=255)
                                                                                            gATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTggtc                                                                                                     >  2:1214539/1‑101 (MQ=255)
                                                                                            gATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCGATCTGCCCCTTGTGGTTGTTGGtgatga                                                      >  2:586456/1‑150 (MQ=255)
                                                                                             aTCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAgg                                                                                                                         <  1:681584/82‑1 (MQ=255)
                                                                                             aTCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAgg                                                                                                                         >  2:681584/1‑82 (MQ=255)
                                                                                             aTCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGgt                                                                  <  2:1299022/137‑1 (MQ=255)
                                                                                             aTCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGgt                                                                  >  1:1299022/1‑137 (MQ=255)
                                                                                             aTCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGtgatgat                                                     >  2:1116340/1‑150 (MQ=255)
                                                                                             aTCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGACGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCa                                                                                       >  1:264150/1‑116 (MQ=255)
                                                                                             aTCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGACGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCa                                                                                       <  2:264150/116‑1 (MQ=255)
                                                                                                   gAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGt                                                                                                              >  1:499837/1‑87 (MQ=255)
                                                                                                   gAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGt                                                                                                              <  2:499837/87‑1 (MQ=255)
                                                                                                               cccATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAg                                   <  1:586456/150‑1 (MQ=255)
                                                                                                                 cATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAGGc                                 >  1:542033/1‑150 (MQ=255)
                                                                                                                       aTTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAGGCCCAGcc                           >  2:781482/1‑150 (MQ=255)
                                                                                                                       aTTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAGGCCCAGcc                           >  2:1151459/1‑150 (MQ=255)
                                                                                                                       aTTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAGGCCCAGcc                           <  1:772547/150‑1 (MQ=255)
                                                                                                                            tCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGtgatga                                                      <  1:391290/118‑1 (MQ=255)
                                                                                                                            tCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGtgatga                                                      >  2:391290/1‑118 (MQ=255)
                                                                                                                            tCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGga                                            >  2:520459/1‑128 (MQ=255)
                                                                                                                            tCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGga                                            <  1:520459/128‑1 (MQ=255)
                                                                                                                              cGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCa                                                                                       >  1:699365/1‑83 (MQ=255)
                                                                                                                              cGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCa                                                                                       <  2:699365/83‑1 (MQ=255)
                                                                                                                               gTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAGGCCCAGCCCGGTAccc                   >  1:1325723/1‑150 (MQ=255)
                                                                                                                                ttCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAGGCCCAGCCCGGTACCCt                  <  1:1147804/150‑1 (MQ=255)
                                                                                                                                   tGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAGGCCCAGCCCGGTACCCTGGc               <  2:542033/150‑1 (MQ=255)
                                                                                                                                   cagCCGGCAGCGGCGGTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGc                                                                                                 <  1:938599/66‑1 (MQ=37)
                                                                                                                                   cagCCGGCAGCGGCGGTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGc                                                                                                 >  2:938599/3‑68 (MQ=37)
                                                                                                                                       cGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAGGCCCAGCCCGGTa                      >  1:1146919/1‑139 (MQ=255)
                                                                                                                                       cGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAGGCCCAGCCCGGTa                      <  2:1146919/139‑1 (MQ=255)
                                                                                                                                       cGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTCGGTGATGATGAAGTAGGAGACCGACAGGCCCAGCCCGGTACCCTGGCCGAt           >  2:705326/1‑150 (MQ=255)
                                                                                                                                           aCCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCGCCTTGTGGTTGGTGGTGATGATGAAGTAGGAGACCGACAGGCGCAGCCCGGTACCCTGGCCGATTTcc       >  2:741112/1‑150 (MQ=255)
                                                                                                                                           aCCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAGGCCCAGCCCGGTACCCTGGCCGATTTcc       >  2:122551/1‑150 (MQ=255)
                                                                                                                                               ccgcTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAGGCCCAGCCCGGTACCCTGGCCGATTTCCTtgg   >  2:1183720/1‑150 (MQ=255)
                                                                                                                                                cgcTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAGGCCCAGCCCGGTACCCTGGCCGATTTCCTtggt  <  1:425500/150‑1 (MQ=255)
                                                                                                                                                cgcTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAGGCCCAGCCCGGTACCCTGGCCGATTTCCTtggt  <  1:781482/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACCTCGCCAATCGCTTCGATCCGCGGGTGGTCCTTGGGCACTCGACGGCCGTCGCCCAGGCCGGTTTCGCCTTTGTCGCCGGTGCGGGTATAGATCTTCGAAAGGCGGTAGCCCATGTCATTGCTCCGTTCTGGCCGGCACCGCCGCTGCCGGCTGGGTGAGGGGCAGGCGTAGGGTAAAGCAGGTGCCTTGGCCTGGCGTGGATTGCACTTCCATCTGCCCCTTGTGGTTGTTGGTGATGATGAAGTAGGAGACCGACAGGCCCAGCCCGGTACCCTGGCCGATTTCCTTGGT  >  NC_002947/1538464‑1538757

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