Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 74,755 T→G 7.1% K24Q (AAG→CAG)  PP_0064 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294774,7550TG 7.1% 280.5 / 7.4 91K24Q (AAG→CAG) PP_0064hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (40/44);  new base G (5/2);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGCTTTGCGCGCCGCCGCCAGCTCGCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATCCACAAC  >  NC_002947/74607‑74901
                                                                                                                                                    |                                                                                                                                                  
cACGCTTTGCGCGCCGCCGCCAGCTCGCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCAAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCtt                                                                                                                                                    <  1:227337/149‑1 (MQ=255)
   gCTTTGCGCGCCGCCGCCAGCTCGCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGcc                                                                                                                                                 <  1:250735/149‑1 (MQ=255)
   gCTTTGCGCGCCGCCGCCAGCTCGCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGcc                                                                                                                                                 <  5:347577/149‑1 (MQ=255)
     tttGCGCGCCGCCGCCAGCTCGCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGcc                                                                                                                                                 >  6:377605/1‑147 (MQ=255)
      ttGCGCGCCGCCGCCAGCTCGCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAg                                                                                                                                              >  8:169328/1‑149 (MQ=255)
        gcgcgcCGCCGCCAGCTCGCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAgg                                                                                                                                             <  5:33364/148‑1 (MQ=255)
         cgcgcCGCCGCCAGCTCGCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGc                                                                                                                                            <  1:81989/148‑1 (MQ=255)
          gcgcCGCCGCCAGCTCGCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGcc                                                                                                                                           >  3:393896/1‑148 (MQ=255)
             ccgccgccAGCTCGCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAAc                                                                                                                                       <  7:169328/149‑1 (MQ=255)
                   ccAGCTCGCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTACCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTcag                                                                                                                                  >  5:284799/1‑148 (MQ=255)
                    cAGCTCGCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTcag                                                                                                                                  >  3:175581/1‑147 (MQ=255)
                        tCGCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCgg                                                                                                                            >  1:280482/1‑149 (MQ=255)
                          gCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCggg                                                                                                                           <  8:12426/148‑1 (MQ=255)
                          gCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCggg                                                                                                                           <  1:318287/148‑1 (MQ=255)
                          gCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCggg                                                                                                                           <  5:354142/148‑1 (MQ=255)
                          gCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCggg                                                                                                                           <  3:408655/148‑1 (MQ=255)
                          gCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCgg                                                                                                                            <  3:132362/147‑1 (MQ=255)
                          gCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGt                                                                                                                          >  5:209673/1‑149 (MQ=255)
                            ccTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTa                                                                                                                        <  5:49422/149‑1 (MQ=255)
                            ccTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTa                                                                                                                        >  6:304323/1‑149 (MQ=255)
                                cAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATc                                                                                                                      <  5:292806/147‑1 (MQ=255)
                                cAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATc                                                                                                                      <  4:393896/147‑1 (MQ=255)
                                cAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATc                                                                                                                      <  2:285617/147‑1 (MQ=255)
                                 aGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCca                                                                                                                    <  7:311918/148‑1 (MQ=255)
                                     gATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCt                                                                                                               >  3:219243/1‑149 (MQ=255)
                                        aGGCCTTACCCAGCAGCTTGCATGCTGGCGATTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCttt                                                                                                             <  7:390663/148‑1 (MQ=255)
                                         ggCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCCGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTg                                                                                                            <  5:266375/148‑1 (MQ=255)
                                           ccTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGcc                                                                                                         <  6:386546/149‑1 (MQ=255)
                                            cTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCa                                                                                                        <  4:175581/149‑1 (MQ=255)
                                            cTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCa                                                                                                        <  7:396870/149‑1 (MQ=255)
                                               aCCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGc                                                                                                      <  5:377605/148‑1 (MQ=255)
                                               aCCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCa                                                                                                     >  7:154788/1‑149 (MQ=255)
                                                 ccAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAgaa                                                                                                     >  6:135534/1‑145 (MQ=255)
                                                  cagcagCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACATAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCAtct                                                                                                  >  6:147225/1‑149 (MQ=255)
                                                   agcagcTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCAtctt                                                                                                 >  1:339823/1‑149 (MQ=255)
                                                    gcagcTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCAtcttc                                                                                                >  4:77379/1‑149 (MQ=255)
                                                      agcTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTtctc                                                                                              <  4:10770/149‑1 (MQ=255)
                                                           ccATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGc                                                                                          >  3:408785/1‑148 (MQ=255)
                                                           ccATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGc                                                                                          >  3:415322/1‑148 (MQ=255)
                                                            cATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGa                                                                                        >  6:80376/1‑149 (MQ=255)
                                                              tgctgcCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCTAGCATCTTCTCCAGCGAtt                                                                                      >  5:161356/1‑149 (MQ=255)
                                                              tgctgcCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGAtt                                                                                      <  5:304323/149‑1 (MQ=255)
                                                                ctgcCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTc                                                                                     <  8:154788/148‑1 (MQ=255)
                                                                   ctgAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACCCCTTTGGCCAGCATCTTCTCCAGCGATTcgcg                                                                                  >  6:230440/3‑148 (MQ=255)
                                                                   ccGAGTATTTCTGATCCTGCTCTACACAGCGTGCCAGGTTCACTGGTGCTTCAGCGCCGTTGCCTTCATCCAGCCATCCCTTGCCCATTCCAATCCTCAGCAGCGGGTTATCCACACCTTTGGACAGCATCTTCTCCAGCGATTcgctc                                                                                 <  7:303574/149‑3 (MQ=255)
                                                                   ccGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTcgcg                                                                                  >  5:429362/1‑148 (MQ=255)
                                                                     gAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCAt                                                                               >  1:103922/1‑149 (MQ=255)
                                                                          tttCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAg                                                                                                                                              <  3:340325/81‑1 (MQ=255)
                                                                          tttCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAg                                                                                                                                              >  4:340325/1‑81 (MQ=255)
                                                                          tttCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGctct                                                                          >  7:222626/1‑149 (MQ=255)
                                                                            tCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAg                                                                                           <  6:79685/130‑1 (MQ=255)
                                                                            tCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAg                                                                                           >  5:79685/1‑130 (MQ=255)
                                                                               gATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCCCACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTtct                                                                     >  8:333423/1‑149 (MQ=255)
                                                                                aTCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTtctc                                                                    <  3:77379/149‑1 (MQ=255)
                                                                                  ccTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTcc                                                                   >  4:195926/1‑148 (MQ=255)
                                                                                    tGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTcc                                                                   >  1:29371/1‑146 (MQ=255)
                                                                                      ctctACATAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTaga                                                               <  5:147225/148‑1 (MQ=255)
                                                                                      ctctACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCACATGCTCTTTTCTCCTag                                                                <  4:415322/147‑1 (MQ=255)
                                                                                      ctctACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCACATGCTCTTTTCTCCTag                                                                <  4:408785/147‑1 (MQ=255)
                                                                                           cacaGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTagaaga                                                            <  2:339823/146‑1 (MQ=255)
                                                                                           cacaGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAAc                                                         <  8:222626/149‑1 (MQ=255)
                                                                                           cacaGCGTACCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTagaaga                                                            <  6:284799/146‑1 (MQ=255)
                                                                                            acaGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACt                                                        >  3:156586/1‑149 (MQ=255)
                                                                                                  tGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGcc                                                  <  2:256183/149‑1 (MQ=255)
                                                                                                    ccAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGa                                                >  5:258125/1‑149 (MQ=255)
                                                                                                    ccAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATACACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGa                                                >  1:311531/1‑149 (MQ=255)
                                                                                                       ggTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGt                                             >  7:273229/1‑149 (MQ=255)
                                                                                                        gTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTg                                            <  8:24349/149‑1 (MQ=255)
                                                                                                               gctgctTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGc                                     <  6:209673/149‑1 (MQ=255)
                                                                                                                   ctTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTtct                                 >  1:25692/1‑149 (MQ=255)
                                                                                                                   ctTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTtc                                  <  5:191990/148‑1 (MQ=255)
                                                                                                                      cAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTcc                               <  2:103922/148‑1 (MQ=255)
                                                                                                                              ttGCCTTCATCCAGCCAGGCCTTGCCTAGGCCAAACCTCAGCATCGTGTTATCCACACCTTTGGCCATCTTCTTCTCGAGCGATTCGCGCTTGGTCTTTTCTCTTAGATTATACTCAAGCCGACGTGGAACAGCTTCTCGACATCGCGa                      <  2:311531/149‑1 (MQ=255)
                                                                                                                              ttGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCacac                                                                                                                 <  2:427507/58‑1 (MQ=255)
                                                                                                                              ttGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCacac                                                                                                                 >  1:427507/1‑58 (MQ=255)
                                                                                                                                gCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGaa                     >  2:344029/1‑148 (MQ=255)
                                                                                                                                  cTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAatat                  >  8:190196/1‑149 (MQ=255)
                                                                                                                                  cTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAatat                  >  4:379021/1‑149 (MQ=255)
                                                                                                                                  cTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAatat                  >  4:381788/1‑149 (MQ=255)
                                                                                                                                  cTTCATCCAGCCAGGCCTGGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAatat                  <  2:280482/149‑1 (MQ=255)
                                                                                                                                       tCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCttt             <  3:195926/149‑1 (MQ=255)
                                                                                                                                       tCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCttt             >  7:62742/1‑149 (MQ=255)
                                                                                                                                         cagccagGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCttttt           <  3:379021/149‑1 (MQ=255)
                                                                                                                                          ggtttggCCTGGGCAAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTtctc                                >  5:32424/6‑127 (MQ=255)
                                                                                                                                          ggtttggCCTGGGCAAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTtctc                                <  6:32424/122‑1 (MQ=255)
                                                                                                                                          agccagGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTa          >  1:86963/1‑149 (MQ=255)
                                                                                                                                            tttggCCTGGGCAAGGCCAAACCTCAGCAGCGGGTTATCCCCCCCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATc        >  1:254980/4‑149 (MQ=255)
                                                                                                                                            tttggCCTGGGCAAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATc        >  5:412277/4‑149 (MQ=255)
                                                                                                                                             ttggCCTGGGCAAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATc        >  7:323231/3‑148 (MQ=255)
                                                                                                                                             ttggCCTGGGCAAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATc        >  4:186555/3‑148 (MQ=255)
                                                                                                                                               ggCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATCcac     <  5:266191/149‑1 (MQ=255)
                                                                                                                                                  cTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGGTTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATCCACAAc  >  6:108199/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CACGCTTTGCGCGCCGCCGCCAGCTCGCCCTGCAGCTGATAGGCCTTACCCAGCAGCTTCCATGCTGCCGAGTATTTCGGATCCTGCTCTACACAGCGTGCCAGGTGCACCGCTGCTTCAGCGCCGTTGCCTTCATCCAGCCAGGCCTTGCCCAGGCCAAACCTCAGCAGCGGGTTATCCACACCTTTGGCCAGCATCTTCTCCAGCGATTCGCGCATGCTCTTTTCTCCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATCCACAAC  >  NC_002947/74607‑74901

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: