Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,748,112 C→A 9.9% V289V (GTG→GTT etfA ← electron transfer flavoprotein subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,748,1120CA 9.9% 595.3 / 40.8 189V289V (GTG→GTTetfAelectron transfer flavoprotein subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (86/84);  new base A (13/6);  total (99/90)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

GTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCAT  >  NC_002947/4747969‑4748261
                                                                                                                                               |                                                                                                                                                     
 tCATCCCTCACCCCTCAAGCCATAAAAATACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGa                                                                                                                                                <  2:4968428/150‑1 (MQ=255)
 tCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGa                                                                                                                                                <  1:3680179/150‑1 (MQ=255)
       cTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCg                                                                                                                                          >  1:4624335/1‑150 (MQ=255)
        tCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAACCGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCgg                                                                                                                                         <  1:4245500/150‑1 (MQ=255)
        tCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCgg                                                                                                                                         >  2:3876643/1‑150 (MQ=255)
              cTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGcttc                                                                                                                                   <  2:5059594/150‑1 (MQ=255)
               tCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGcttct                                                                                                                                  <  1:4103856/150‑1 (MQ=255)
                cAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGcttctt                                                                                                                                 >  2:2085857/1‑150 (MQ=255)
                 aaGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCTGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCttcttc                                                                                                                                <  2:1010685/150‑1 (MQ=255)
                 aaGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCttcttc                                                                                                                                <  2:4463919/150‑1 (MQ=255)
                 aaGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCttcttc                                                                                                                                <  1:2545141/150‑1 (MQ=255)
                      aTAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCt                                                                                                                           >  2:4353699/1‑150 (MQ=255)
                       tAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCtt                                                                                                                          >  1:3859034/1‑150 (MQ=255)
                       tAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCtt                                                                                                                          >  2:3758870/1‑150 (MQ=255)
                        aaaaaaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCttg                                                                                                                         >  1:514679/1‑150 (MQ=255)
                         aaaaaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCttgt                                                                                                                        <  1:411760/150‑1 (MQ=255)
                         aaaaaaCCCGACCCTCAAAACGAGGGCCGGGTTCTGTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGAACGGCGCTACATCGTCCtagt                                                                                                                        >  1:2303009/1‑150 (MQ=255)
                              aCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATc                                                                                                                   <  1:3303541/150‑1 (MQ=255)
                                ccGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGc                                                                                                                 <  1:470555/150‑1 (MQ=255)
                                 cGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGcc                                                                                                                <  2:1886971/150‑1 (MQ=255)
                                  gACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCa                                                                                                               >  2:4879082/1‑150 (MQ=255)
                                     ccTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGa                                                                                                            >  2:898659/1‑150 (MQ=255)
                                     ccTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGa                                                                                                            <  2:1729075/150‑1 (MQ=255)
                                     ccTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGa                                                                                                            >  2:2114820/1‑150 (MQ=255)
                                       tCAAACCGCGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCATAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATc                                                                                                          <  2:1791641/150‑1 (MQ=255)
                                       tCAAAACGAGGGCCGGGTTCTTTATAAGTGGAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATc                                                                                                          >  1:3145/1‑150 (MQ=255)
                                       tCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCttcttc                                                                                                                                <  2:522719/128‑1 (MQ=255)
                                       tCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATc                                                                                                          >  1:3440750/1‑150 (MQ=255)
                                       tCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATc                                                                                                          >  1:1779259/1‑150 (MQ=255)
                                        cAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCa                                                                                                         >  2:2339980/1‑150 (MQ=255)
                                            aCGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACttt                                                                                                     <  1:1307150/150‑1 (MQ=255)
                                              gAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTgg                                                                                                   <  1:3943617/150‑1 (MQ=255)
                                              gAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTgg                                                                                                   <  1:1725908/150‑1 (MQ=255)
                                               aGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                                                                  <  1:1606763/150‑1 (MQ=255)
                                                  gCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTc                                                                                               <  2:3153827/150‑1 (MQ=255)
                                                   ccGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCt                                                                                              >  2:832267/1‑150 (MQ=255)
                                                     gggTTAATTATAAGTGAAGCAGGCGGATCAGACCACCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCttt                                                                                            <  2:222723/150‑1 (MQ=255)
                                                      ggTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTc                                                                                           <  1:4162581/150‑1 (MQ=255)
                                                      ggTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTc                                                                                           <  2:3856829/150‑1 (MQ=255)
                                                          cTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGc                                                                                       <  1:1092738/150‑1 (MQ=255)
                                                           tttATAAGTGAAGCAGGCGGATCAGACCGGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGcc                                                                                      >  1:4386714/1‑150 (MQ=255)
                                                           tttATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGcc                                                                                      >  1:5078795/1‑150 (MQ=255)
                                                             tataAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCgg                                                                                    >  2:4485557/1‑150 (MQ=255)
                                                             tataAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCgg                                                                                    >  1:5050614/1‑150 (MQ=255)
                                                                    gAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTg                                                                             >  1:3355694/1‑150 (MQ=255)
                                                                       gCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTg                                                                          <  2:2203523/150‑1 (MQ=255)
                                                                         aGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGa                                                                        <  1:2079403/150‑1 (MQ=255)
                                                                           tataaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATc                                                                      <  1:76004/146‑1 (MQ=255)
                                                                           gCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGACGATCGGCGCTTCTTCGTCCTTGTTGATCGCCCCGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATc                                                                      >  2:2255009/1‑150 (MQ=255)
                                                                           gCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATc                                                                      <  2:4227710/150‑1 (MQ=255)
                                                                                tCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCAcc                                                                 >  2:4579444/1‑150 (MQ=255)
                                                                                 cAGACCAGCTTTTCCAGCTCCGGAACCGCTTAGACCAGGTCTGCGACCAGGCCGTAGTCGGCCACCTGGATGATCGGCGCTTATTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATcgcgaag                                                                <  2:462529/150‑4 (MQ=255)
                                                                                   gACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGga                                                              >  2:1438837/1‑150 (MQ=255)
                                                                                   gACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGga                                                              >  2:1520723/1‑150 (MQ=255)
                                                                                       aGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATa                                                          <  2:3200797/150‑1 (MQ=255)
                                                                                        gCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATAc                                                         <  2:1168408/150‑1 (MQ=255)
                                                                                              ccAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATCCCGACgg                                                   >  2:1266325/1‑150 (MQ=255)
                                                                                                aGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCg                                                 >  1:657689/1‑150 (MQ=255)
                                                                                                aGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCg                                                 <  1:3284640/150‑1 (MQ=255)
                                                                                                  cTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGAt                                               <  2:4780768/150‑1 (MQ=255)
                                                                                                   tCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATg                                              <  1:1520093/150‑1 (MQ=255)
                                                                                                        aaCCGCTTCGAACAGGTCCGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTCCAg                                         >  2:3799926/1‑150 (MQ=255)
                                                                                                           cGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCtg                                      >  2:1049245/1‑150 (MQ=255)
                                                                                                             cggtgAAGAGATCACCGACGAGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCtccg                                    <  2:1452006/146‑4 (MQ=255)
                                                                                                              ttCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGGAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTgg                                   >  1:4085517/1‑150 (MQ=255)
                                                                                                              ttCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTgg                                   >  1:339159/1‑150 (MQ=255)
                                                                                                                cGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGgcg                                 >  1:1189259/1‑150 (MQ=255)
                                                                                                                cGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGgcg                                 >  2:3379856/1‑150 (MQ=255)
                                                                                                                 gAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATAGGAGCTTCTTCTTCCTAGTTGATCGCCACTATCACTTTGGAGTCTTTCATGCAGTCCAGGTGCTGGATCTCGCCGGAGATACCGACGGCGATTAACAGCTGTGgcgc                                >  2:637514/1‑150 (MQ=255)
                                                                                                                      gatCACCGACGAGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAAcca                           <  2:93387/148‑3 (MQ=255)
                                                                                                                      gatCACCGACGAGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAAcca                           >  2:2436956/3‑148 (MQ=255)
                                                                                                                       gTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGAt                          >  2:4689874/1‑150 (MQ=255)
                                                                                                                         cAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCt                        >  1:3714954/1‑150 (MQ=255)
                                                                                                                          gccgacgaGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCtt                       <  1:4657204/148‑1 (MQ=255)
                                                                                                                          aGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCtt                       <  1:1411642/150‑1 (MQ=255)
                                                                                                                           gCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTg                      >  1:52479/1‑150 (MQ=255)
                                                                                                                           ccgacgaGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGAca                                                    <  2:2957345/119‑2 (MQ=255)
                                                                                                                             gACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGcc                    <  1:4879082/150‑1 (MQ=255)
                                                                                                                              acgaGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGccc                   >  1:2959808/4‑149 (MQ=255)
                                                                                                                                gaGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGCCCGt                 >  2:2808735/2‑150 (MQ=255)
                                                                                                                                cAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGt                 >  1:2199459/1‑150 (MQ=255)
                                                                                                                                 aGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATTCCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTc                <  1:2717409/150‑1 (MQ=255)
                                                                                                                                  ggCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCt               <  2:2741935/150‑1 (MQ=255)
                                                                                                                                  ggCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCt               >  2:4508514/1‑150 (MQ=255)
                                                                                                                                  ggCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGCCCGTCt               >  1:1124473/1‑150 (MQ=255)
                                                                                                                                   gCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATc                         >  1:4575505/1‑139 (MQ=255)
                                                                                                                                    ccataatCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGCCCGTCTgg             <  1:2867795/144‑1 (MQ=255)
                                                                                                                                    ccGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTgg             <  1:1063073/150‑1 (MQ=255)
                                                                                                                                     cataatCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGCCCGTCTGGc            >  2:3027425/6‑150 (MQ=255)
                                                                                                                                     cataatCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGCCCGTCTGGc            >  1:3412783/6‑150 (MQ=255)
                                                                                                                                     cataatCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGCCCGTCTGGc            >  1:4803423/6‑150 (MQ=255)
                                                                                                                                      gTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGcc           <  2:3242108/150‑1 (MQ=255)
                                                                                                                                      ataatCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGCCCGTCTGGcc           >  1:1149621/5‑150 (MQ=255)
                                                                                                                                         gTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGAc        >  1:3542281/1‑150 (MQ=255)
                                                                                                                                           cGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGTCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCt      <  2:2303009/150‑1 (MQ=255)
                                                                                                                                           cGGCCACCTGGAAGATCGGCGATTCTTCGTCCTTGTTGATCACCACGATCACTTTGGAGTCTTTAATTCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGAGCTACGATCTTGCCGGTCTGGCCGACCt      <  1:305380/150‑1 (MQ=255)
                                                                                                                                             gCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGCCGGTCTGGCCAACCTGc    <  1:3265242/150‑1 (MQ=255)
                                                                                                                                              ccACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGCCGGTCTGGCCAACCTGcc   >  2:3624878/1‑149 (MQ=255)
                                                                                                                                               cACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGCCGGTCTGGCCAACCTGcca  <  2:3968120/150‑3 (MQ=255)
                                                                                                                                               cACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCGGGTTGGAGTCTTTCATGCCGGCCAGGTGCTGGAGCGCGCCGGAGATACCGACGGCGAGGTACTGCTGTGGCGCAACGGTCTTGCCGGTCTGGCCGAGCTGCag  >  2:3345429/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCAT  >  NC_002947/4747969‑4748261

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