Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,891,849 A→T 10.5% L12Q (CTG→CAG)  PP_2546 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,891,8490AT10.5% 226.9 / 15.6 76L12Q (CTG→CAG) PP_2546membrane protein
Reads supporting (aligned to +/- strand):  ref base A (31/37);  new base T (3/5);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

CCCAGGTGGTGAGCGGGCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGCCACCACCCGGCGCGTACGCT  >  NC_002947/2891703‑2891996
                                                                                                                                                  |                                                                                                                                                   
cccAGGTGGTGAGCGGGCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGTCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTgcagcagca                                                                                                                                                  >  1:1043937/1‑150 (MQ=255)
 ccAGGTGGTGAGCGGGCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCAc                                                                                                                                                 >  1:408429/1‑150 (MQ=255)
 ccAGGTGGTGAGCGGGCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCAc                                                                                                                                                 <  2:670383/150‑1 (MQ=255)
     gtggtgAGCGGGCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAgg                                                                                                                                             <  2:1866941/150‑1 (MQ=255)
      tggtgAGCGGGCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCTGCTGCACCTggc                                                                                                                                            <  2:1252253/150‑2 (MQ=255)
        gtgAGCGGGCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGgtg                                                                                                                                          >  1:13787/1‑150 (MQ=255)
              gggCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTcggcg                                                                                                                                    <  2:2214231/150‑1 (MQ=255)
                gCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCggcggc                                                                                                                                  >  1:1752467/1‑150 (MQ=255)
                   gCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAg                                                                                                                               <  2:391524/150‑1 (MQ=255)
                     ccAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGc                                                                                                                             >  2:805993/1‑150 (MQ=255)
                       aGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCgg                                                                                                                           >  2:2192736/1‑150 (MQ=255)
                         gTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCaa                                                                                                                                >  1:1534045/1‑143 (MQ=255)
                         gTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCaa                                                                                                                                <  2:1534045/143‑1 (MQ=255)
                          tGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTag                                                                                                                        >  1:353760/1‑150 (MQ=255)
                              ggCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTaga                                                                                                                       <  2:248580/147‑1 (MQ=255)
                              ggCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTaga                                                                                                                       >  1:248580/1‑147 (MQ=255)
                               gCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAggacggc                                                                                                                         >  1:1741990/1‑139 (MQ=255)
                               gCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAggacggc                                                                                                                         <  2:1741990/144‑6 (MQ=255)
                                cACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCATGTCGAACAGCAACCGGAAGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCAtg                                                                                                                  <  1:1500083/150‑1 (MQ=255)
                                 aCTCCACCGCCAGGCCCCAGTCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCAtgg                                                                                                                 <  2:1043937/150‑1 (MQ=255)
                                   tCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCCtgggt                                                                                                               >  2:869065/1‑150 (MQ=255)
                                      aCCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCAtgggttgg                                                                                                            <  1:2192736/150‑1 (MQ=255)
                                               ccccaggccAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCggg                                                                                                   >  2:1983541/1‑150 (MQ=255)
                                               ccccaggccAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCggg                                                                                                   <  2:2160760/150‑1 (MQ=255)
                                                cccaggccAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCgggg                                                                                                  >  1:1784375/1‑150 (MQ=255)
                                                 ccaggccaggACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGa                                                                                                 >  2:513919/1‑150 (MQ=255)
                                                   cggccaggCCCCGCAAGTCACGCAAACGGCCGCGGCGGCAGCGCAGCTGGATCATCCCGCCGAGCATTACCAGCATGAGGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGAtt                                                                                               >  2:1233374/2‑150 (MQ=255)
                                                      ccaggACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTggg                                                                                            <  1:1838102/150‑1 (MQ=255)
                                                      ccaggACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTggg                                                                                            <  1:2314611/150‑1 (MQ=255)
                                                      ccaggACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTggg                                                                                            <  1:1983541/150‑1 (MQ=255)
                                                       caggACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAg                                                                                                                                              <  2:90495/99‑1 (MQ=255)
                                                       caggACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAg                                                                                                                                              >  1:90495/1‑99 (MQ=255)
                                                        aggACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGt                                                                                          >  1:1527305/1‑150 (MQ=255)
                                                                 aaGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGc                                                                                 <  2:408429/150‑1 (MQ=255)
                                                                         caACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTgg                                                                         <  2:614663/150‑1 (MQ=255)
                                                                           aCAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCgg                                                                                                    <  1:1248127/121‑1 (MQ=255)
                                                                           aCAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCgg                                                                                                    >  2:1248127/1‑121 (MQ=255)
                                                                            cagggcagggCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCGTAGTTGGGtt                                                                      <  2:2127262/150‑1 (MQ=255)
                                                                                     gcgcCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGccggtcct                                                             <  2:1527305/150‑7 (MQ=255)
                                                                                      cgcCAGCGCAGCACCAGCAGCCCGCCGAACAGCAACAGGAGGAGGAAGGACAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTg                                                            <  2:193896/150‑1 (MQ=255)
                                                                                       gcCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTgg                                                           >  2:1806605/1‑150 (MQ=255)
                                                                                         cAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCTGGTGCTGCTGCACCTggc                                                                                                                                            <  2:454734/67‑2 (MQ=37)
                                                                                             gcagcaCCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTaca                                                     >  2:2385821/1‑150 (MQ=255)
                                                                                               agcaCCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGt                                                                                          <  2:1758700/111‑1 (MQ=255)
                                                                                               agcaCCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGt                                                                                          >  1:1758700/1‑111 (MQ=255)
                                                                                                 caCCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACAcccc                                                 <  1:1806605/150‑1 (MQ=255)
                                                                                                   ccAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCtt                                               >  2:479385/1‑150 (MQ=255)
                                                                                                      gcagcCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCt                                            >  1:462545/1‑150 (MQ=255)
                                                                                                          cccGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTAGAc                                        >  1:1809304/1‑150 (MQ=255)
                                                                                                          cccGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGACGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGAc                                        <  1:869065/150‑1 (MQ=255)
                                                                                                           ccgccgAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTacac                                                    >  2:1713583/1‑137 (MQ=255)
                                                                                                           ccgccgAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTacac                                                    <  1:1713583/137‑1 (MQ=255)
                                                                                                           ccgccgAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGAcc                                       <  2:353760/150‑1 (MQ=255)
                                                                                                            cgccgAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGAGTTGGGTGCCCGGGGGTTGGGGTGATCCCAGTAGAGGTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCa                                      >  2:880033/1‑150 (MQ=255)
                                                                                                              ccgAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATc                                    <  2:895733/150‑1 (MQ=255)
                                                                                                               cgAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCt                                   <  1:805993/150‑1 (MQ=255)
                                                                                                               cgAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCt                                   >  2:2403690/1‑150 (MQ=255)
                                                                                                                   cAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAg                               <  1:1156032/150‑1 (MQ=255)
                                                                                                                       agcagcagGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCGTCCTCGACCATCTGCAGGaaa                           <  2:564637/148‑1 (MQ=255)
                                                                                                                       aCCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGaaa                           <  2:1193290/150‑1 (MQ=255)
                                                                                                                       aCCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGaaa                           <  1:479385/150‑1 (MQ=255)
                                                                                                                       aCCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGaaa                           <  2:462545/150‑1 (MQ=255)
                                                                                                                       aCCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGaaa                           <  2:1752467/150‑1 (MQ=255)
                                                                                                                       aCCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGCTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGaaa                           <  2:71234/150‑1 (MQ=255)
                                                                                                                         cagcagGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCa                                                                                   >  2:590601/1‑92 (MQ=255)
                                                                                                                         cagcagGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCa                                                                                   <  1:590601/92‑1 (MQ=255)
                                                                                                                            cagGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGCGAATGAGATGGGTGCCTGGGGGTGGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGc                      >  1:2100830/1‑150 (MQ=255)
                                                                                                                             gttaTGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGcc                     >  1:1893158/4‑150 (MQ=255)
                                                                                                                               gATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGccac                   >  1:1742188/1‑150 (MQ=255)
                                                                                                                                  gAAGGCCAGGTGCAACAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGCcaccac                >  2:2022857/1‑150 (MQ=255)
                                                                                                                                       cctggTGCTGCTGCACCTGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGaaa                           <  1:716375/131‑1 (MQ=255)
                                                                                                                                       cctggTGCTGCTGCACCTGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGACTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGaaa                           >  2:716375/4‑134 (MQ=255)
                                                                                                                                       cctggTGCTGCTGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTggg                                                                                            <  1:1116397/66‑1 (MQ=255)
                                                                                                                                       cctggTGCTGCTGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTggg                                                                                            >  2:1116397/4‑69 (MQ=255)
                                                                                                                                       cctggTGCTGCTGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGc                      <  2:101565/136‑1 (MQ=255)
                                                                                                                                       cctggTGCTGCTGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGc                      >  1:101565/4‑139 (MQ=255)
                                                                                                                                       ccAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGCCACCACCCGgc           <  2:43618/150‑1 (MQ=255)
                                                                                                                                             gcagcagcaCCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGCCACCACCCGGCGCGTAc     >  1:1944653/1‑150 (MQ=255)
                                                                                                                                               agcggcaCCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCGCAGCATACTTGGGTGGCCTGATCTGGCCTACACGCCTTGCTCGACTATCTGCACGAAAGTCGCCACCACCCGGCGCGTACgc   >  2:2370233/1‑150 (MQ=255)
                                                                                                                                               agcagcaCCAGGGTGTCGGCGGCAGGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGCCACCACCCGGCGCGTACgc   <  1:880033/150‑1 (MQ=255)
                                                                                                                                               agcagcaCCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGCCACCACCCGGCGCGTACgc   >  1:949254/1‑150 (MQ=255)
                                                                                                                                                gcagcaCCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATTTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGCCACCACCCGGCGCGTACgct  >  1:407657/1‑150 (MQ=255)
                                                                                                                                                gcagcaCCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTAGGGTGCCCGGGGATTGGGATTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGCCACCACCCGGCGCGTACgct  >  2:935301/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCCAGGTGGTGAGCGGGCAGCCCAGGTGCAGGCACTCCACCGCCAGGCCCCAGGCCAGGACCGGCAAGTGCAGCAACAGGGCAGGGCGCCAGCGCAGCACCAGCAGCCCGCCGAACAGCACCAGCAGGATGAAGGCCAGGTGCAGCAGCACCAGGGTGTCGGCGGCAAGGCGGTAGAGCATGGGTTGGGTGCCCGGGGATTGGGGTTATCCCAGCATAGTTGGGTTGCCTGATCTGGCCTACACCCCTTCCTCGACCATCTGCAGGAAAGTCGCCACCACCCGGCGCGTACGCT  >  NC_002947/2891703‑2891996

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