Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,209,677 T→A 5.9% F301I (TTC→ATC) ‡ PP_2813 → BNR domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,209,6770TA 5.9% 209.3 / 2.6 68F301I (TTC→ATC) ‡PP_2813BNR domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base A (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

GGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGTTTTCTTGCCTGGGCGAATGGGTGGTCGT  >  NC_002947/3209536‑3209820
                                                                                                                                             |                                                                                                                                                
ggCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGc                                                                                                                                           <  7:134198/149‑1 (MQ=255)
 gCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTc                                                                                                                                                <  4:167957/143‑1 (MQ=255)
 gCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTc                                                                                                                                                >  3:167957/1‑143 (MQ=255)
  ctctGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGcc                                                                                                                                          >  6:124509/1‑148 (MQ=255)
  ctctGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGAACAACGCCg                                                                                                                                         <  5:271920/149‑1 (MQ=255)
     tGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAg                                                                                                                                      >  6:177563/1‑149 (MQ=255)
     tGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAg                                                                                                                                      >  4:200138/1‑149 (MQ=255)
                        gCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGtttt                                                                                                                    <  5:124509/148‑1 (MQ=255)
                             ctctCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGcc                                                                                                              <  5:3241/149‑1 (MQ=255)
                                  gtgtGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGAGGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAggg                                                                                                         <  4:203190/149‑1 (MQ=255)
                                   tgtgCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAgggg                                                                                                        <  5:153404/149‑1 (MQ=255)
                                           cGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGccc                                                                                               <  1:116836/149‑1 (MQ=255)
                                           cGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGccc                                                                                               >  2:550/1‑149 (MQ=255)
                                                gggCTGGATGGGCTTAGCGTGCAGAGCAGCGAGGGCGGGCACAGATTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCa                                                                                          <  4:8245/149‑1 (MQ=255)
                                                   cTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGa                                                                                       <  4:202667/149‑1 (MQ=255)
                                                     ggATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGt                                                                                     >  6:158812/1‑149 (MQ=255)
                                                      gATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTg                                                                                    <  1:550/149‑1 (MQ=255)
                                                               aGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTgg                                                                           >  8:101861/1‑149 (MQ=255)
                                                               aGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTgg                                                                           >  1:126038/1‑149 (MQ=255)
                                                               aGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCATTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTgg                                                                           >  7:156314/1‑149 (MQ=255)
                                                                gCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGc                                                                          <  5:158812/149‑1 (MQ=255)
                                                                 cGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGc                                                                          >  3:58999/1‑148 (MQ=255)
                                                                     cAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTc                                                                     <  7:101861/149‑1 (MQ=255)
                                                                        agcagcGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAg                                                                                 >  7:195653/1‑134 (MQ=255)
                                                                        agcagcGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAg                                                                                 <  8:195653/134‑1 (MQ=255)
                                                                         gcagcGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTc                                                                 <  4:58999/149‑1 (MQ=255)
                                                                             cGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTTCCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGaa                                                             <  5:177563/149‑1 (MQ=255)
                                                                             cGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAg                                                               >  5:74626/1‑147 (MQ=255)
                                                                             cGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAg                                                               >  7:92849/1‑147 (MQ=255)
                                                                             cGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGCCCTTTCAATGACGGCGGCGTTGTTCAACGCCGAACGTGTACCTGTATTGTTTTCACGCCGAGGAATAAAATACCAAAAACGAGTGAAGGGTTTGACACTCGCTCAg                                                               >  1:108082/1‑148 (MQ=255)
                                                                                  gCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCa                                                         <  6:51925/148‑1 (MQ=255)
                                                                                  gCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCa                                                         <  3:200138/148‑1 (MQ=255)
                                                                                  gCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAg                                                        <  6:74626/149‑1 (MQ=255)
                                                                                      gCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCAGAAGGTGTACCTGTATTGTATTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCAt                                                    >  7:207977/1‑149 (MQ=255)
                                                                                        acaGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAg                                                        <  4:50278/143‑1 (MQ=255)
                                                                                        acaGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAg                                                        >  3:50278/1‑143 (MQ=255)
                                                                                            tttCCAGGAAACGCAGACCGCGGATGACCTTTCACTGACGGCGGCGTTGTTAAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCaaa                                              >  5:327846/1‑149 (MQ=255)
                                                                                              tCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAAcc                                            >  6:157821/1‑149 (MQ=255)
                                                                                               ccAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAAcc                                            >  3:229365/1‑148 (MQ=255)
                                                                                                 aGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGAt                                         >  7:110048/1‑149 (MQ=255)
                                                                                                 aGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGAt                                         >  3:8645/1‑149 (MQ=255)
                                                                                                   gAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTc                                                                                                                   <  8:229460/74‑1 (MQ=255)
                                                                                                     aaCGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATgcct                                     <  4:229365/149‑1 (MQ=255)
                                                                                                        gCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATgcctggc                                  <  7:188576/149‑1 (MQ=255)
                                                                                                           gACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATgcctggcctg                               >  1:306582/1‑149 (MQ=255)
                                                                                                            aCCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGt                              >  7:43793/1‑149 (MQ=255)
                                                                                                            aCCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGt                              >  4:240267/1‑149 (MQ=255)
                                                                                                             ccGCGGATGGCCTTTCACTGAGTGCGGCGTTGTTCATGGCCGAATGTTTACCTGTATTGTTTTCAAGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGGCACTCTGTCAGAAGCCAGGCATCGCAAAGCGATGCCTGGCCTGtt                             <  2:108082/149‑1 (MQ=255)
                                                                                                             ccGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGtt                             <  8:92849/149‑1 (MQ=255)
                                                                                                                        cTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGtttt                           <  8:176846/140‑1 (MQ=255)
                                                                                                                        cTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGtttt                           >  7:176846/1‑140 (MQ=255)
                                                                                                                        cTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGTTTTCTTGCCTgg                  >  5:78381/1‑149 (MQ=255)
                                                                                                                             aCTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGAttc                                       >  4:189883/1‑121 (MQ=255)
                                                                                                                              cTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGTTTTCTTGCCTGGGCGAAt            <  6:327846/149‑1 (MQ=255)
                                                                                                                               tGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGtttt                           <  5:238039/133‑1 (MQ=255)
                                                                                                                               tGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGtttt                           >  6:238039/1‑133 (MQ=255)
                                                                                                                               tGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGTTTTCTTGCCTGGGCGAATg           >  3:88084/1‑149 (MQ=255)
                                                                                                                                accttcggcgTTGAACAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGTTTTCTTGCCt                    >  2:234335/6‑139 (MQ=255)
                                                                                                                                accttcggcgTTGAACAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGTTTTCTTGCCt                    <  1:234335/134‑1 (MQ=255)
                                                                                                                                 aCGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGTTTTCTTGCCTGGGCGAATgg          <  5:157821/148‑1 (MQ=255)
                                                                                                                                 aCGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGTTTTCTTGCCTGGGCGAATgg          <  8:110048/148‑1 (MQ=255)
                                                                                                                                 aCGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGTTTTCTTGCCTGGGCGAATgg          <  2:145592/148‑1 (MQ=255)
                                                                                                                                  cggcggCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGTTTTCTTGCCTGGGCGAATGggt        >  8:165395/1‑149 (MQ=255)
                                                                                                                                   ttcggcgTTGAACAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGTTTTCTTGCCTGGGCGAATGggtg       >  2:252002/3‑149 (MQ=255)
                                                                                                                                     cggcgTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGTTTTCTTGCCTGGGCGAATGggtggt     >  4:311796/1‑149 (MQ=255)
                                                                                                                                        cgTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAg                                                        >  5:107824/1‑95 (MQ=255)
                                                                                                                                        cgTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAg                                                        <  6:107824/95‑1 (MQ=255)
                                                                                                                                        cgTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGTTTTCTTGCCTGGGCGAATGGGTGGTCGt  <  2:306582/149‑1 (MQ=255)
                                                                                                                                            gttCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAAc                                             >  2:270107/1‑102 (MQ=255)
                                                                                                                                            gttCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAAc                                             <  1:270107/102‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
GGCTCTGTCACTTGCGTAGCGGTAGCCGGCTCTCGTGTGCTGGCGGTCGGGCTGGATGGGCTTAGCGTGCAGAGCAGCGATGGCGGGCACAGTTTCCAGGAAACGCAGACCGCGGATGGCCTTTCACTGACGGCGGCGTTGTTCAACGCCGAAGGTGTACCTGTATTGTTTTCACGCCGAGGGGTAATGCCCTACCAGGAGTGAAGGGTTGGCCACTCGCTCAGAAGCCAGGCATCGCAAACCGATGCCTGGCCTGTTTTCTTGCCTGGGCGAATGGGTGGTCGT  >  NC_002947/3209536‑3209820

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