Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,118,720 T→A 8.1% L799Q (CTG→CAG)  PP_1884 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,118,7200TA 8.1% 211.9 / 8.3 74L799Q (CTG→CAG) PP_1884hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/32);  new base A (4/2);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

GGCAGATCGATATCGATGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTACCAATT  >  NC_002947/2118579‑2118862
                                                                                                                                             |                                                                                                                                              
ggCAGATCGATATCGATGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGccgctgccgctg                                                                                                                                         <  8:129884/149‑1 (MQ=255)
 gCAGATCGATATCGTTGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCATACCAGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTgg                                                                                                                                        >  1:304127/1‑149 (MQ=255)
  cAGATCGATATCGATGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTggg                                                                                                                                       <  3:40890/149‑1 (MQ=255)
       cgatatccatGAAGACGGCAAGCCTTACTACTAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTAACCCATCACGCTCGAGGGTTTCATCGTGCTGCAGCTGCCGCTGGGgcag                                                                                                                                   >  6:268642/1‑148 (MQ=255)
        gatatggatGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGAATTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGgcagca                                                                                                                                 <  6:280754/149‑1 (MQ=255)
          tatcgatGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGgcagca                                                                                                                                 >  6:243912/1‑147 (MQ=255)
          tatcgatGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGcaca                                                                                                                               <  3:58466/149‑1 (MQ=255)
           atcgatGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAg                                                                                                                              >  4:347597/1‑149 (MQ=255)
              gatGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTc                                                                                                                            <  4:67480/148‑1 (MQ=255)
                  aaGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGc                                                                                                                        <  4:51172/148‑1 (MQ=255)
                   aGACGGCAATCCTTACTAGCAGGCCGACGCCGCGAGGCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTgg                                                                                                                                        <  4:224481/131‑1 (MQ=255)
                   aGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTgg                                                                                                                                        >  3:224481/1‑131 (MQ=255)
                   aGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCt                                                                                                                      <  7:31107/149‑1 (MQ=255)
                            gCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAG‑CTTTGAGGCCCt                                                                                                            >  3:277227/1‑149 (MQ=255)
                               ttACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGcccgtt                                                                                                          >  4:245700/1‑146 (MQ=255)
                                 actacCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGa                                                                                                         <  3:96373/148‑1 (MQ=255)
                                  ctacCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAg                                                                                                        <  3:123583/148‑1 (MQ=255)
                                   tacCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCaa                                                                                                                          <  3:107414/129‑1 (MQ=255)
                                   tacCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCaa                                                                                                                          >  4:107414/1‑129 (MQ=255)
                                        ggCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAgcaggcag                                                                                                 >  4:229615/1‑149 (MQ=255)
                                            gACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAgcgc                                                                                              >  3:355844/1‑148 (MQ=255)
                                               gccgcGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACg                                                                                          <  7:167250/149‑1 (MQ=255)
                                                  gcgAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACggg                                                                                        >  5:375736/1‑148 (MQ=255)
                                                         ccGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGcac                                                                                                                                <  7:250736/101‑1 (MQ=255)
                                                         ccGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGcac                                                                                                                                >  8:250736/1‑101 (MQ=255)
                                                         ccGCACATTTCTCCCTGGTTCGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCAACACGCTCGAGGCTTTCATCGTGCTGCCGCTGCCGCAGGGGCAGCACAGTCAACCCTATCAGGCCCTAGACCACCCAGCGCAACCCGCAAcgcg                                                                                 >  2:107192/1‑148 (MQ=255)
                                                          cGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAAcgcg                                                                                 >  1:69048/1‑147 (MQ=255)
                                                             acaTTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAgg                                                                             >  5:168292/1‑148 (MQ=255)
                                                                tttCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGAAGGCAGAGCAACGGGCAACGCGCAGGTAGt                                                                         >  8:249843/1‑149 (MQ=255)
                                                                      cGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCg                                                                   >  6:159006/1‑149 (MQ=255)
                                                                           ttggtGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGa                                                                                                         >  5:178462/1‑106 (MQ=255)
                                                                           ttggtGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGa                                                                                                         <  6:178462/106‑1 (MQ=255)
                                                                              gtGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAg                                                                              >  5:207242/1‑130 (MQ=255)
                                                                              gtGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAg                                                                              <  6:207242/130‑1 (MQ=255)
                                                                              gtGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCa                                                            <  4:355844/148‑1 (MQ=255)
                                                                              gtGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCAc                                                           >  5:28165/1‑149 (MQ=255)
                                                                                 cGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCt                                                         >  1:284302/1‑148 (MQ=255)
                                                                                   gTCTTGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTaaa                                                      <  3:229615/149‑1 (MQ=255)
                                                                                   gTCTGTCGGCTCTCCCCCTACTAACCAATCCCGCGCTAGGGTTTAATCGTCCTGCCCCTGCCGCTGGGGAAGAAAAGTCATG‑CTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAAc                                                     <  4:277227/149‑1 (MQ=255)
                                                                                      tGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAAcc                                                    >  1:346125/1‑148 (MQ=255)
                                                                                        gAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAAcc                                                    <  3:347597/146‑1 (MQ=255)
                                                                                         aGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATa                                                <  5:243912/149‑1 (MQ=255)
                                                                                             tctcCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAAcgc                                                                                  <  4:53947/111‑1 (MQ=255)
                                                                                             tctcCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAAcgc                                                                                  >  3:53947/1‑111 (MQ=255)
                                                                                             tctcCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCtt                                            >  5:5931/1‑149 (MQ=255)
                                                                                                 ccccGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAAc                                                     >  1:237590/1‑136 (MQ=255)
                                                                                                 ccccGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAAc                                                     <  2:237590/136‑1 (MQ=255)
                                                                                                 ccccGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTc                                          <  2:69048/147‑1 (MQ=255)
                                                                                                 ccccGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAg                                        >  3:243100/1‑149 (MQ=255)
                                                                                                  cccGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAg                                        >  3:129193/1‑148 (MQ=255)
                                                                                                    cGACTTACCCATGTCGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCAc                                     <  5:268642/149‑1 (MQ=255)
                                                                                                    cGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCAc                                     >  5:226104/1‑149 (MQ=255)
                                                                                                           ccccTCACGCTCGAGGGGTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGtt                              <  4:20435/149‑1 (MQ=255)
                                                                                                           cccATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTGAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAAGGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGt                               <  1:107192/148‑1 (MQ=255)
                                                                                                           cccATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGt                               <  6:5931/148‑1 (MQ=255)
                                                                                                              aTCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTGAGCACTCCAGTTGAg                           >  8:332664/1‑149 (MQ=255)
                                                                                                              aTCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCCATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCCTACCTTGAAGCACTCCAGTTGAg                           >  8:81586/1‑149 (MQ=255)
                                                                                                                  cGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATa                                                >  5:46872/1‑124 (MQ=255)
                                                                                                                  cGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATa                                                <  6:46872/124‑1 (MQ=255)
                                                                                                                  cGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCggg                       <  7:249843/149‑1 (MQ=255)
                                                                                                                    cTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGc                     <  6:28165/149‑1 (MQ=255)
                                                                                                                     tCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCa                                         >  3:350847/1‑128 (MQ=255)
                                                                                                                     tCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCa                                         <  4:350847/128‑1 (MQ=255)
                                                                                                                         gggTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCAcat                >  5:242413/1‑149 (MQ=255)
                                                                                                                            tttCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCGTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCAcatcat             >  8:9546/1‑149 (MQ=255)
                                                                                                                            tttCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGGAGTGAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCAcatcat             <  6:375736/149‑1 (MQ=255)
                                                                                                                             ttCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATg            >  2:276705/1‑149 (MQ=255)
                                                                                                                               cATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGc           <  5:159006/148‑1 (MQ=255)
                                                                                                                                aTCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGcc          >  7:139913/1‑148 (MQ=255)
                                                                                                                                  tgTGCTGCCCCAGTGGCAGCGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTAc       >  2:345190/2‑149 (MQ=255)
                                                                                                                                  tgTGCTGCCCCAGCGGCAGCGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCa                                                            <  7:82339/95‑1 (MQ=38)
                                                                                                                                  tgTGCTGCCCCAGCGGCAGCGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCa                                                            >  8:82339/2‑96 (MQ=38)
                                                                                                                                   gTGCTGCCCCAGCGGCAGCGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTa        <  5:52441/147‑1 (MQ=255)
                                                                                                                                   gTGCTGCCCCAGCGGCAGCGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTa        >  6:52441/1‑147 (MQ=255)
                                                                                                                                    tgctgcCCCAGCGGCAGCGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTAcc      >  5:372517/1‑148 (MQ=255)
                                                                                                                                      cttccgctgccgctgGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTACCaa    <  4:243100/145‑1 (MQ=255)
                                                                                                                                      ctgccgctgccgCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCATGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTACCaa    <  6:226104/148‑1 (MQ=255)
                                                                                                                                      ctgccgctgccgCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTACCAAt   <  2:346125/149‑1 (MQ=255)
                                                                                                                                       tgccgctgccgcTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTACCAAt   >  3:53331/1‑148 (MQ=255)
                                                                                                                                        tgtgctgccccagCGGCAGCACAGTCAAGCCTATCAGGCGCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCAc                                     >  5:78172/4‑113 (MQ=255)
                                                                                                                                        tgtgctgccccagCGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCAc                                     <  6:78172/110‑1 (MQ=255)
                                                                                                                                        gccgctgccgctGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTACCAAtt  >  4:31236/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                              
GGCAGATCGATATCGATGAAGACGGCAAGCCTTACTACCAGGCCGACGCCGCGAGCCCCGCACATTTCTGCGTGGTTGGTGCGGTCTGGAGGCTCTCCCCCGACTTACCCATCACGCTCGAGGGTTTCATCGTGCTGCCGCTGCCGCTGGGGCAGCACAGTCAAGCCTATCAGGCCCTTGAGCAGGCAGCGCAACGGGCAACGCGCAGGTAGTCAACCGGCGTTCACCTAAACCCATACCTTTCAGCACTCCAGTTGAGCGGGGCCACATCATGCCTACCAATT  >  NC_002947/2118579‑2118862

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