Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,547,946 C→A 7.5% L71M (CTG→ATG) ‡ PP_1358 → heme iron utilization protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,547,9460CA 7.5% 292.8 / 9.9 93L71M (CTG→ATG) ‡PP_1358heme iron utilization protein
Reads supporting (aligned to +/- strand):  ref base C (46/40);  new base A (3/4);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

Only 100 of 104 total aligned reads displayed.

TGAGTACCAACGCCATTCGCCCCGCCCGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAGCGCACAAGCTGGTCGAAGAGGC  >  NC_002947/1547800‑1548091
                                                                                                                                                  |                                                                                                                                                 
tGAGTACCAACGCCATTCGCCCCGCCCGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGctgatcctga                                                                                                                                                >  2:2149751/1‑150 (MQ=255)
      ccAACGCCATTCGCCCCGCCCGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGctgatcctga                                                                                                                                                <  2:238927/144‑1 (MQ=255)
      ccAACGCCATTCGCCCCGCCCGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGctgatcctga                                                                                                                                                >  1:238927/1‑144 (MQ=255)
      ccAACGCCATTCGCCCCGCCCGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGcc                                                                                                                                          >  2:1550174/1‑150 (MQ=255)
       cAACGCCATTCGCCCCGCCCGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCg                                                                                                                                         >  1:406447/1‑150 (MQ=255)
        aaCGCCATTCGCCCCGCCCGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGc                                                                                                                                        >  1:778081/1‑150 (MQ=255)
        aaCGCCATTCGCCCCGCCCGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCAGGATCAGcggg                                                                                                                                        <  2:1705031/150‑4 (MQ=255)
              aTTCGCCCCGCCCGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGcc                                                                                                                                  <  2:2264438/150‑1 (MQ=255)
              aTTCGCCCCGCCCGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGcc                                                                                                                                  >  2:135624/1‑150 (MQ=255)
              aTTCGCCCCGCCCGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCAGGATCAGcgg                                                                                                                                         <  2:2407176/143‑3 (MQ=255)
              aTTCGCCCCGCCCGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCAGGATCAGcgg                                                                                                                                         >  1:2407176/1‑141 (MQ=255)
               ttCGCCCCGCCCGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGccc                                                                                                                                 <  1:501632/150‑1 (MQ=255)
                  gccccgcccGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGc                                                                                                                              >  2:483712/1‑150 (MQ=255)
                  gccccgcccGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGc                                                                                                                              <  2:1581240/150‑1 (MQ=255)
                    cccgcccgGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAg                                                                                                                               <  1:844884/147‑1 (MQ=255)
                    cccgcccgGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAg                                                                                                                               >  2:844884/1‑147 (MQ=255)
                     ccgcccgGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCAGGATCAGcgg                                                                                                                                         <  2:1524624/136‑3 (MQ=255)
                     ccgcccgGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCAGGATCAGcgg                                                                                                                                         >  1:1524624/1‑134 (MQ=255)
                        cccgGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACAccc                                                                                                                        >  2:2056015/1‑150 (MQ=255)
                         ccgGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGCGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCca                                                                                                                       <  2:176025/150‑1 (MQ=255)
                          cgGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGcc                                                                                                                                          <  2:1857505/130‑1 (MQ=255)
                          cgGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGcc                                                                                                                                          >  1:1857505/1‑130 (MQ=255)
                           gggAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCTTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCAGTGCCCAGCACCCCaaca                                                                                                                     >  2:1706909/1‑150 (MQ=255)
                               aCTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCt                                                                                                                 <  1:2402183/150‑1 (MQ=255)
                                 tgctgcTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCGGGTCGATGCCTGGCTACCCCTTCGGTTCGGCCATGCTGGACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGc                                                                                                               >  1:161198/1‑150 (MQ=255)
                                 tgctgcTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGc                                                                                                               >  1:2076169/1‑150 (MQ=255)
                                   ctgctCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACACACAACCTGCAg                                                                                                             >  2:1371220/1‑150 (MQ=255)
                                       tCAAGGAATACCGCGGCGTGCTTGCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAgg                                                                                                         >  2:2395327/1‑150 (MQ=255)
                                         aaGGAATACCGCGGCGTGCTTTCGACGCATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGAc                                                                                                       >  1:1037314/1‑150 (MQ=255)
                                          aGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGAcc                                                                                                      >  2:1916462/1‑150 (MQ=255)
                                          aGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGAcc                                                                                                      >  1:367589/1‑150 (MQ=255)
                                           ggAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGAccc                                                                                                     >  1:1996870/1‑150 (MQ=255)
                                                aCCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTATTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAgt                                                                                                <  1:1369092/150‑1 (MQ=255)
                                                aCCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATCGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAgt                                                                                                <  1:1947260/150‑1 (MQ=255)
                                                            tttCGACACATTCCAGGTCGATGCCTGCATACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCtgg                                                                                    >  1:1921238/1‑150 (MQ=255)
                                                                gACACATTACAGGTCGATGCCTGGCTACCCCGTCGGTTCGGTCGTGCCGGACTGCCTGGACGCCCAGCGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTggg                                                                                >  2:642810/1‑150 (MQ=255)
                                                                     aTTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACAcc                                                                                                                         <  2:1805423/104‑1 (MQ=255)
                                                                     aTTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACAcc                                                                                                                         >  1:1805423/1‑104 (MQ=255)
                                                                     aTTACAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAgc                                                                           >  1:2055849/1‑150 (MQ=255)
                                                                      ttCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAgcg                                                                          >  1:496047/1‑150 (MQ=255)
                                                                            gTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGc                                                                    >  2:1489804/1‑150 (MQ=255)
                                                                            gTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGc                                                                    <  2:496047/150‑1 (MQ=255)
                                                                               gATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGa                                                                 <  1:1550174/150‑1 (MQ=255)
                                                                               aaTGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGa                                                                 <  1:2395327/149‑1 (MQ=255)
                                                                                aTGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAg                                                                >  1:1829602/1‑150 (MQ=255)
                                                                                 tGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAgg                                                               >  2:33615/1‑150 (MQ=255)
                                                                                    cTGGCTACCCCTTCGGTTCGGTCGTGCCGTAGTGTTTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGAtg                                                            >  1:1172243/1‑150 (MQ=255)
                                                                                    cTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGAtg                                                            >  1:1479864/1‑150 (MQ=255)
                                                                                       gCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCAtt                                                                                                                                     >  2:1348920/1‑74 (MQ=255)
                                                                                       gCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCAtt                                                                                                                                     <  1:1348920/74‑1 (MQ=255)
                                                                                         tACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCtggtgg                                                                                 <  2:1190387/124‑1 (MQ=255)
                                                                                         tACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCtggtgg                                                                                 >  1:1190387/1‑124 (MQ=255)
                                                                                         tACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAg                                                       <  1:2056015/150‑1 (MQ=255)
                                                                                          aCCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAgg                                                      <  2:2376121/150‑1 (MQ=255)
                                                                                           ccccTTCGGTTCGGTCGTGCCGTACTGCCTGGCAGACTAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGc                                                     <  2:1172243/150‑1 (MQ=255)
                                                                                           ccccTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGc                                                     <  2:778081/150‑1 (MQ=255)
                                                                                           ccccTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGc                                                     >  1:118601/1‑150 (MQ=255)
                                                                                            cccTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCg                                                    >  1:563839/1‑150 (MQ=255)
                                                                                               ttcggttcggTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTa                                                 <  1:1850250/150‑1 (MQ=255)
                                                                                                tcggttcggtCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAg                                                <  2:1829602/150‑1 (MQ=255)
                                                                                                 cggttcggtCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAgg                                               <  1:33615/150‑1 (MQ=255)
                                                                                                 cggttcggtCGTGCAGTACTGCCTGAACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGACGTGTGCTCCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAgg                                               <  1:2332759/150‑1 (MQ=255)
                                                                                                   gttcggtCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGgc                                             >  2:2227723/1‑150 (MQ=255)
                                                                                                     tcggtCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCACTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGgcgc                                           <  2:2197696/150‑1 (MQ=255)
                                                                                                        gtCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAg                                                                <  1:714501/126‑1 (MQ=255)
                                                                                                        gtCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAg                                                                >  2:714501/1‑126 (MQ=255)
                                                                                                        gtCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTg                                        >  2:535643/1‑150 (MQ=255)
                                                                                                             gCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGt                                   >  2:1050050/1‑150 (MQ=255)
                                                                                                              ccGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACACACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTa                                  <  1:1371220/150‑1 (MQ=255)
                                                                                                                gTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGa                                                                       >  2:2276341/1‑111 (MQ=255)
                                                                                                                gTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGa                                                                       <  1:2276341/111‑1 (MQ=255)
                                                                                                                     gCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCg                           <  1:642810/150‑1 (MQ=255)
                                                                                                                      ccTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGa                          >  2:1104695/1‑150 (MQ=255)
                                                                                                                       cTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGaa                         <  1:1104695/150‑1 (MQ=255)
                                                                                                                       cTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGaa                         <  2:563839/150‑1 (MQ=255)
                                                                                                                        tGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAgt                                                                                                >  2:1067772/1‑78 (MQ=255)
                                                                                                                        tGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAgt                                                                                                <  1:1067772/78‑1 (MQ=255)
                                                                                                                        tGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAg                        >  1:1528219/1‑150 (MQ=255)
                                                                                                                         ggACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAgc                       >  1:2093403/1‑150 (MQ=255)
                                                                                                                          gACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAgcg                      <  1:1050050/150‑1 (MQ=255)
                                                                                                                            cGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAGCGca                    <  1:1706909/150‑1 (MQ=255)
                                                                                                                               ccAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACCAACCGCAGACGGACCCCAAGTGTCGGATGCCGGTGGGGGAGCGCGAGGGTGAGGATGTGCAGGCGGGCGGGCGCCTGACGGTCATGGCCGAAGCGCCCaa                 >  1:906696/1‑150 (MQ=255)
                                                                                                                               ccAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAGCGCACaa                 >  2:585493/1‑150 (MQ=255)
                                                                                                                                     cAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAGCGCACAAGCTGGt           <  1:483712/150‑1 (MQ=255)
                                                                                                                                      aaCCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAGCGCACAAGCTGGTc          <  2:367589/150‑1 (MQ=255)
                                                                                                                                        cccGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAgt                                                                                                >  1:306882/1‑62 (MQ=255)
                                                                                                                                        cccGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAgt                                                                                                <  2:306882/62‑1 (MQ=255)
                                                                                                                                         ccGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAGCGCACAAGCTGGTCGaa       <  2:406447/150‑1 (MQ=255)
                                                                                                                                           gCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAGCGCACAAGCTGGTCGAaga     >  1:926496/1‑150 (MQ=255)
                                                                                                                                           gCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAGCGCACAAGCTGGTCGAaga     <  1:535643/150‑1 (MQ=255)
                                                                                                                                           gCTGATCAGGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAg                                                <  2:2365998/107‑1 (MQ=255)
                                                                                                                                           gCTGATCAGGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAg                                                >  1:2365998/1‑107 (MQ=255)
                                                                                                                                            ctgatcctgatcAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAGCGCACAAGCTGGTCGAagag    <  1:2227723/150‑1 (MQ=255)
                                                                                                                                             tgatcctgatcAGCCGCATTGCCCAGCACACCCACAACCTGCCGAAGGACCCCCCGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAGCGCACAAGCTGGTCGAAGAgg   <  2:2076169/150‑1 (MQ=255)
                                                                                                                                             tgatcctgatcAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAGCGCACAAGCTGGTCGAAGAgg   <  1:2149751/150‑1 (MQ=255)
                                                                                                                                              gatcctgatcAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAGCGCACAAGCTGGTCGAAGAGGc  <  2:1921238/150‑1 (MQ=255)
                                                                                                                                               atcttgagcATCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAgg                                                               <  1:1353930/88‑1 (MQ=255)
                                                                                                                                               atcctgatcAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAgg                                                               >  2:1353930/1‑88 (MQ=255)
                                                                                                                                                tcctgatcAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAGCGCACAAGCTGGTCGaa       <  2:1548431/143‑1 (MQ=255)
                                                                                                                                                tcctgatcAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAGCGCACAAGCTGGTCGaa       >  1:1548431/1‑143 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGAGTACCAACGCCATTCGCCCCGCCCGGGAACTGCTGCTCAAGGAATACCGCGGCGTGCTTTCGACACATTCCAGGTCGATGCCTGGCTACCCCTTCGGTTCGGTCGTGCCGTACTGCCTGGACGCCCAGGGCAACCCGCTGATCCTGATCAGCCGCATTGCCCAGCACACCCACAACCTGCAGAAGGACCCCAAGTGTTCGCTGCTGGTGGGGGAGCGCGAGGCTGAGGATGTGCAGGCGGTAGGGCGCCTGACGGTAATGGCCGAAGCGCACAAGCTGGTCGAAGAGGC  >  NC_002947/1547800‑1548091

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