Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,440,866 A→G 5.2% L247P (CTG→CCG) ‡ fabB ← 3‑oxoacyl‑[acyl‑carrier‑protein] synthase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,440,8660AG 5.2% 346.8 / 6.9 116L247P (CTG→CCG) ‡fabB3‑oxoacyl‑[acyl‑carrier‑protein] synthase I
Reads supporting (aligned to +/- strand):  ref base A (60/49);  new base G (3/3);  total (64/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e-01

Only 100 of 117 total aligned reads displayed.

GTTCAGGTAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCC  >  NC_000913/2440718‑2441015
                                                                                                                                                    |                                                                                                                                                     
gTTCAGGTAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCGTCAGAGGTTGCGCCGTAGCCCACGATTTCAGCATAGATGTGAGCACCACGCGCCAg                                                                                                                                                      >  1:122993/1‑150 (MQ=255)
gTTCAGGTAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAg                                                                                                                                                      >  1:1811416/1‑150 (MQ=255)
  gcAGGTAATCGATTGGGGTATAAACGCAATGCCTCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCCTGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGAAGTGAGCACCACGCGCCAgcg                                                                                                                                                    <  2:905986/149‑1 (MQ=255)
    aGGTAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGGACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCAACACGCGCTGGcgcg                                                                                                                                                  <  1:416018/150‑1 (MQ=255)
     ggTACTCGATTGGGGTATCAACGCCCTGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGGCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCgt                                                                                                                                                 <  2:1714183/150‑1 (MQ=255)
      gTAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCgtg                                                                                                                                                <  2:1094110/150‑1 (MQ=255)
         aTCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCTGGCGCgtggtg                                                                                                                                             <  2:617588/150‑4 (MQ=255)
         aTCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTc                                                                                                                                             >  1:1441746/1‑150 (MQ=255)
          tCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTcc                                                                                                                                            <  2:70185/150‑1 (MQ=255)
          gcGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCTTGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCTGGCGCgtggtgc                                                                                                                                            <  2:1640174/149‑5 (MQ=255)
           cGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCa                                                                                                                                           <  1:1490827/150‑1 (MQ=255)
            gATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAg                                                                                                                                          <  1:1846189/150‑1 (MQ=255)
            gATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGAGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCTCGCGCCAGCGCGTGTTCCAg                                                                                                                                          <  2:853072/150‑1 (MQ=255)
             aTTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGc                                                                                                                                         <  2:478963/150‑1 (MQ=255)
               tGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGctc                                                                                                                                       >  2:889571/1‑150 (MQ=255)
                ggggTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGctct                                                                                                                                      >  2:719127/1‑150 (MQ=255)
                ggggTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCACAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGAGCGAGTTCCAGctct                                                                                                                                      >  1:191692/1‑150 (MQ=255)
                 gggTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGCCGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGAGCCAGCGCGTGTTCCAGCTCtt                                                                                                                                     >  2:1686707/1‑150 (MQ=255)
                  ggTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTc                                                                                                                                    <  1:1291251/150‑1 (MQ=255)
                      tCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCACCCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTGCCAGCTCTTCAAcc                                                                                                                                >  2:1210258/1‑150 (MQ=255)
                      tCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAccc                                                                                                                                >  2:1280026/1‑150 (MQ=255)
                       cAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTATCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCa                                                                                                                               <  2:142643/150‑1 (MQ=255)
                       cAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTGAACCa                                                                                                                               >  1:1386381/1‑150 (MQ=255)
                        aaCGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCAc                                                                                                                              >  2:1059501/1‑150 (MQ=255)
                        aaCGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCAc                                                                                                                              <  2:394127/150‑1 (MQ=255)
                         aCGCCATGCATCGCCCTCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGCTGTGCGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACt                                                                                                                             >  1:722763/1‑150 (MQ=255)
                           gCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTAc                                                                                                                           <  2:1819681/150‑1 (MQ=255)
                                 cATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATAcc                                                                                                                     >  1:346/1‑150 (MQ=255)
                                 cATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATAcc                                                                                                                     <  2:18342/150‑1 (MQ=255)
                                    cGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATAccgcc                                                                                                                  <  2:1218055/150‑1 (MQ=255)
                                    cGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATAccgcc                                                                                                                  >  2:1703059/1‑150 (MQ=255)
                                     gCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCAGAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATAccgccg                                                                                                                 >  1:965011/1‑150 (MQ=255)
                                       cATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATAccgccgcc                                                                                                               >  2:684704/1‑150 (MQ=255)
                                             cATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGc                                                                                                         >  1:1036761/1‑150 (MQ=255)
                                               tGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGa                                                                                                       >  1:821678/1‑150 (MQ=255)
                                                 ccgtgTACTGGGACGTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGGTGCGCCGTAGCCAACGATTTCAGCATAGATGTTAGCACCACGCGCCAGCGCGTGTTCCATCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATa                                                                                                     <  2:324977/146‑1 (MQ=255)
                                                 ccgCGTACTGCGCCTTCGCCAGCCGGAGCAACCATGTCTGCCCCATCCGAGGTTGCGCCGTAGCCAACGAGTTCCGCATAGCTGTGCGCACCACGCGCCCGCGCGTGGTCCCGCTCTTCCCCCCCTCCCATACCGCCGCCGCCAGCGaga                                                                                                     >  1:1412323/3‑148 (MQ=255)
                                                  aGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATaa                                                                                                    <  2:587334/150‑1 (MQ=255)
                                                      tACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCATCGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGaa                                                                                                >  2:918147/1‑150 (MQ=255)
                                                        ccgcgcCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCAACACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAAc                                                                                              <  2:1741446/148‑1 (MQ=255)
                                                          gcgcCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCg                                                                                            >  1:635419/1‑150 (MQ=255)
                                                             ccTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGCTGTGAGCACCACGCGCCAGCGCGTGGTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGCTAAAGCACCCGtgg                                                                                         >  2:347830/1‑148 (MQ=255)
                                                             ccTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCa                                                                                         >  1:1866694/1‑150 (MQ=255)
                                                              cTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTCCCCAACGATTTCAGCATAGATGTGAGCACCACAAGCCAGCGCGTGTTACAGCTCTTCAACCACTACCATACCGCCTCCGCCAGCTAAAACGAATCCGTCag                                                                                        >  2:1006357/1‑149 (MQ=255)
                                                               ttCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACg                                                                                       <  1:1562601/150‑1 (MQ=255)
                                                                tCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACgg                                                                                      <  1:587249/150‑1 (MQ=255)
                                                                  gCCAGACGGAGCAACCAAGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCCCTACCATACCGCCGCCGCCAGCGATAACGAACCCGTCACGGTg                                                                                    >  1:1701880/1‑150 (MQ=255)
                                                                   ccAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGa                                                                                   >  2:1755476/1‑150 (MQ=255)
                                                                         ggAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAACCCGTCACGGTGAGcgtcg                                                                             <  2:71989/150‑1 (MQ=255)
                                                                             aaagaaTGTTTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCCTAGATGTGAGCACCACGCGCCAGCGCGTTTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAg                                                                         <  1:73209/145‑1 (MQ=255)
                                                                                ccATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTAc                                                                      <  2:1522777/150‑1 (MQ=255)
                                                                                 cATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACg                                                                     >  2:1039129/1‑150 (MQ=255)
                                                                                     tCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAg                                                                 <  2:1498750/150‑1 (MQ=255)
                                                                                      cTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAgg                                                                >  2:669161/1‑150 (MQ=255)
                                                                                       tGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCACCCACTACCATACCGCCGCCGCCCGCGAAAACGAAACCGTCACGGAGAGCGTCGTAAGTACGGGAGGc                                                               >  1:70273/1‑150 (MQ=255)
                                                                                             aTCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGCAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTc                                                         >  1:39299/1‑150 (MQ=255)
                                                                                               cAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCCCCGCCGCCAGCGATAACGAAACCGTTACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCg                                                       >  1:1862506/1‑150 (MQ=255)
                                                                                                 gagGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCggg                                                     <  1:1224786/150‑1 (MQ=255)
                                                                                                 gagGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCggg                                                     <  1:828334/150‑1 (MQ=255)
                                                                                                  agGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCCACCACTACCATACCGCCGCCGCCCGCGCTCACGAAACAGACACGGTGAGCGTCGTACGTAACGGAGGCTTTATcgaggg                                                    >  1:300272/1‑145 (MQ=255)
                                                                                                  agGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCgggg                                                    <  2:1504819/150‑1 (MQ=255)
                                                                                                  agGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCgggg                                                    <  2:1036761/150‑1 (MQ=255)
                                                                                                   ggTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGgt                                                   >  2:545396/1‑150 (MQ=255)
                                                                                                      tGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTc                                                >  2:578413/1‑150 (MQ=255)
                                                                                                       gcgcCGTAGCCAACGATTTCAGCATAGATGGGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCg                                               <  1:919961/150‑1 (MQ=255)
                                                                                                          ccGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTg                                            <  1:848203/150‑1 (MQ=255)
                                                                                                              aGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGACGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTAtt                                        <  2:1068058/150‑1 (MQ=255)
                                                                                                               gCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTAttt                                       <  2:1441746/150‑1 (MQ=255)
                                                                                                                ccAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTta                                      <  1:1273396/150‑1 (MQ=255)
                                                                                                                ccAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTta                                      >  1:603947/1‑150 (MQ=255)
                                                                                                                 cAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTtag                                     >  2:612212/1‑150 (MQ=255)
                                                                                                                 cAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTtag                                     >  1:903717/1‑150 (MQ=255)
                                                                                                                   aCGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGCGCGTCGTAAGTACGGGAGGCTATTTCCGGGGTGACGGTGTATTtagtt                                   >  1:173702/1‑149 (MQ=255)
                                                                                                                   aCGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTtagta                                   <  1:1059501/150‑1 (MQ=255)
                                                                                                                          cAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAATACGTGAGGGAGAGCGTCGTAAGTAAGTCAGGCTTTTTCCGTGGAGTCGTTGTAATTAGTAGaaaccg                            >  2:1333829/1‑145 (MQ=255)
                                                                                                                           aGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGATAACGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCTTTGTATTTAGTAGACAgcgc                           >  2:393966/1‑150 (MQ=255)
                                                                                                                            gCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCGCGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAgcgcc                          >  1:1649244/1‑149 (MQ=255)
                                                                                                                                aGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAACCCGTGACGGTGAGCGTCGTTAGTACGGGCGGCTTTTTCCGGGGTGTCGTGGTAGTTAGTAGACAGCGCACCCa                      >  1:353513/1‑150 (MQ=255)
                                                                                                                                     tgAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGATAGCGTGACGGGGAAGGTGGTAAGTTGGGGAGGCGTTTACCCTGGTGACGGTGGAATTAGAAGACAGCGCACGCATTGCg                 >  2:321475/1‑150 (MQ=255)
                                                                                                                                     tgAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAGCCGTCAGGGGGCGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTGGGTGGATTTTGTAGACAGGGCGCGCATTGCg                 >  2:916599/1‑150 (MQ=255)
                                                                                                                                       aGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTc               <  1:551682/150‑1 (MQ=255)
                                                                                                                                       aGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTc               <  2:1265145/150‑1 (MQ=255)
                                                                                                                                       aGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGGCGTTGTATTTAGTAGACAGCGCACCCATTGCGTc               <  1:1686707/150‑1 (MQ=255)
                                                                                                                                       aGCACCACGCGCCAGCGCGTGTTCCAGCGCTTCCCCCCCTCCCATACCGCCGCCGCCAGCGATACCGAAACCGGCCCGGTGAGCGGCGTCAGGACGGGCGGCTTGTTCCGGGGTGTCGGTGTAGTTAGTCGCCCGCGCACCCCTTGCggc               >  1:148789/1‑148 (MQ=255)
                                                                                                                                       aGCACCACGCGCCAGCGCGTGTTCAAGCTCTTCAACCACTACCATACCGCCGCCGCCCGCGATAACGAAACCGTAACGGTGAGCGTAGTAAGTACGGGAGGCTTTTTCCTGGGTGTCGTTGTATTTAGTTGCAAGCGCACCCATTGCGta               >  1:143124/1‑149 (MQ=255)
                                                                                                                                        ggaacacGCGCTGGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCg              >  2:1822774/5‑150 (MQ=255)
                                                                                                                                        gctccacGCGCCAGCGCGGGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCg              <  2:1811416/147‑1 (MQ=255)
                                                                                                                                        gCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCg              >  1:955949/1‑150 (MQ=255)
                                                                                                                                        gCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCg              <  2:1603892/150‑1 (MQ=255)
                                                                                                                                         gaacacGCGCTGGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATGTAGTAGACAGCGCACCCATTGCGTAGa             >  1:651961/4‑150 (MQ=255)
                                                                                                                                          accacGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGaa            >  1:789496/1‑150 (MQ=255)
                                                                                                                                           ccacGCGCCAGCGCGTGTTCCAGCTCTTCCACCACTACCATACCGCCGCCGCCAGCGAAAACGAAACCGTCACGGAGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTGGTAATTAGTAGACAGCGCTCCCATTGCGTAGatg           >  1:923190/1‑148 (MQ=255)
                                                                                                                                           ccacGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAAt           <  1:1755476/150‑1 (MQ=255)
                                                                                                                                           ccacGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGAAAACGAAACCGTCAAGGTGAGCGTCGTAAGTACGGGAGGCTTTTTACGGGGGGTCGTTGTATTTAGTAGACAGCGCCCACATTGCGTCGAAt           >  1:1634630/1‑150 (MQ=255)
                                                                                                                                             gcgcgCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCCGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTc         >  1:381785/2‑150 (MQ=255)
                                                                                                                                              cgcgCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGATGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCg        <  1:918147/150‑1 (MQ=255)
                                                                                                                                              cgcgCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCg        <  2:1447226/150‑1 (MQ=255)
                                                                                                                                               gcgcCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGc       <  2:346/150‑1 (MQ=255)
                                                                                                                                                   cAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGc   <  2:1866694/150‑1 (MQ=255)
                                                                                                                                                    aGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGcc  >  2:624092/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTCAGGTAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCC  >  NC_000913/2440718‑2441015

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