Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,110,087 T→A 6.3% F622Y (TTC→TAC)  PP_3616 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,110,0870TA 6.3% 177.2 / 3.5 63F622Y (TTC→TAC) PP_3616membrane protein
Reads supporting (aligned to +/- strand):  ref base T (25/34);  new base A (2/2);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TTACGGAAAACCCCTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGG  >  NC_002947/4109943‑4110236
                                                                                                                                                |                                                                                                                                                     
ttACGGAAAACCCCTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTg                                                                                                                                                  <  1:1464283/150‑1 (MQ=32)
   cGGAAAACCCCTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCAGCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGt                                                                                                                                               <  2:889690/150‑1 (MQ=17)
         aCCCCTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCa                                                                                                                                         >  2:716503/1‑150 (MQ=31)
                   ccAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTAACCTCcggcgg                                                                                                                               >  1:230347/1‑150 (MQ=17)
                    cAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggt                                                                                                                              <  1:1311538/150‑1 (MQ=32)
                     aggaggTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtg                                                                                                                             <  1:716503/150‑1 (MQ=32)
                         ggTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtgg                                                                                                                         <  1:309653/150‑1 (MQ=32)
                          gTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggt                                                                                                                        <  1:653167/150‑1 (MQ=32)
                           tGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGAGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggtg                                                                                                                       <  2:490092/150‑1 (MQ=17)
                             gatgttgaTTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggtggt                                                                                                                     <  2:588995/143‑1 (MQ=17)
                             cAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggtggt                                                                                                                     <  1:490241/150‑1 (MQ=32)
                                cGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCAATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTa                                                                                                                  >  1:148504/1‑150 (MQ=17)
                                 ggTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTag                                                                                                                 <  2:1509598/150‑2 (MQ=17)
                                  gTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggtg                                                                                                                       >  1:907485/1‑143 (MQ=32)
                                  gTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggtg                                                                                                                       <  2:907485/143‑1 (MQ=32)
                                       tGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCgg                                                                                                           <  1:1483740/150‑1 (MQ=32)
                                       tGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCgg                                                                                                           <  2:230347/150‑1 (MQ=32)
                                       tGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCgg                                                                                                           <  1:516868/150‑1 (MQ=32)
                                        gAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCggg                                                                                                          >  1:771376/1‑150 (MQ=32)
                                               ggCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCtg                                                                                                   <  1:26214/150‑1 (MQ=32)
                                                           gTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggt                                                                                                                        <  1:6172/117‑1 (MQ=33)
                                                           gTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCggtggtggt                                                                                                                        >  2:6172/1‑117 (MQ=33)
                                                           gTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATc                                                                                       <  2:941859/150‑1 (MQ=32)
                                                                agagATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATc                                                                                  <  1:720723/150‑1 (MQ=32)
                                                                   gaTGAAAGCGTTGCGCATCTCGGGTAGCCGCATAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCt                                                                               >  1:279165/1‑150 (MQ=17)
                                                                    aTGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTaaa                                                                                                                >  1:120067/1‑116 (MQ=33)
                                                                    aTGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTaaa                                                                                                                <  2:120067/116‑1 (MQ=33)
                                                                    aTGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTc                                                                              >  1:18844/1‑150 (MQ=32)
                                                                             ttGCGCATCTCGGGTAGCCGCATAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGcc                                                                     <  2:279165/150‑1 (MQ=17)
                                                                               gcgcATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTa                                                                   <  2:18844/150‑1 (MQ=32)
                                                                                   aTCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTacaac                                                               >  1:1040016/1‑150 (MQ=32)
                                                                                    tctcGGATAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTacaaca                                                              >  2:602404/1‑150 (MQ=17)
                                                                                         ggTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCgg                                                                                                           >  2:284784/1‑100 (MQ=34)
                                                                                         ggTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCgg                                                                                                           <  1:284784/100‑1 (MQ=34)
                                                                                                gCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCagaaga                                                  >  1:1268128/1‑150 (MQ=35)
                                                                                                gCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCagaaga                                                  <  2:1040016/150‑1 (MQ=35)
                                                                                                  acaGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGc                                                >  1:1458902/1‑150 (MQ=35)
                                                                                                         ccATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGTTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCgg                                                                                                           <  1:1330393/84‑1 (MQ=25)
                                                                                                         ccATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCgg                                                                                                           >  2:1330393/1‑84 (MQ=35)
                                                                                                         ccATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCAc                                         >  1:1612695/1‑150 (MQ=35)
                                                                                                          cATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCAcc                                        >  1:672044/1‑150 (MQ=35)
                                                                                                           aTCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCg                                       <  2:1458902/150‑1 (MQ=35)
                                                                                                                  tGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGc                                                            <  1:156551/122‑1 (MQ=33)
                                                                                                                  tGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGc                                                            >  2:156551/1‑122 (MQ=32)
                                                                                                                    accaccACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACGGCCTGCAGAAGAGCGATGCACCGGAATACAcc                              <  2:1268128/150‑1 (MQ=25)
                                                                                                                     ccaccaCTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACAccc                             >  2:59465/1‑150 (MQ=35)
                                                                                                                     ccaccaCTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCGCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGATCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACAccc                             <  1:1446760/150‑1 (MQ=21)
                                                                                                                      caccgccggaggtgagtAACCCAAGGAAGTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGc                                                        >  2:1455195‑M1/18‑122 (MQ=255)
                                                                                                                      caccgccggaggtgagtAACCCAAGGAAGTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGc                                                        <  1:1455195‑M1/105‑1 (MQ=255)
                                                                                                                       accgccggaggtgagtAACCCAAGGAAGTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGa                                                                          <  1:1609045‑M1/87‑1 (MQ=255)
                                                                                                                       accgccggaggtgagtAACCCAAGGAAGTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGa                                                                          >  2:1609045‑M1/17‑103 (MQ=255)
                                                                                                                       accacTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTg                           >  2:1577703/1‑150 (MQ=35)
                                                                                                                               gcgTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGc                   >  1:1386341/1‑150 (MQ=35)
                                                                                                                                 gTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCtt                 >  1:1058588/1‑150 (MQ=35)
                                                                                                                                     gAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTac                                                                  >  1:1465648/1‑97 (MQ=33)
                                                                                                                                     gAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTac                                                                  <  2:1465648/97‑1 (MQ=34)
                                                                                                                                        ACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCtg                                                                                                   <  2:272975‑M1/61‑1 (MQ=255)
                                                                                                                                        ACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCtg                                                                                                   >  1:272975‑M1/1‑61 (MQ=255)
                                                                                                                                          CCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATc                                                                                  >  2:891902‑M1/1‑76 (MQ=255)
                                                                                                                                          CCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATc                                                                                  <  1:891902‑M1/76‑1 (MQ=255)
                                                                                                                                           CAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGtt       <  1:1450306‑M1/150‑1 (MQ=255)
                                                                                                                                             ACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCg     <  2:487406‑M1/150‑1 (MQ=255)
                                                                                                                                              cTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCTTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCgg    <  1:34593/150‑1 (MQ=38)
                                                                                                                                               ttCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCggg   <  1:1577703/150‑1 (MQ=38)
                                                                                                                                               ttCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCggg   <  2:771376/150‑1 (MQ=38)
                                                                                                                                                tCCTTGGGTTACGCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATCCACCGGAATACACCCTGGGGTTCTCTTCATCTGTGTTCgggg  <  2:148504/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TTACGGAAAACCCCTTTAGCCAGGAGGTGCAGCGGTATTTGAAGAGCGGCATGGCGCAGGTTGAAGAGATGAAAGCGTTGCGCATCTCGGGTAGCCGCACAGTGTCCATCGAGCTGACCACCACTGCGCGTGCGAAACCCAACTTCCTTGGGTTACTCACCTCCGGCGGTGGTGGTGGTTTAAACGCGGGGATGCTGTGGTTCAACATCTTATCGCTCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGG  >  NC_002947/4109943‑4110236

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: