Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,408,406:1 +T 6.2% coding (202/1566 nt) PP_2111 ← aerotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,408,4061.T 6.2% 383.7 / 18.6 97coding (202/1566 nt)PP_2111aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base . (44/47);  new base T (3/3);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-01

CCGCACCGACTCGAAGCCGACAACCTGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTGGAGATCAACCGC  >  NC_002947/2408263‑2408554
                                                                                                                                                |                                                                                                                                                    
ccGCACCGACTCGAAGCCGACAACCTGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCaaa                                                                                                                                                >  2:2096642/1‑150 (MQ=255)
ccGCACCGACTCGAAGCCGACAACCTGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACCATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCaaa                                                                                                                                                >  2:2096702/1‑150 (MQ=255)
      cGACTCGAAGCCGACAACCTGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGc                                                                                                                                          <  1:2134777/150‑1 (MQ=255)
         cTCGAAGCCGACAACCTGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAgg                                                                                                                                       >  2:2100163/1‑150 (MQ=255)
              aGCCGACAACCTGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTa                                                                                                                                  <  2:1937085/150‑1 (MQ=255)
                       ccTGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGCGGCCTATCCACGGCTGGCCCTGCTTGAGGGCTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGt                                                                                                                         >  1:2223018/1‑150 (MQ=255)
                       ccTGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGt                                                                                                                         >  2:2353610/1‑150 (MQ=255)
                        cTGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTg                                                                                                                        <  2:177638/150‑1 (MQ=255)
                         tGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTgg                                                                                                                       >  2:2297249/1‑150 (MQ=255)
                            ttgttATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACgg                                                                                                                                           <  2:1990574/127‑1 (MQ=255)
                            ttgttATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACgg                                                                                                                                           >  1:1990574/1‑127 (MQ=255)
                            ttgttATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGgcgc                                                                                                                    <  1:234349/150‑1 (MQ=255)
                             tgttATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCa                                                                                                                   <  1:1558151/150‑1 (MQ=255)
                               gcaTCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCAcc                                                                                                                 >  1:1892729/3‑150 (MQ=255)
                                 aTCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCggg                                                                                                                           >  2:1194110/1‑138 (MQ=255)
                                 aTCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCggg                                                                                                                           <  1:1194110/138‑1 (MQ=255)
                                  tCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAgg                                                                                                              >  1:1675677/1‑150 (MQ=255)
                                    gAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGtt                                                                                                            >  2:2062979/1‑150 (MQ=255)
                                    gAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGtt                                                                                                            <  2:539372/150‑1 (MQ=255)
                                       aaTGGGGGTGACATAGGCGTTGACCCCGTAATGGTAGCATGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTCAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTtgtg                                                                                                         <  1:1705995/150‑1 (MQ=37)
                                        aTGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTgtgc                                                                                                        <  1:1389428/150‑1 (MQ=255)
                                           ggggTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTTCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTgtgcggt                                                                                                     <  1:2100163/150‑1 (MQ=255)
                                             ggTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACg                                                                                                                                            >  2:1314317/1‑109 (MQ=255)
                                             ggTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACg                                                                                                                                            <  1:1314317/109‑1 (MQ=255)
                                             ggTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTAGGGGTGGCGCACCAGGTTgtgcggtgc                                                                                                   >  1:1103748/1‑150 (MQ=255)
                                              gTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATgg                                                                                                                                                     <  1:435524/99‑1 (MQ=255)
                                              gTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATgg                                                                                                                                                     >  2:435524/1‑99 (MQ=255)
                                               tGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTgcgc                                                                                                 <  1:2353610/150‑1 (MQ=255)
                                                   aTAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAac                                                                                                                                               <  2:2229029/100‑1 (MQ=255)
                                                   aTAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAac                                                                                                                                               >  1:2229029/1‑100 (MQ=255)
                                                   aTAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGt                                                                                             >  2:951918/1‑150 (MQ=255)
                                                   aTAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGCCAGTCACGGCATGCGGTACGTCGGGGTGGCGGACCAGGTTGTGCGGTAAGCCCGt                                                                                             >  1:1353276/1‑150 (MQ=255)
                                                     aGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCa                                                                                           >  1:1084557/1‑150 (MQ=255)
                                                       gCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGc                                                                                         >  1:2303537/1‑150 (MQ=255)
                                                       gCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGc                                                                                         >  1:1848935/1‑150 (MQ=255)
                                                              cccAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGt                                                                                                                                   >  2:2343796/1‑101 (MQ=255)
                                                              cccAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGt                                                                                                                                   <  1:2343796/101‑1 (MQ=255)
                                                               ccAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTgtgcggt                                                                                                     <  2:303268/130‑1 (MQ=255)
                                                               ccAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTgtgcggt                                                                                                     >  1:303268/1‑130 (MQ=255)
                                                               ccAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTcgcg                                                                                 <  2:2223018/150‑1 (MQ=255)
                                                               ccAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGCGTGGCGGAACAGGTTGTGCGGTGCGCCCGTCAGCTCCTcgcg                                                                                 >  1:2354412/1‑150 (MQ=255)
                                                                cAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCgg                                                                                <  1:2096642/150‑1 (MQ=255)
                                                                cAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCgg                                                                                <  2:2274917/150‑1 (MQ=255)
                                                                cAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACCATACCCATCCACGGCTGGCCCTGCTTGCGGGCTTGCCACATACGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCgg                                                                                <  1:2096702/150‑1 (MQ=255)
                                                                cAGTAATGGTCGCCGGACCTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCgg                                                                                <  1:1898507/150‑1 (MQ=255)
                                                                  gTAATGGTCGCCGGACTTGCAGCGGTTCTTGTCGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGt                                                                              <  1:749729/150‑1 (MQ=255)
                                                                  gTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGt                                                                              >  1:1445910/1‑150 (MQ=255)
                                                                      tGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGaaa                                                                          <  1:1194731/150‑1 (MQ=255)
                                                                        gTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAAcc                                                                        >  1:882137/1‑150 (MQ=255)
                                                                          cGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTTCGCCCGTCAGCTCCTCGCGGGTGAAACCAc                                                                      <  1:316768/150‑1 (MQ=255)
                                                                           gCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACt                                                                     <  2:1848935/150‑1 (MQ=255)
                                                                           gCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACt                                                                     <  1:951918/150‑1 (MQ=255)
                                                                           gCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACt                                                                     >  2:2092500/1‑150 (MQ=255)
                                                                                 cTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTc                                                               >  1:1038117/1‑150 (MQ=255)
                                                                                   tGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGa                                                             <  1:2092500/150‑1 (MQ=255)
                                                                                    gCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGAt                                                            >  2:97694/1‑150 (MQ=255)
                                                                                    gCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGAt                                                            >  1:609547/1‑150 (MQ=255)
                                                                                      aGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATaa                                                          >  1:764704/1‑150 (MQ=255)
                                                                                         ggTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAgg                                                       <  2:1852044/150‑1 (MQ=255)
                                                                                            tCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCAt                                                    <  2:1084557/150‑1 (MQ=255)
                                                                                            tCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCAt                                                    <  2:1675677/150‑1 (MQ=255)
                                                                                             cTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATc                                                   <  2:764704/150‑1 (MQ=255)
                                                                                              ttGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTgg                                                                                                                       <  1:2256815/81‑1 (MQ=255)
                                                                                              ttGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTgg                                                                                                                       >  2:2256815/1‑81 (MQ=255)
                                                                                              ttGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTc                                                                                            >  1:1013864/1‑108 (MQ=255)
                                                                                              ttGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTc                                                                                            <  2:1013864/108‑1 (MQ=255)
                                                                                              ttGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCg                                                  <  2:1445910/150‑1 (MQ=255)
                                                                                              ttGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGCGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCg                                                  >  1:1042993/1‑150 (MQ=255)
                                                                                                gACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGtt                                                >  2:1521617/1‑150 (MQ=255)
                                                                                                       cccATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTg                                               <  1:2037749/144‑1 (MQ=255)
                                                                                                       cccATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTg                                               >  2:2037749/1‑144 (MQ=255)
                                                                                                        ccATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCGGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAgg                                        <  2:1892729/150‑1 (MQ=255)
                                                                                                        ccATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAgg                                        >  1:1734055/1‑150 (MQ=255)
                                                                                                         cATCGACGTCTAGCCTTGCTTTAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTGGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGt                                       <  2:882137/150‑1 (MQ=38)
                                                                                                         cATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGccc                                                                                               <  2:1426474/94‑1 (MQ=255)
                                                                                                         cATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGccc                                                                                               >  1:1426474/1‑94 (MQ=255)
                                                                                                         cATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCgg                                                                                >  2:1344727/1‑109 (MQ=255)
                                                                                                         cATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCgg                                                                                <  1:1344727/109‑1 (MQ=255)
                                                                                                         cATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGt                                       >  1:2153162/1‑150 (MQ=255)
                                                                                                          aTCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTg                                      >  2:1046900/1‑150 (MQ=255)
                                                                                                           tCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGa                                     <  2:609547/150‑1 (MQ=255)
                                                                                                             cACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATc                                   <  2:1038117/150‑1 (MQ=255)
                                                                                                                            ttGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTg                    <  1:722556/150‑1 (MQ=255)
                                                                                                                              ccgtgTTTGCCCATATGTGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGt                  >  1:2333280/5‑150 (MQ=255)
                                                                                                                                 tgTTTGCCCATATGTGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCgg                                                                                <  1:1690446/84‑1 (MQ=38)
                                                                                                                                 tgTTTGCCCATATGTGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCgg                                                                                >  2:1690446/2‑85 (MQ=37)
                                                                                                                                 tgTTTGCCCATATGTGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAg                                                        >  2:425126/2‑109 (MQ=39)
                                                                                                                                 tgTTTGCCCATATGTGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAg                                                        <  1:425126/108‑1 (MQ=39)
                                                                                                                                  gTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTg              <  2:2303537/150‑1 (MQ=255)
                                                                                                                                   tttGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGCGATCACACCCTTGGCGTTGGTAGTgg             <  2:2354412/150‑1 (MQ=255)
                                                                                                                                    ttGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTGga            >  1:2031433/1‑150 (MQ=255)
                                                                                                                                     tGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTGgag           <  1:2062979/150‑1 (MQ=255)
                                                                                                                                       gcccaTATGTGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTGgaga          <  2:2333280/147‑1 (MQ=255)
                                                                                                                                       ccACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATaaa                                                         >  2:950136/1‑102 (MQ=255)
                                                                                                                                       ccACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATaaa                                                         <  1:950136/102‑1 (MQ=255)
                                                                                                                                       ccACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTGGAGAt         >  2:1688362/1‑150 (MQ=255)
                                                                                                                                       ccACATATGGGCAAACACGGCAGCCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGCTAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTGGAGAt         >  2:993793/1‑150 (MQ=255)
                                                                                                                                        cacaTATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTCGCATTAGGTGATCACACCCTTGGCGTTGGTAGTGGAGATc        <  1:2297249/150‑1 (MQ=255)
                                                                                                                                              tGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTGGAGATCAACCgc  >  2:1913310/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCGCACCGACTCGAAGCCGACAACCTGGTTGTTATCGAAAATGGGGGTGACATAGGCGTTGACCCAGTAATGGTCGCCGGACTTGCAGCGGTTCTTGACGATACCCATCCACGGCTGGCCCTGCTTGAGGGTTTGCCACATATGGGCAAACACGGCAGGCGGTACGTCGGGGTGGCGCACCAGGTTGTGCGGTGCGCCCGTCAGCTCCTCGCGGGTGAAACCACTGATTTCGATAAAGGCATCGTTGCAGTAGGTGATCACACCCTTGGCGTTGGTAGTGGAGATCAACCGC  >  NC_002947/2408263‑2408554

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