Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,592,873 T→C T189A (ACT→GCT)  livF ← branched‑chain amino acid ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,592,8730TC98.4% 226.9 / ‑6.5 62T189A (ACT→GCT) livFbranched‑chain amino acid ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (25/36);  total (25/37)
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

ACCCGCCTAAATACGCACTTCTCACCGCTTCATTCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGTCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAGCGCACGACCAATCGCCAGCATCTGCTGTTCACCGC  >  NC_000913/3592726‑3593019
                                                                                                                                                   |                                                                                                                                                  
aCCCGCCTAAATACGCACTTCTCACCGCTTCATTCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCAt                                                                                                                                                 >  1:620218/1‑151 (MQ=255)
aCCCGCCTAAATACGCACTTCTCACCGCTTCATTCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCAt                                                                                                                                                 >  1:39232/1‑151 (MQ=255)
 cccGCCTAAATACGCACTTCTCACCGCTTCATTCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATc                                                                                                                                                >  2:105328/1‑151 (MQ=255)
 cccGCCTAAATACGCACTTCTCACCGCTTCATTCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATc                                                                                                                                                <  1:497967/151‑1 (MQ=255)
  ccGCCTAAATACGCACTTCTCACCGCTTCATTCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATcc                                                                                                                                               >  2:510795/1‑151 (MQ=255)
  ccGCCTAAATACGCACTTCTCACCGCTTCATTCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATcc                                                                                                                                               <  1:696734/151‑1 (MQ=255)
    gCCTAAATACGCACTTCTCACCGCTTCATTCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATcccc                                                                                                                                             >  1:347390/1‑151 (MQ=255)
         aaTACGCACTTCTCACCGCTTCATTCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTc                                                                                                                                        <  2:934334/151‑1 (MQ=255)
            aCGCACTTCTCACCGCTTCATTCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTcgcg                                                                                                                                     >  2:657543/1‑151 (MQ=255)
             cGCACTTCTCACCGCTTCATTCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCgcgc                                                                                                                                    >  1:1053685/1‑151 (MQ=255)
                   tctcACCGCTTCATTCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAgctgc                                                                                                                              >  1:943150/1‑151 (MQ=255)
                      cACCGCTTCATTCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAgctgc                                                                                                                              <  1:312975/148‑1 (MQ=255)
                      cACCGCTTCATTCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCg                                                                                                                           >  1:782882/1‑151 (MQ=255)
                                 tCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGa                                                                                                                 <  1:352272/150‑1 (MQ=255)
                                 tCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGa                                                                                                                 <  2:910290/150‑1 (MQ=255)
                                  cGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGaaa                                                                                                               <  2:108944/151‑1 (MQ=255)
                                   gCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGaaaa                                                                                                              <  2:47729/151‑1 (MQ=255)
                                   gCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGaaa                                                                                                               >  1:214980/1‑150 (MQ=255)
                                         agcGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTTCTCGATGGTGTCGAAAATTTGCt                                                                                                        <  2:295734/151‑1 (MQ=255)
                                         agcGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCt                                                                                                        >  2:230493/1‑151 (MQ=255)
                                           cgcATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGTCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTgg                                                                                                      <  1:755285/151‑1 (MQ=255)
                                           cgcATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTgg                                                                                                      <  2:667436/151‑1 (MQ=255)
                                             cATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGgat                                                                                                    <  1:709747/151‑1 (MQ=255)
                                                  ccAGTATCGGAAGGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATaa                                                                                               <  2:45132/151‑1 (MQ=255)
                                                     gTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCg                                                                                            <  2:1026783/151‑1 (MQ=255)
                                                       aTCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCgg                                                                                           <  1:463400/150‑1 (MQ=255)
                                                       aTCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGgc                                                                                          <  2:554479/151‑1 (MQ=255)
                                                         cGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGgcgc                                                                                        <  1:440832/151‑1 (MQ=255)
                                                           gAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCaa                                                                                      <  2:719566/151‑1 (MQ=255)
                                                             aaGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGa                                                                                    <  1:166218/151‑1 (MQ=255)
                                                                    aCATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGc                                                                              <  2:234574/150‑1 (MQ=255)
                                                                    aCATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGc                                                                              <  1:80557/150‑1 (MQ=255)
                                                                     cATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGa                                                                            <  1:510809/151‑1 (MQ=255)
                                                                        ggCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATg                                                                          <  1:288215/150‑1 (MQ=255)
                                                                         gCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGc                                                                        <  2:494694/151‑1 (MQ=255)
                                                                                 ccAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCtcatca                                                                  >  1:365205/1‑149 (MQ=255)
                                                                                  cAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGATCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGc                                                               >  1:36193/1‑151 (MQ=255)
                                                                                   aGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCa                                                              <  2:620201/151‑1 (MQ=255)
                                                                                           gCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCaa                                                       >  1:474949/1‑150 (MQ=255)
                                                                                             cgcgATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAAcg                                                    >  2:589696/1‑151 (MQ=255)
                                                                                             cgcgATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAAcg                                                    <  1:1054871/151‑1 (MQ=255)
                                                                                              gcgATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAAcgc                                                   <  1:303443/151‑1 (MQ=255)
                                                                                                    cGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTg                                             <  2:819594/151‑1 (MQ=255)
                                                                                                      ccAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCt                                           >  1:345008/1‑151 (MQ=255)
                                                                                                      ccAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCt                                           <  2:74351/151‑1 (MQ=255)
                                                                                                      ccAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCt                                           >  1:776458/1‑151 (MQ=255)
                                                                                                       cAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCtc                                          >  1:743543/1‑151 (MQ=255)
                                                                                                       cAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGAGAGCCAGCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCtc                                          <  2:411329/151‑1 (MQ=255)
                                                                                                             tAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAgc                                    >  1:526643/1‑151 (MQ=255)
                                                                                                                    tGGTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAGCGCACGAc                             >  2:149322/1‑151 (MQ=255)
                                                                                                                     ggTTGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAGCGCACGAcc                            >  2:264301/1‑151 (MQ=255)
                                                                                                                        tGGCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAGCGCACGACCAAt                         >  2:459336/1‑151 (MQ=255)
                                                                                                                          gCGTTCTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAGCGCACGACCAATCg                       <  2:145814/151‑1 (MQ=255)
                                                                                                                               cTGCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAGCGCACGACCAATCGCCAGc                  >  2:737349/1‑151 (MQ=255)
                                                                                                                                 gCTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAGCGCACGACCAATCGCCAGc                  <  2:39232/149‑1 (MQ=255)
                                                                                                                                  cTCGACGAGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAGCGCACGACCAATCGCCAGCATc               >  1:238273/1‑151 (MQ=255)
                                                                                                                                    cgacgaGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAGCGCACGACCAATCGCCAGCATct              <  1:821459/150‑1 (MQ=255)
                                                                                                                                    cgacgaGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAGCGCACGACCAATCGCCAGCATct              <  2:482585/150‑1 (MQ=255)
                                                                                                                                     gacgaGAAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAGCGCACGACCAATCGCCAGCATctgc            <  2:1062624/151‑1 (MQ=255)
                                                                                                                                        gagaAAGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAGCGCACGACCAATCGCCAGCATCTGCTGt         <  1:864621/151‑1 (MQ=255)
                                                                                                                                            aaGATAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAGCGCACGACCAATCGCCAGCATCTGCTGTTCAc     >  1:312775/1‑151 (MQ=255)
                                                                                                                                               aTAGCCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAGCGCACGACCAATCGCCAGCATCTGCTGTTCAccgc  <  2:526341/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACCCGCCTAAATACGCACTTCTCACCGCTTCATTCGCCAGCAGCGCATCACCAGTATCGGAAAGCACTACATGGCCGTTTTCCAGCACGTAGCCGCGATCCGCCAGCTTTAGCGCCTGGTTGGCGTTCTGCTCGACGAGAAAGATAGTCATCCCCTGCTCGCGCAGCTGCTCGATGGTGTCGAAAATTTGCTGGATGATAATCGGCGCAAGACCGAGCGATGGCTCATCAAGCAGTAGCAAACGCGGGTTGCTCATCAGCGCACGACCAATCGCCAGCATCTGCTGTTCACCGC  >  NC_000913/3592726‑3593019

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