Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,731,892 C→T 5.7% A464A (GCG→GCA)  PP_3299 ← lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,731,8920CT 5.7% 264.7 / 4.0 89A464A (GCG→GCA) PP_3299lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (44/39);  new base T (3/2);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

TTGGATTACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTT  >  NC_002947/3731744‑3732040
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                                             <  1:971615/131‑1 (MQ=255)
                           tCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCa                                                                                                                                             >  2:971615/1‑131 (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                                                                                                                             <  1:1648672/123‑1 (MQ=255)
                                                   gCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTc                                                                                                                             >  2:1648672/1‑123 (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                                                                                                             <  1:388630/117‑1 (MQ=255)
                                                                         ccATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCt                                                                                                             >  2:388630/1‑117 (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                                                                      <  1:1988051/146‑1 (MQ=255)
                                                                                   tgtgCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACt                                                                      >  2:1988051/1‑146 (MQ=255)
                                                                                    gtgCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTga                                                                                                           >  2:1307504/1‑108 (MQ=255)
                                                                                    gtgCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTga                                                                                                           <  1:1307504/108‑1 (MQ=255)
                                                                                       ccACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCa                                                                                                                                             <  1:1791497/71‑1 (MQ=255)
                                                                                       ccACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCa                                                                                                                                             >  2:1791497/1‑71 (MQ=255)
                                                                                       ccACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGc                                                                   >  2:2195976/1‑145 (MQ=255)
                                                                                       ccACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGc                                                                   <  1:2195976/145‑1 (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                                                                                                                            <  2:1850848/82‑1 (MQ=255)
                                                                                             cTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCg                                                                                                                            >  1:1850848/1‑82 (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                                                                                        >  2:1420671/1‑110 (MQ=255)
                                                                                                     gTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGtt                                                                                        <  1:1420671/110‑1 (MQ=255)
                                                                                                        gACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCCGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACg                                             >  1:2302202/1‑150 (MQ=255)
                                                                                                          ccaccCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGaa                                           >  1:1551682/1‑150 (MQ=255)
                                                                                                           caccCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGc                                                                                           >  1:117430/1‑101 (MQ=255)
                                                                                                           caccCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGc                                                                                           <  2:117430/101‑1 (MQ=255)
                                                                                                           caccCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAg                                          >  1:1472090/1‑150 (MQ=255)
                                                                                                           caccCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAg                                          >  2:718861/1‑150 (MQ=255)
                                                                                                            acacGCCACCCAGGGCGAGGTAGGCAGCAATAGAAGAACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCTATGTCTTGAACAGCTGGTGATCTGCGTTCAGAACGTCGATCAAACGGGCGACACCGTTCTGGAATTCGACGAAgg                                         >  2:190836/4‑150 (MQ=255)
                                                                                                              ccGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGAc                                       >  2:2258787/1‑150 (MQ=255)
                                                                                                              ccGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGAc                                       >  1:1917944/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                                <  2:499509/131‑1 (MQ=255)
                                                                                                                                        aaTAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAAcgtc                                >  1:499509/1‑131 (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)
                                                                                                                                                  cgcgCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGt   <  1:630122/150‑1 (MQ=255)
                                                                                                                                                   gcgcTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGAGGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGtt  >  1:1951730/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTGGATTACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTT  >  NC_002947/3731744‑3732040

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: ‑