Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,041,948 C→A 9.4% T143N (ACC→AAC)  colR → two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,041,9480CA 9.4% 288.8 / 13.0 95T143N (ACC→AAC) colRtwo‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base C (46/40);  new base A (4/5);  total (50/45)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TCGGGGGCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACAGCCTGCGCA  >  NC_002947/1041800‑1042094
                                                                                                                                                    |                                                                                                                                                  
tCGGGGGCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTAcc                                                                                                                                                   >  2:340010/1‑150 (MQ=255)
tCGGGGGCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTAcc                                                                                                                                                   >  1:2014996/1‑150 (MQ=255)
     ggCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCa                                                                                                                                              <  2:2080520/150‑1 (MQ=255)
     ggCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCa                                                                                                                                              <  2:2033957/150‑1 (MQ=255)
      gCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAg                                                                                                                                             <  2:1907613/150‑1 (MQ=255)
      gCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAg                                                                                                                                             <  2:471185/150‑1 (MQ=255)
       cTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGCGGGTAACCTCCAgg                                                                                                                                            <  1:601189/150‑1 (MQ=255)
        tggcgacTACCTGCGCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTACGCCGTGCCCAGGGCGGTGTGCGGCGCACGCAGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAggg                                                                                                                                           >  2:605348/4‑150 (MQ=255)
        tGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAggg                                                                                                                                           >  1:2346316/1‑150 (MQ=255)
        tGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAggg                                                                                                                                           >  1:1680095/1‑150 (MQ=255)
         gacgacTACCTGCTCAAGCCATTTGCCATGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGt                                                                                                                                          <  2:1270693/150‑1 (MQ=255)
            gacTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGTCTTACATGAGCTACTACGTCGATACCCTGGAGGTTACCCGCCAGGGgcgc                                                                                                                                       >  2:2410610/1‑150 (MQ=255)
               tACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGCGGGTAACCTCCAGGGtatcga                                                                                                                                     <  2:1447868/149‑6 (MQ=255)
               tACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGCGGGTAACCTCCAGGGtatcga                                                                                                                                     >  1:1447868/1‑144 (MQ=255)
                aCCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTACTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGTCGGTTACCCGCCAGTGGCGccggg                                                                                                                                   >  1:499974/1‑147 (MQ=255)
                 cctcctCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGtt                                                                                                                                  <  1:231533/150‑1 (MQ=255)
                  cTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGcc                                                                                                                                      <  2:461361/145‑1 (MQ=255)
                  cTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGcc                                                                                                                                      >  1:461361/1‑145 (MQ=255)
                  cTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGt                                                                                                                                   >  1:1399196/1‑148 (MQ=255)
                  cTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGt                                                                                                                                   <  2:1399196/148‑1 (MQ=255)
                    gCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGGGCGACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGaa                                                                                                                               <  1:1023719/150‑1 (MQ=255)
                      tCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTg                                                                                                                                    <  1:1373407/143‑1 (MQ=255)
                      tCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGt                                                                                                                                   <  2:1197292/144‑1 (MQ=255)
                      tCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGt                                                                                                                                   >  1:1197292/1‑144 (MQ=255)
                      tCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTAGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTg                                                                                                                                    >  2:1373407/1‑143 (MQ=255)
                        aaGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGCGGGTAACCTCCAGGGtatcgagg                                                                                                                                   >  2:1805507/1‑135 (MQ=255)
                        aaGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGCGGGTAACCTCCAGGGtatcgagg                                                                                                                                   <  1:1805507/142‑8 (MQ=255)
                        aaGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTc                                                                                                                           >  2:451625/1‑150 (MQ=255)
                        aaGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTc                                                                                                                           >  2:1751985/1‑150 (MQ=255)
                          gCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATCCCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCaa                                                                                                                         <  2:1552117/150‑1 (MQ=255)
                          gCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCaa                                                                                                                         <  2:203583/150‑1 (MQ=255)
                                 gCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTc                                                                                                                  <  1:1118419/150‑1 (MQ=255)
                                 gCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTAGCCTGTTGAAGCTCAACCCgggg                                                                                                                  >  2:2396999/1‑148 (MQ=255)
                                    cTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGc                                                                                                               >  1:489820/1‑150 (MQ=255)
                                         ggAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCaa                                                                                                          >  2:2005536/1‑150 (MQ=255)
                                              tGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAgctgc                                                                                                     <  1:2383420/150‑1 (MQ=255)
                                                 cggcgCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTgg                                                                                                  <  2:2346316/150‑1 (MQ=255)
                                                  ggcgCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGc                                                                                                 <  2:1051157/150‑1 (MQ=255)
                                                       gcATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAGCCCGGTCGGCCTCAAGCTGCTGGCGGTAc                                                                                            <  1:2355990/150‑1 (MQ=255)
                                                       gcATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTAc                                                                                            <  2:2291545/150‑1 (MQ=255)
                                                       gcATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTAc                                                                                            >  2:2196224/1‑150 (MQ=255)
                                                         aTCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTg                                                                                          >  2:901630/1‑150 (MQ=255)
                                                          tCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGa                                                                                         >  1:2062900/1‑150 (MQ=255)
                                                          tCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGa                                                                                         >  1:2248690/1‑150 (MQ=255)
                                                          tCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGa                                                                                         >  1:2405301/1‑150 (MQ=255)
                                                            gAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATg                                                                                       >  2:2418361/1‑150 (MQ=255)
                                                             aaGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGCGGGTAACCTCCAggg                                                                                                                                           >  2:2420811/1‑97 (MQ=255)
                                                             aaGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGCGGGTAACCTCCAggg                                                                                                                                           <  1:2420811/97‑1 (MQ=255)
                                                              aGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCa                                                                                     <  1:758709/150‑1 (MQ=255)
                                                                     ctgcGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGtt                                                                                                                                  <  1:2210645/98‑1 (MQ=255)
                                                                     ctgcGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGtt                                                                                                                                  >  2:2210645/1‑98 (MQ=255)
                                                                     ctgcGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGc                                                                              <  2:2049271/150‑1 (MQ=255)
                                                                      tgcGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGc                                                                                                                                       <  1:1039958/92‑1 (MQ=255)
                                                                      tgcGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGc                                                                                                                                       >  2:1039958/1‑92 (MQ=255)
                                                                       gcgcCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTc                                                                                                                  >  2:2014452/1‑112 (MQ=255)
                                                                       gcgcCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTc                                                                                                                  <  1:2014452/112‑1 (MQ=255)
                                                                                 cAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCGCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTg                                                                  <  1:2396999/150‑1 (MQ=255)
                                                                                          gggCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAg                                                         <  1:2418361/150‑1 (MQ=255)
                                                                                                 gcACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACGCGCCAGGGGCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGCGGTCCTgg                                                  >  2:230671/1‑150 (MQ=255)
                                                                                                   aCCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGgaa                                                <  2:489820/150‑1 (MQ=255)
                                                                                                      cTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAg                                                                                                                              <  2:2091062/69‑1 (MQ=255)
                                                                                                      cTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAg                                                                                                                              >  1:2091062/1‑69 (MQ=255)
                                                                                                       tGCAGGTGGCTGACCTGAGCTACGACCTAGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGCACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGaagaag                                            >  2:1248227/1‑150 (MQ=255)
                                                                                                        gCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGACTGTTGAAGGACAACAAGGTCGGCCTCAAGATCCTGGCGGTACTGAGGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGc                                           <  2:1465408/150‑1 (MQ=255)
                                                                                                         cAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGc                                                     <  2:1622785/139‑1 (MQ=255)
                                                                                                         cAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGc                                                     <  2:1622778/139‑1 (MQ=255)
                                                                                                         cAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGc                                                     >  1:1622785/1‑139 (MQ=255)
                                                                                                         cAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGc                                                     >  1:1622778/1‑139 (MQ=255)
                                                                                                            gTGGCTGACCTGAGCTACGACCTCGCTGGCCTCGAGGTTACCCGCCAGGGTCGCCTGTTGCCGCTCAACCCGCTCGACCTCAACCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTtg                                       >  1:663389/1‑150 (MQ=255)
                                                                                                            gTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCACGCTGCTGGCGGTAGTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTtg                                       >  2:1302886/1‑150 (MQ=255)
                                                                                                                    ccTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGcga                               >  2:274394/1‑150 (MQ=255)
                                                                                                                     cTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCATAAAAGCCCGCACGTACTGCTCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGcgac                              <  1:1128212/150‑1 (MQ=255)
                                                                                                                     cTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGcgac                              >  1:1104834/1‑150 (MQ=255)
                                                                                                                     cTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGcgac                              <  1:2005536/150‑1 (MQ=255)
                                                                                                                     cTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGcgac                              >  2:1364256/1‑150 (MQ=255)
                                                                                                                          cTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCATTGTGGGGCGACGACTg                         >  1:1606042/1‑150 (MQ=255)
                                                                                                                            aCGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGAGGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGcc                       <  2:459040/150‑1 (MQ=255)
                                                                                                                             cGACCTCGGTACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGccc                      >  1:1888188/1‑150 (MQ=255)
                                                                                                                             cGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAAcc                                                                                                                       <  1:1484719/53‑1 (MQ=255)
                                                                                                                             cGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAAcc                                                                                                                       >  2:1484719/1‑53 (MQ=255)
                                                                                                                              gACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCt                     >  1:2042258/1‑150 (MQ=255)
                                                                                                                              gACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCt                     <  2:286835/150‑1 (MQ=255)
                                                                                                                                 cTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGAt                  >  2:231968/1‑150 (MQ=255)
                                                                                                                                 cTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGAt                  >  2:231998/1‑150 (MQ=255)
                                                                                                                                  tCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATa                 <  1:1364256/150‑1 (MQ=255)
                                                                                                                                  tCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATa                 >  2:1534880/1‑150 (MQ=255)
                                                                                                                                    gcgaCCCTGGCGGGTAACCTCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGc               >  1:2068693/4‑150 (MQ=255)
                                                                                                                                    gATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGc               <  2:2062900/150‑1 (MQ=255)
                                                                                                                                    gATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGc               >  2:2278715/1‑150 (MQ=255)
                                                                                                                                    gATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCACGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGc               >  2:1831099/1‑150 (MQ=255)
                                                                                                                                        cccTGGCGGGTAACCTCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACa           <  2:2068693/150‑1 (MQ=255)
                                                                                                                                        cccTGGAGGTTACCCGCCAGGGTCGCCTGTTGACGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACa           <  2:1772665/150‑1 (MQ=255)
                                                                                                                                        cccTGGAGGTTACCCGACAGGGTAGCCTGTTGAAGATCATCTAGGTCGAACCCTAGCGGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACa           <  2:1104834/150‑1 (MQ=255)
                                                                                                                                          cTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACAGc         >  2:1257800/1‑150 (MQ=255)
                                                                                                                                          cTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACAGc         >  1:872302/1‑150 (MQ=255)
                                                                                                                                           tGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACAGcc        >  1:885033/1‑150 (MQ=255)
                                                                                                                                                 ttACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACAGCCTGCGCa  >  1:921863/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCGGGGGCTGACGACTACCTGCTCAAGCCATTTGCCCTGTCGGAGCTGGCGGCGCGCATCGAAGCGGTGCTGCGCCGTGCCCAGGGCGGTGGGCGGCGCACCCTGCAGGTGGCTGACCTGAGCTACGACCTCGATACCCTGGAGGTTACCCGCCAGGGTCGCCTGTTGAAGCTCAACCCGGTCGGCCTCAAGCTGCTGGCGGTACTGATGCAGAAAAGCCCGCACGTACTGCGCCGCGAGGTGCTGGAAGAAGCCTTGTGGGGCGACGACTGCCCTGATAGCGACAGCCTGCGCA  >  NC_002947/1041800‑1042094

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: