Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 82,161 C→A 7.0% G457C (GGT→TGT)  leuA ← 2‑isopropylmalate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091382,1610CA 7.0% 307.0 / 8.2 97G457C (GGT→TGT) leuA2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base C (41/49);  new base A (4/3);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

TCGACTTCTGCGGCACGCCAGATATTGTTCAGAACGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGCCACAGGCCAGTTTGACGGCGG  >  NC_000913/82013‑82307
                                                                                                                                                    |                                                                                                                                                  
tCGACTTCTGCGGCACGCCAGATATTGTTCAGAACGTGCACCATGGCTTTGGCAGATGACTCAACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATAACCACCTGACCCAGCGCATCTTTAcc                                                                                                                                                    >  4:440428/1‑149 (MQ=255)
 cGACTTCTGCGGCACGCCAGATATTGTTCAGAACGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCg                                                                                                                                                   >  3:427054/1‑149 (MQ=255)
  gACTTCTGCGGCACGCCAGATATTGTTCAGAACGTGCTCCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCgg                                                                                                                                                  <  7:352834/149‑2 (MQ=255)
   aCTTCTGCGGCACGCCAGATTTTGTTCAGAACGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTg                                                                                                                                                 <  4:56939/149‑1 (MQ=255)
   aCTTCTGCGGCACGCCAGATATTGTTCAGAACGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTg                                                                                                                                                 <  7:275014/149‑1 (MQ=255)
    cTTCTGCGGCACGCCAGATATTGTTCAGAACGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCCTTACCGTgg                                                                                                                                                >  8:348099/1‑149 (MQ=255)
     ttCTGCGGCACGCCAGATATTGTTCAGACCGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGGTAGCGACGATATCCACCTGACCCCGCGCATCTTTACCGTGGc                                                                                                                                               >  5:422146/1‑149 (MQ=255)
                  cAGATATTGTTCAGAACGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAg                                                                                                                                  >  2:485364/1‑149 (MQ=255)
                  cAGATATTGTGCAGATCGGGCACCATGGATTTGGCAGAGGACTCGACAATATCGGTAGCCAGGCCGACGCCGGGTAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAg                                                                                                                                  <  6:67311/149‑1 (MQ=255)
                    gATATTGTTCAGAACGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGc                                                                                                                                <  2:40369/149‑1 (MQ=255)
                     atatTGTTCAGAACGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCt                                                                                                                               <  3:108784/149‑1 (MQ=255)
                        ttGTTCAGAACGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGt                                                                                                                             <  5:157404/148‑1 (MQ=255)
                        ttGTTCAGAACGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCCTTACCGTGGCCTTTGGCGGTCAGGCTGTa                                                                                                                            <  7:348099/149‑1 (MQ=255)
                          gTTCAGAACGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTAt                                                                                                                           <  8:78361/148‑1 (MQ=255)
                              aGACCGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCAc                                                                                                                      <  8:468309/149‑1 (MQ=255)
                              aGAACGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCAc                                                                                                                      >  1:427306/1‑149 (MQ=255)
                              aGAACGTGCACCATGGCGTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCa                                                                                                                       <  6:123772/148‑1 (MQ=255)
                               gAACGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCAc                                                                                                                      >  4:443047/1‑148 (MQ=255)
                               gAACGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCAc                                                                                                                      <  3:443047/148‑1 (MQ=255)
                                 aCGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAg                                                                                                                   <  4:206176/149‑1 (MQ=255)
                                     gCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTc                                                                                                                >  2:33695/1‑148 (MQ=255)
                                        ccATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACg                                                                                                            <  6:194658/149‑1 (MQ=255)
                                             gCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTtat                                                                                                        <  3:350454/148‑1 (MQ=255)
                                                   gCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAg                                                                                                  <  1:166883/148‑1 (MQ=255)
                                                         gACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAg                                                                                                                                  >  1:57777/1‑110 (MQ=255)
                                                         gACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAg                                                                                                                                  <  2:57777/110‑1 (MQ=255)
                                                         gACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAg                                                                                                  <  8:77792/142‑1 (MQ=255)
                                                         gACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAg                                                                                                  >  7:77792/1‑142 (MQ=255)
                                                           cTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACAGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTtata                                                                                                       >  7:141950/1‑135 (MQ=255)
                                                           cTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACAGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTtata                                                                                                       <  8:141950/135‑1 (MQ=255)
                                                           cTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCATTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGt                                                                                         <  8:466183/149‑1 (MQ=255)
                                                            tCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGTGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACAGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGt                                                                                         >  7:93478/1‑148 (MQ=255)
                                                              gACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTaa                                                                                      >  8:102001/1‑149 (MQ=255)
                                                              gACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTaa                                                                                      >  1:283982/1‑149 (MQ=255)
                                                               aCAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAAt                                                                                     >  4:331507/1‑149 (MQ=255)
                                                               aCAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAAt                                                                                     >  2:434129/1‑149 (MQ=255)
                                                                cAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAAtt                                                                                    <  3:401174/149‑1 (MQ=255)
                                                                  atatCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGc                                                                                  <  1:485364/149‑1 (MQ=255)
                                                                   tatCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTg                                                                                                                                          >  8:255933/1‑92 (MQ=255)
                                                                   tatCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTg                                                                                                                                          <  7:255933/92‑1 (MQ=255)
                                                                   tatCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGcc                                                                                 <  8:101821/149‑1 (MQ=255)
                                                                   tatCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGcc                                                                                 >  1:472201/1‑149 (MQ=255)
                                                                     tCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCt                                                                                <  1:434129/148‑1 (MQ=255)
                                                                     tCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTg                                                                               <  7:102001/149‑1 (MQ=255)
                                                                       ggTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGAt                                                                             <  2:427306/149‑1 (MQ=255)
                                                                        gTATCCAGCCCCCCGCCGTGCAAGCGCCGACCGGTCTACCTAGCGACGCTAACCACCTGACCCAGCGCAGCTGTACCGTGGCCTTTGGCGGTCAGGCTGTTTTTCACCACTTCTCCGTTATATTCAGTGATCCGGTTAATTGCCTGATa                                                                            <  6:355995/149‑1 (MQ=255)
                                                                        gTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTa                                                                                                                            <  5:96782/101‑1 (MQ=255)
                                                                        gTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTa                                                                                                                            >  6:96782/1‑101 (MQ=255)
                                                                        gTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATa                                                                            <  5:129362/149‑1 (MQ=255)
                                                                             cAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACgg                                                                       >  6:64203/1‑149 (MQ=255)
                                                                             cAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACgg                                                                       >  6:121758/1‑149 (MQ=255)
                                                                                gCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCAt                                                                    >  2:56635/1‑149 (MQ=255)
                                                                                 ccGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATc                                                                   <  8:88900/149‑1 (MQ=255)
                                                                                   gACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGcc                                                                                 >  8:498157/1‑133 (MQ=255)
                                                                                   gACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGcc                                                                                 <  7:498157/133‑1 (MQ=255)
                                                                                     cGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGCCGATATCCACCTGACCCAGCGCATCTTTACCGTGGCGTTTGGCGGTCAGGCTGTATTTCCCCAGTTGGACGTTATATTCGGTGATGCGGTTAATTGCCTGATAGCCGGCATCgccc                                                               >  4:328506/1‑149 (MQ=255)
                                                                                         gTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCgaccggac                                                           >  4:162868/1‑149 (MQ=255)
                                                                                          tGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCgaccggacc                                                          >  5:363288/1‑149 (MQ=255)
                                                                                                  gcgACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGaccgttacc                                                    <  1:450001/147‑1 (MQ=255)
                                                                                                  gcgACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGaccgttacc                                                    >  2:450001/1‑147 (MQ=255)
                                                                                                   cgACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGCTGCGGTTAATTGCCTGATAGACGGCATCGACCGGaccgttaccgtt                                                 <  6:363288/149‑1 (MQ=255)
                                                                                                        gTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCCGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCgacc                                                               >  5:452893/1‑130 (MQ=255)
                                                                                                        gTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCCGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCgacc                                                               <  6:452893/130‑1 (MQ=255)
                                                                                                         ttGTAGTTAGCGACGATATCCACCTGACCCAGCGCATTTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCgac                                                                <  7:327977/128‑1 (MQ=255)
                                                                                                         ttGTAGTTAGCGACGATATCCACCTGACCCAGCGCATTTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCgac                                                                >  8:327977/1‑128 (MQ=255)
                                                                                                         ttGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTggcggc                                           >  3:460074/1‑149 (MQ=255)
                                                                                                            tAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGtt                                                                                                          >  3:278626/1‑83 (MQ=255)
                                                                                                            tAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGtt                                                                                                          <  4:278626/83‑1 (MQ=255)
                                                                                                            tAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGaccgtta                                                      >  1:187334/1‑135 (MQ=255)
                                                                                                            tAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGaccgtta                                                      <  2:187334/135‑1 (MQ=255)
                                                                                                             aGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGcc                                                                                 <  3:463051/107‑1 (MQ=255)
                                                                                                             aGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGcc                                                                                 >  4:463051/1‑107 (MQ=255)
                                                                                                              gTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGGTTTATGCTGGGATGGGGTTATTGGCGGGATTGAGGGCTGCGCCCGGCCGTTTACGGTTggcgggtgatg                                      >  3:126605/1‑143 (MQ=255)
                                                                                                               ttAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGTCATCGACCGGACCGTTACCGTTGGCGGCTTCTGc                                     <  2:384657/149‑1 (MQ=255)
                                                                                                                 aGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCtt                                   >  8:447513/1‑149 (MQ=255)
                                                                                                                 aGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCt                                    <  2:283982/148‑1 (MQ=255)
                                                                                                                 aGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCt                                    <  7:405240/148‑1 (MQ=255)
                                                                                                                  gCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCttt                                  >  7:370376/1‑149 (MQ=255)
                                                                                                                  gCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTGCTGCttt                                  <  6:422146/149‑1 (MQ=255)
                                                                                                                           tCCACCTGACCCAGCGGATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCTCCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGTCGGCATCGACCGGACCGTTAGCGTTGGCGGCTTCTGCTTTGACttcttc                         <  4:338265/149‑1 (MQ=255)
                                                                                                                           tCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACttcttc                         <  3:4072/149‑1 (MQ=255)
                                                                                                                             cACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGc                       <  1:33695/149‑1 (MQ=255)
                                                                                                                             cACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCAACAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCg                        <  8:206045/148‑1 (MQ=255)
                                                                                                                              aCCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGcc                      >  7:474788/1‑149 (MQ=255)
                                                                                                                                   aCCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGCcaca                   <  5:121758/147‑1 (MQ=255)
                                                                                                                                   aCCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGCCACAgg                 <  8:370376/149‑1 (MQ=255)
                                                                                                                                     ccAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGCCACAGGc                >  2:460422/1‑148 (MQ=255)
                                                                                                                                     ccAGCGCATCTTTACAGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGCCACAGGc                >  6:146798/1‑148 (MQ=255)
                                                                                                                                      cAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGCCACAGGc                >  2:111140/1‑147 (MQ=255)
                                                                                                                                        gcgcATCTTTACAGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGCCACAGGCCAGt            <  3:79985/149‑1 (MQ=255)
                                                                                                                                           cATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGCCACAGGCCAGttt          >  3:305929/1‑148 (MQ=255)
                                                                                                                                               tttACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGCCACAGGCCAGTTTGAcgg     >  6:410474/1‑149 (MQ=255)
                                                                                                                                               tttACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGCCACAGGCCAGTTTGAcgg     >  2:286732/1‑149 (MQ=255)
                                                                                                                                                ttACAGTGGCCTTTGGCGGTCAGGCTGAATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGTAGGATTCTGCTTTGACTTCTTCGCCACAGGCCAGTTTGAcggc    >  6:161638/1‑149 (MQ=255)
                                                                                                                                                 tACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGCCACAGGCCAGTTTGAcg      <  5:64203/146‑1 (MQ=255)
                                                                                                                                                 tACAGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGCCACAGGCCAGTTTGAc       <  8:93478/145‑1 (MQ=255)
                                                                                                                                                  aCCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGCCACAGGCCAGTTTGAcggcgg  >  4:295093/1‑149 (MQ=255)
                                                                                                                                                  aCCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGCCACAGGCCAGTTTGAcggcgg  >  6:159433/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCGACTTCTGCGGCACGCCAGATATTGTTCAGAACGTGCACCATGGCTTTGGCAGATGACTCGACAATATCGGTAGCCAGGCCGACGCCGTGGAAGCGGCGACCGTTGTAGTTAGCGACGATATCCACCTGACCCAGCGCATCTTTACCGTGGCCTTTGGCGGTCAGGCTGTATTTCACCAGTTCGACGTTATATTCAGTGATGCGGTTAATTGCCTGATAGACGGCATCGACCGGACCGTTACCGTTGGCGGCTTCTGCTTTGACTTCTTCGCCACAGGCCAGTTTGACGGCGG  >  NC_000913/82013‑82307

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