Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,731,889 A→G 5.8% R465R (CGT→CGC PP_3299 ← lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,731,8890AG 5.8% 265.4 / 4.3 88R465R (CGT→CGCPP_3299lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (45/37);  new base G (3/2);  total (48/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

ACTTGGATTACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCG  >  NC_002947/3731742‑3732036
                                                                                                                                                   |                                                                                                                                                   
aCTTGGATTACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCAcg                                                                                                                                                   >  2:1205183/1‑150 (MQ=255)
  ttGGATTACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCAcgcg                                                                                                                                                 >  2:1820619/1‑150 (MQ=255)
     gATTACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCgctg                                                                                                                                              >  2:121070/1‑150 (MQ=255)
       ttACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCgctgct                                                                                                                                            >  1:2315870/1‑150 (MQ=255)
          cgcgGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCGCGTgctgctatt                                                                                                                                         <  1:2336651/150‑4 (MQ=255)
           gcgGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAg                                                                                                                                        >  2:2229360/1‑150 (MQ=255)
               tATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTg                                                                                                                                    <  1:121070/150‑1 (MQ=255)
                aTGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTAACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTgc                                                                                                                                   <  1:350375/150‑1 (MQ=255)
                aTGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGCCCTCTTTGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTgc                                                                                                                                   <  1:149975/150‑1 (MQ=255)
                     cccAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCTCAAc                                                                                                                              >  1:2210591/1‑150 (MQ=255)
                     cccAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAAc                                                                                                                              <  2:2195261/150‑1 (MQ=255)
                        aaTCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGt                                                                                                                                       <  2:2235057/138‑1 (MQ=255)
                        aaTCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGt                                                                                                                                       >  1:2235057/1‑138 (MQ=255)
                           ccTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCg                                                                                                                        >  1:1053683/1‑150 (MQ=255)
                           ccTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCg                                                                                                                        >  1:373966/1‑150 (MQ=255)
                             tCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCa                                                                                                                                         >  2:971615/1‑131 (MQ=255)
                             tCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCa                                                                                                                                         <  1:971615/131‑1 (MQ=255)
                                  gAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCTCAACTGCTCGGAAATcc                                                                                                                 <  2:2210591/150‑1 (MQ=255)
                                    aGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGa                                                                                                               <  2:181496/150‑1 (MQ=255)
                                           tACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTg                                                                                                        >  1:2045477/1‑150 (MQ=255)
                                           tACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAAATACAAGACGGGGTGCTTTAGTTTTGGCAACTTCCTGGAAATCCGGTGTCCtt                                                                                                        >  1:1282388/1‑149 (MQ=255)
                                                cttctGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAg                                                                                                   <  1:659187/150‑1 (MQ=255)
                                                    tGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAgctgc                                                                                               <  2:2315870/150‑1 (MQ=255)
                                                     gCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTc                                                                                                                         >  2:1648672/1‑123 (MQ=255)
                                                     gCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTc                                                                                                                         <  1:1648672/123‑1 (MQ=255)
                                                         cAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATc                                                                                          <  2:928455/150‑1 (MQ=255)
                                                               ttttACGCAGCCCCATCATCGCTGGGCCAAACCTGAGTCAGTTGCCCACCGGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGtt                                                                                    <  2:1282388/149‑1 (MQ=255)
                                                               aTTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGtt                                                                                    <  2:1053683/150‑1 (MQ=255)
                                                               aTTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGtt                                                                                    >  2:1518467/1‑150 (MQ=255)
                                                                   aCGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGa                                                                                <  1:1820619/150‑1 (MQ=255)
                                                                    cGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGaa                                                                               <  1:1205183/150‑1 (MQ=255)
                                                                     gCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCGCGTgctgctat                                                                                                                                          >  1:82990/1‑88 (MQ=255)
                                                                     gCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCGCGTgctgctat                                                                                                                                          <  2:82990/90‑3 (MQ=255)
                                                                       aGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGt                                                                            <  2:1362889/150‑1 (MQ=255)
                                                                           ccATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCt                                                                                                         >  2:388630/1‑117 (MQ=255)
                                                                           ccATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCt                                                                                                         <  1:388630/117‑1 (MQ=255)
                                                                           ccATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTAt                                                                        >  1:1867179/1‑150 (MQ=255)
                                                                            catcatCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATc                                                                       <  1:2229360/150‑1 (MQ=255)
                                                                             atcatcGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCa                                                                      >  2:2037890/1‑150 (MQ=255)
                                                                             atcatcGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCa                                                                      >  2:2037960/1‑150 (MQ=255)
                                                                                    cTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACCGGc                                                               <  2:2253220/150‑1 (MQ=255)
                                                                                     tgtgCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACt                                                                  >  2:1988051/1‑146 (MQ=255)
                                                                                     tgtgCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACt                                                                  <  1:1988051/146‑1 (MQ=255)
                                                                                      gtgCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTga                                                                                                       <  1:1307504/108‑1 (MQ=255)
                                                                                      gtgCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTga                                                                                                       >  2:1307504/1‑108 (MQ=255)
                                                                                         ccACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCa                                                                                                                                         >  2:1791497/1‑71 (MQ=255)
                                                                                         ccACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCa                                                                                                                                         <  1:1791497/71‑1 (MQ=255)
                                                                                         ccACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGc                                                               <  1:2195976/145‑1 (MQ=255)
                                                                                         ccACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGc                                                               >  2:2195976/1‑145 (MQ=255)
                                                                                         ccACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTacac                                                          >  2:2043368/1‑150 (MQ=255)
                                                                                           acaTCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCg                                                        >  2:630122/1‑150 (MQ=255)
                                                                                           acaTCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCg                                                        >  2:2271275/1‑150 (MQ=255)
                                                                                             aTCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGtt                                                      >  1:85408/1‑150 (MQ=255)
                                                                                             aTCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTACTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACAGCGtt                                                      >  1:534538/1‑150 (MQ=255)
                                                                                              tCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGCTCTGCGTTCAGAACGTCTATCAAACTTGCTACACCGTTc                                                     <  2:799361/150‑1 (MQ=255)
                                                                                              tCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTc                                                     >  1:816679/1‑150 (MQ=255)
                                                                                               cTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCg                                                                                                                        >  1:1850848/1‑82 (MQ=255)
                                                                                               cTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCg                                                                                                                        <  2:1850848/82‑1 (MQ=255)
                                                                                                tGAGTCAGTTGACCACCCGCCGCCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTg                                                   <  2:2045477/150‑1 (MQ=255)
                                                                                                tGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTg                                                   <  2:1867179/150‑1 (MQ=255)
                                                                                                   gTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGaa                                                <  1:2271275/150‑1 (MQ=255)
                                                                                                    tCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTc                                                                           >  2:2151942/1‑122 (MQ=255)
                                                                                                    tCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTc                                                                           <  1:2151942/122‑1 (MQ=255)
                                                                                                     cAGTTGACTACCCTCCACCCAGGGCGCGGTAGGCAGCACTAGCGGTACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGCTGTCCTGAGACGCTGCTGATTTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATg                                              <  1:1121964/150‑1 (MQ=255)
                                                                                                       gTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGtt                                                                                    <  1:1420671/110‑1 (MQ=255)
                                                                                                       gTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGtt                                                                                    >  2:1420671/1‑110 (MQ=255)
                                                                                                          gACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCCGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACg                                         >  1:2302202/1‑150 (MQ=255)
                                                                                                            ccaccCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGaa                                       >  1:1551682/1‑150 (MQ=255)
                                                                                                             caccCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGc                                                                                       <  2:117430/101‑1 (MQ=255)
                                                                                                             caccCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGc                                                                                       >  1:117430/1‑101 (MQ=255)
                                                                                                             caccCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAg                                      >  2:718861/1‑150 (MQ=255)
                                                                                                             caccCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAg                                      >  1:1472090/1‑150 (MQ=255)
                                                                                                              acacGCCACCCAGGGCGAGGTAGGCAGCAATAGAAGAACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCTATGTCTTGAACAGCTGGTGATCTGCGTTCAGAACGTCGATCAAACGGGCGACACCGTTCTGGAATTCGACGAAgg                                     >  2:190836/4‑150 (MQ=255)
                                                                                                                ccGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGAc                                   >  1:1917944/1‑150 (MQ=255)
                                                                                                                ccGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGAc                                   >  2:2258787/1‑150 (MQ=255)
                                                                                                                 cGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACg                                  >  1:828791/1‑150 (MQ=255)
                                                                                                                        cAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAAcgtcg                           <  2:828791/150‑1 (MQ=255)
                                                                                                                          gggCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTg                         >  1:879480/1‑150 (MQ=255)
                                                                                                                               aGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTgctcgc                    >  1:1462907/1‑150 (MQ=255)
                                                                                                                                  tAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTcc                 >  2:917895/1‑150 (MQ=255)
                                                                                                                                     gcggcaATAGCGGCTCGCGCTGCTTCAGTTTGTGCAACTGCTCGGAAATCCGATGTTCTGAGAAGCTGCTCATCTTCGTTCAGTACGTCTATAAACCTGGCTCCACCGTTCTGTAATGCTACGAAGGACGAATGTCGTTCTCGCTCCAgg              <  2:624189/147‑1 (MQ=255)
                                                                                                                                     gcagcaATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAgg              <  1:718861/150‑1 (MQ=255)
                                                                                                                                          tggagcagcGCGTGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATAc         >  1:2404444/4‑150 (MQ=255)
                                                                                                                                          aaTAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAAcgtc                            >  1:499509/1‑131 (MQ=255)
                                                                                                                                          aaTAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAAcgtc                            <  2:499509/131‑1 (MQ=255)
                                                                                                                                           gaagcagcGCGTGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATAcc        >  2:2350518/3‑150 (MQ=255)
                                                                                                                                            tAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTg                                                   <  2:609920/106‑1 (MQ=255)
                                                                                                                                            tAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTg                                                   >  1:609920/1‑106 (MQ=255)
                                                                                                                                                 gcaCGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACt                                                                  <  2:69779/86‑1 (MQ=255)
                                                                                                                                                 gcaCGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACt                                                                  >  1:69779/1‑86 (MQ=255)
                                                                                                                                                 gcaCGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACcg  >  1:55918/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACTTGGATTACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCG  >  NC_002947/3731742‑3732036

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