Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,550,731 A→G V189A (GTA→GCA)  gntP ← fructuronate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,550,7310AG99.2% 477.5 / NA 122V189A (GTA→GCA) gntPfructuronate transporter
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (60/61);  total (60/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 123 total aligned reads displayed.

ATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTC  >  NC_000913/4550584‑4550877
                                                                                                                                                   |                                                                                                                                                  
aTGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTgcac                                                                                                                                                 >  4:714531/1‑151 (MQ=255)
  gATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACtt                                                                                                                                               >  1:1156313/1‑151 (MQ=255)
  gATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACtt                                                                                                                                               >  1:1681391/1‑151 (MQ=255)
   aTCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGTGGAACTTCCGCATTATCACACCTGCGCACATTACACttc                                                                                                                                              <  1:1077695/151‑2 (MQ=255)
   aTCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTg                                                                                                                                              <  1:600733/151‑1 (MQ=255)
     cGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTggg                                                                                                                                            <  4:1341507/151‑1 (MQ=255)
      gCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTggg                                                                                                                                            >  4:1519869/1‑150 (MQ=255)
        cGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATc                                                                                                                                         <  1:1913680/151‑1 (MQ=255)
            aTCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTc                                                                                                                                      >  4:621504/1‑150 (MQ=255)
            aTCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCa                                                                                                                                     <  4:1829788/151‑1 (MQ=255)
               aGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCAcc                                                                                                                                   >  4:1299164/1‑150 (MQ=255)
               aGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCAcc                                                                                                                                   >  4:1721177/1‑150 (MQ=255)
                gCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAg                                                                                                                                 <  4:662036/151‑1 (MQ=255)
                 cGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGt                                                                                                                                >  4:1700362/1‑151 (MQ=255)
                   gCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTAc                                                                                                                              <  4:298989/151‑1 (MQ=255)
                     aCCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGc                                                                                                                            >  4:1349180/1‑151 (MQ=255)
                          aaTCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCAt                                                                                                                         >  1:1793947/1‑149 (MQ=255)
                          aaTCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAg                                                                                                                       <  1:1504337/151‑1 (MQ=255)
                           aTCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGa                                                                                                                      <  2:192099/151‑1 (MQ=255)
                             cGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGata                                                                                                                    <  1:1265808/151‑1 (MQ=255)
                             cGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGat                                                                                                                     <  4:1785556/150‑1 (MQ=255)
                             cGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGat                                                                                                                     <  1:1075733/150‑1 (MQ=255)
                               aaaCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTAGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGAtata                                                                                                                  <  4:633643/151‑1 (MQ=255)
                                aaCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGAtata                                                                                                                  >  1:295645/1‑150 (MQ=255)
                                  cGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAAc                                                                                                               >  1:689101/1‑151 (MQ=255)
                                     cGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCAt                                                                                                            <  1:1075455/151‑1 (MQ=255)
                                     cGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCAt                                                                                                            <  4:1185847/151‑1 (MQ=255)
                                     cGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCAt                                                                                                            <  1:1598568/151‑1 (MQ=255)
                                      gAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATc                                                                                                           >  4:1019371/1‑151 (MQ=255)
                                      gAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATc                                                                                                           <  2:222021/151‑1 (MQ=255)
                                       acagatgtcAGAATATTCATATCTACCCGTTGAGCTGCTTACAGGAATGATGGCGTCGGGCTCTTAATATTTCCGACGACCTCCTGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACTCCTTAGATATAAACCATcc                                                                                                          <  3:541953/143‑1 (MQ=255)
                                        aaGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATccc                                                                                                         >  1:1158405/1‑151 (MQ=255)
                                         agagGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATcccc                                                                                                        <  1:898902/151‑1 (MQ=255)
                                             ggCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCatat                                                                                                    <  4:1508404/151‑1 (MQ=255)
                                                aGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCcg                                                                                                 <  1:1540298/151‑1 (MQ=255)
                                                aGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCcg                                                                                                 >  1:117889/1‑151 (MQ=255)
                                                 gATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCcgc                                                                                                <  1:262053/151‑1 (MQ=255)
                                                  attattCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCcgcg                                                                                               <  1:1817366/151‑1 (MQ=255)
                                                  attattCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCcgcg                                                                                               <  1:189287/151‑1 (MQ=255)
                                                       tCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCAt                                                                                           <  1:1313602/150‑1 (MQ=255)
                                                        cATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATa                                                                                          >  2:98061/1‑150 (MQ=255)
                                                           atCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGc                                                                                       <  2:492463/150‑1 (MQ=255)
                                                           atCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGc                                                                                       >  1:1655328/1‑150 (MQ=255)
                                                            tCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCAt                                                                                     >  1:870841/1‑151 (MQ=255)
                                                             cTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCAtt                                                                                    >  1:1658929/1‑151 (MQ=255)
                                                              tACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTc                                                                                   >  1:504997/1‑151 (MQ=255)
                                                                  ggTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCa                                                                               <  2:250186/151‑1 (MQ=255)
                                                                   gTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAg                                                                              <  4:1801623/151‑1 (MQ=255)
                                                                    ttGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAgc                                                                             >  3:356692/1‑151 (MQ=255)
                                                                      gATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCCGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAgcgc                                                                           >  4:1172086/1‑151 (MQ=255)
                                                                      gATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAgcgc                                                                           >  4:1437985/1‑151 (MQ=255)
                                                                          tGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCAcc                                                                       >  2:483250/1‑151 (MQ=255)
                                                                           gCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCAcc                                                                       <  4:907017/150‑1 (MQ=255)
                                                                           gCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCAc                                                                        <  4:1053066/149‑1 (MQ=255)
                                                                           gCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCg                                                                      >  1:211200/1‑151 (MQ=255)
                                                                              ttCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGa                                                                    <  2:278359/150‑1 (MQ=255)
                                                                              ttCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGAc                                                                   >  1:1273459/1‑151 (MQ=255)
                                                                               tCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGAc                                                                   <  2:356703/150‑1 (MQ=255)
                                                                                cAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGAccc                                                                 <  1:963797/151‑1 (MQ=255)
                                                                                cAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGAccc                                                                 <  1:576737/151‑1 (MQ=255)
                                                                                 aGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCg                                                                >  4:1496574/1‑151 (MQ=255)
                                                                                        atgGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCg                                                          >  4:88562/1‑150 (MQ=255)
                                                                                          ggCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGt                                                        <  1:621243/150‑1 (MQ=255)
                                                                                            cGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTggg                                                     <  1:1462558/151‑1 (MQ=255)
                                                                                            cGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTggg                                                     >  3:519365/1‑151 (MQ=255)
                                                                                            cGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTggg                                                     >  1:909932/1‑151 (MQ=255)
                                                                                            cGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTggg                                                     >  4:1144847/1‑151 (MQ=255)
                                                                                             gTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGa                                                    <  4:490325/151‑1 (MQ=255)
                                                                                               tGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAAc                                                  >  1:57615/1‑151 (MQ=255)
                                                                                                gggCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAAc                                                  <  4:1829659/150‑1 (MQ=255)
                                                                                                gggCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAAc                                                  >  1:1750130/1‑150 (MQ=255)
                                                                                                gggCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACa                                                 <  4:1662799/151‑1 (MQ=255)
                                                                                                     cTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTg                                              >  1:1501070/1‑149 (MQ=255)
                                                                                                     cTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGa                                             >  1:1880822/1‑150 (MQ=255)
                                                                                                      tCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAAt                                           >  1:318184/1‑151 (MQ=255)
                                                                                                      tCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAAt                                           >  1:25069/1‑151 (MQ=255)
                                                                                                       cAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGaa                                            <  4:1865473/149‑1 (MQ=255)
                                                                                                             tGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCgg                                     <  1:77804/150‑1 (MQ=255)
                                                                                                             tGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGt                                    >  2:402370/1‑151 (MQ=255)
                                                                                                                cGAGGAACTTCGGCAGGGTCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAg                                  >  3:75051/1‑150 (MQ=255)
                                                                                                                 gAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTg                                <  4:1697523/151‑1 (MQ=255)
                                                                                                                 gAGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCAGTGGGAACAGTGAATGTGCGGTAGTg                                >  1:378882/1‑151 (MQ=255)
                                                                                                                  aGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTgg                               <  3:308616/151‑1 (MQ=255)
                                                                                                                  aGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTgg                               <  2:475871/151‑1 (MQ=255)
                                                                                                                  aGGAACTTCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTgg                               <  1:1437405/151‑1 (MQ=255)
                                                                                                                        ttCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAgct                          >  4:935723/1‑150 (MQ=255)
                                                                                                                         tCGGCAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAgctgc                        >  2:42846/1‑151 (MQ=255)
                                                                                                                             cAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACt                    <  4:733331/151‑1 (MQ=255)
                                                                                                                             cAGGATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACt                    >  1:577796/1‑151 (MQ=255)
                                                                                                                               ggATCAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACTGc                  >  1:216714/1‑151 (MQ=255)
                                                                                                                                   cAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCggg              >  2:116238/1‑151 (MQ=255)
                                                                                                                                   cAGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCgg               <  4:1895247/150‑1 (MQ=255)
                                                                                                                                    aGACCTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCggg              >  4:490249/1‑150 (MQ=255)
                                                                                                                                       ccTGCGCAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCg          <  1:297560/151‑1 (MQ=255)
                                                                                                                                          gcgcAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGcc        >  1:944585/1‑150 (MQ=255)
                                                                                                                                             cAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGtt    <  2:410320/151‑1 (MQ=255)
                                                                                                                                             cAGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGtt    >  4:1491608/1‑151 (MQ=255)
                                                                                                                                              aGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGttt   <  1:1519269/151‑1 (MQ=255)
                                                                                                                                              aGATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGttt   >  1:351230/1‑151 (MQ=255)
                                                                                                                                               gATTGCACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTc  <  1:689597/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTC  >  NC_000913/4550584‑4550877

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑