Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,488,983 C→A 6.2% V400L (GTG→TTG)  argD ← bifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,488,9830CA 6.2% 226.3 / 4.2 65V400L (GTG→TTG) argDbifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (32/29);  new base A (2/2);  total (34/31)
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.82e-01

GCGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTG  >  NC_000913/3488839‑3489131
                                                                                                                                                |                                                                                                                                                    
gCGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACcgcg                                                                                                                                                  >  2:531030/1‑149 (MQ=255)
 cGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGt                                                                                                                                                 <  2:235377/149‑1 (MQ=255)
 cGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGt                                                                                                                                                 <  3:94383/149‑1 (MQ=255)
  gACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTg                                                                                                                                                <  3:274267/149‑1 (MQ=255)
  gACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTg                                                                                                                                                <  3:36466/149‑1 (MQ=255)
   aCCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTgg                                                                                                                                               <  2:63032/149‑1 (MQ=255)
   aCCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTgg                                                                                                                                               <  2:142383/149‑1 (MQ=255)
   aCCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTgg                                                                                                                                               >  1:74815/1‑149 (MQ=255)
           aCGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaaggt                                                                                                                                       <  1:293925/149‑4 (MQ=255)
           aCGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACg                                                                                                                                       >  2:331336/1‑149 (MQ=255)
            cGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaaggtg                                                                                                                                      <  4:294159/149‑5 (MQ=255)
            cGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGt                                                                                                                                      <  1:531030/149‑1 (MQ=255)
             gggATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGtt                                                                                                                                     >  1:253620/1‑149 (MQ=255)
             gggATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGtt                                                                                                                                     <  1:155554/149‑1 (MQ=255)
                          gTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATc                                                                                                                        <  4:240510/149‑1 (MQ=255)
                          gTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATc                                                                                                                        <  2:74815/149‑1 (MQ=255)
                           tCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTg                                                                                                                                    >  4:167532/1‑136 (MQ=255)
                           tCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTg                                                                                                                                    <  3:167532/136‑1 (MQ=255)
                           tCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCg                                                                                                                       >  4:349190/1‑149 (MQ=255)
                             gCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGat                                                                                                                     <  4:99763/149‑1 (MQ=255)
                               aaGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGatat                                                                                                                   >  2:549990/1‑149 (MQ=255)
                                   tGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGc                                                                                                               <  4:289224/149‑1 (MQ=255)
                                                ggCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTccaccacc                                                                                                  <  1:380059/149‑1 (MQ=255)
                                                     aaTGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGa                                                                                             <  2:183993/149‑1 (MQ=255)
                                                        gCAGGCTAGACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACgg                                                                                          <  2:20432/149‑1 (MQ=255)
                                                         cAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGt                                                                                         <  4:418054/149‑1 (MQ=255)
                                                               aaaCCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCaaaa                                                                                   <  2:596044/149‑1 (MQ=255)
                                                                aaCCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAAc                                                                                  >  2:126697/1‑149 (MQ=255)
                                                                 aCCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACg                                                                                 >  3:520744/1‑149 (MQ=255)
                                                                   cATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCCACGCGGTGGCGaa                                                                                                                                          >  4:249150/1‑90 (MQ=38)
                                                                   cATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCa                                                                               >  3:38322/1‑149 (MQ=255)
                                                                    aTCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCAt                                                                              >  3:282405/1‑149 (MQ=255)
                                                                        gATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATcaca                                                                          >  1:153186/1‑149 (MQ=255)
                                                                        gATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATcaca                                                                          >  4:450080/1‑149 (MQ=255)
                                                                          tAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtc                                                                        <  3:455153/149‑1 (MQ=255)
                                                                                    aCGCCGCATCCGACATCTTTTGCCTGAGGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCa                                                              >  2:434872/1‑149 (MQ=255)
                                                                                       ccgcATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTc                                                           <  4:282405/149‑1 (MQ=255)
                                                                                        cgcATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGc                                                                                                               >  4:526101/1‑96 (MQ=255)
                                                                                        cgcATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGc                                                                                                               <  3:526101/96‑1 (MQ=255)
                                                                                               gACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCAGCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCAt                                                   >  3:482527/1‑149 (MQ=255)
                                                                                               gACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCAt                                                   >  3:298377/1‑149 (MQ=255)
                                                                                                     ttttGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtccggtccgg                                                                <  1:257049/130‑1 (MQ=255)
                                                                                                     ttttGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACAtccggtccgg                                                                >  2:257049/1‑130 (MQ=255)
                                                                                                     ttttGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGcc                                             >  3:326484/1‑149 (MQ=255)
                                                                                                           cTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCctc                                       <  3:450080/149‑1 (MQ=255)
                                                                                                           cTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCctc                                       >  3:583549/1‑149 (MQ=255)
                                                                                                            tGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCctct                                      >  3:581435/1‑149 (MQ=255)
                                                                                                            tGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCctct                                      >  1:501628/1‑149 (MQ=255)
                                                                                                                       attgaggCCCCAACCACCTTCGCCACCCCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCAtt                           <  4:506369/143‑2 (MQ=255)
                                                                                                                          tACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAg                        <  1:549990/149‑1 (MQ=255)
                                                                                                                              ccccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGaaa                    >  3:353049/1‑149 (MQ=255)
                                                                                                                               cccAACCACCTTCGCCACCGCGTGGGAGAAACGTTGCAGCCCTTCATCGATATCCGCATCTTCCACTACCAGCGACGGTGCAAAACGCATCACATCCGATCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAAt                   >  3:363501/1‑149 (MQ=255)
                                                                                                                                ccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCCTACAGGAAATc                  <  1:443452/149‑1 (MQ=255)
                                                                                                                                ccAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATc                  >  3:446769/1‑149 (MQ=255)
                                                                                                                                 cAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCa                 <  4:520744/149‑1 (MQ=255)
                                                                                                                                 cAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCa                 >  1:280128/1‑149 (MQ=255)
                                                                                                                                  aaCCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAc                >  1:598003/1‑149 (MQ=255)
                                                                                                                                   accaccTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAcg               >  4:206454/1‑149 (MQ=255)
                                                                                                                                     caccTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCAcgcg             >  4:247780/1‑149 (MQ=255)
                                                                                                                                      accTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACgcgc            >  2:375596/1‑149 (MQ=255)
                                                                                                                                        tttCGCCCACGCGGTGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGGCCGGCATTCAGCACCAGTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGccc          >  2:564984/2‑149 (MQ=255)
                                                                                                                                         ttCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCg         >  3:287751/1‑149 (MQ=255)
                                                                                                                                             ccACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCt     <  4:298377/149‑1 (MQ=255)
                                                                                                                                               aCCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCttt   >  3:240682/1‑149 (MQ=255)
                                                                                                                                               aCCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCttt   <  3:296555/149‑1 (MQ=255)
                                                                                                                                                ccGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTg  <  2:598003/149‑1 (MQ=255)
                                                                                                                                                ccGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTg  <  2:501628/149‑1 (MQ=255)
                                                                                                                                                ccGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTg  <  4:581435/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCGACCGCATCACGGGATTTGGCTGAGTCGCAAGTTGCGGGATTCGAAGGCGTAATGCAGGCTAAACCATCGGATAAGGCATTCACGCCGCATCCGACATCTTTTGCCTGATGCTTCGTATTTACGCCCCAACCACCTTCGCCACCGCGTGGGCGAAACGTTGCATCCCTTCATCGATATCCGCATCTTCCACCACCAGCGACGGTGCAAAACGCATCACATCCGGTCCGGCATTCAGCACCATTACGCCAGCCTCTGCGCCCGCATACAGGAAATCACGCGCCCGACCTTTG  >  NC_000913/3488839‑3489131

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