Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 452,102 G→A 7.2% R7H (CGC→CAC)  aruC → acetylornithine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947452,1020GA 7.2% 204.1 / 3.3 69R7H (CGC→CAC) aruCacetylornithine aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base A (3/2);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATGATTCGGTCGATTGTTACCCCGGATGAAAAACACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGA  >  NC_002947/451961‑452249
                                                                                                                                             |                                                                                                                                                   
gCGATGATGCGGTCGATTGTTACCCAGGATGAAAAACCCTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCg                                                                                                                                              <  8:58926/149‑1 (MQ=255)
         cGGTCGATTGTTACCCCGGATGAAAAACACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTgcc                                                                                                                                     >  8:187647/1‑149 (MQ=255)
           gTCGATTGTTACCCCGGATGAAAAACACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTgccgc                                                                                                                                   >  6:272243/1‑149 (MQ=255)
                ttGTTACCCCGGATGAAAAACACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTgccgccgccg                                                                                                                              >  5:157882/1‑149 (MQ=255)
                ttGTTACCCCGGATGAAAAACACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGACACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCACCGCTCGGCTCCTgccgccgccg                                                                                                                              >  6:375657/1‑149 (MQ=255)
                      ccccGGATGAAAAACACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGc                                                                                                                        <  2:128714/149‑1 (MQ=255)
                        ccGGATGAAAAACACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGcc                                                                                                                       >  8:159359/1‑148 (MQ=255)
                         cGGATGAAAAACACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGccc                                                                                                                      >  6:50920/1‑148 (MQ=255)
                         cGGATGAAAAACACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGccc                                                                                                                      >  1:149354/1‑148 (MQ=255)
                            aTGAAAAACACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGc                                                                                                                  <  6:302203/149‑1 (MQ=255)
                            aTGAAAAACACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGc                                                                                                                  <  7:187647/149‑1 (MQ=255)
                                 aaaCACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTgcg                                                                                                             >  4:110943/1‑149 (MQ=255)
                                 aaaCACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTgcg                                                                                                             >  4:15699/1‑149 (MQ=255)
                                 aaaCACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTgcg                                                                                                             >  7:295466/1‑149 (MQ=255)
                                  aaCACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTgcgc                                                                                                            <  7:121799/149‑1 (MQ=255)
                                     acTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCgg                                                                                                         <  8:305396/149‑1 (MQ=255)
                                     acTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCgg                                                                                                         <  5:272243/149‑1 (MQ=255)
                                        tttGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGa                                                                                                      >  7:267797/1‑149 (MQ=255)
                                         ttGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGa                                                                                                      <  1:75487/148‑1 (MQ=255)
                                                gcaaggcaaAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCACCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCg                                                                                              >  6:336081/1‑149 (MQ=255)
                                                gcaaggcaaAATATCGAGTGCAGGCACCAAGACCGCACCGTGACTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCg                                                                                              >  1:72889/1‑149 (MQ=255)
                                                 caaggcaaAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTg                                                                                                 <  3:253763/145‑1 (MQ=255)
                                                 caaggcaaAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTg                                                                                                 >  4:253763/1‑145 (MQ=255)
                                                 caaggcaaAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGc                                                                                             <  5:50920/149‑1 (MQ=255)
                                                   aggcaaAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCAc                                                                                           <  8:38675/149‑1 (MQ=255)
                                                         aaTATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCt                                                                                     <  2:72889/149‑1 (MQ=255)
                                                          atatCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCt                                                                                     >  5:277357/1‑148 (MQ=255)
                                                             tCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCaa                                                                                                                    >  5:324265/1‑114 (MQ=255)
                                                             tCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCaa                                                                                                                    <  6:324265/114‑1 (MQ=255)
                                                             tCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAg                                                                                  <  7:183594/148‑1 (MQ=255)
                                                                 gTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTc                                                                             >  7:93787/1‑149 (MQ=255)
                                                                    cAGGCACCACGACCGCTCCGGGTCTGCCCTTTACCGTGATTGGTGGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCGCCTGCCGCCGCCGCCGAGCCCAAGCGAGCGCCGGTGATCGAGCCGGACGCGTCTGAGCAGTCGc                                                                           <  6:356254/148‑1 (MQ=255)
                                                                              gACCGCACCGTGACTGCCCTTTCCCGTGAGTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTg                                                                <  4:50715/149‑1 (MQ=255)
                                                                               aCCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGa                                                                                   <  1:239034/129‑1 (MQ=255)
                                                                               aCCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGa                                                                                   >  2:239034/1‑129 (MQ=255)
                                                                                    aCCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGATCGCCTGATGccc                                                          <  2:22070/149‑1 (MQ=255)
                                                                                                cTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGt                                              <  8:267797/149‑1 (MQ=255)
                                                                                                cTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCACCGCTCGGCTCCTGCCGCCGCCGCGGAGCCCAAGCGTGCGCCCGTGATCGTGCCGCACGCGTCTGAGCAGACGCCTGAGCGCCTGATGCCCGGCACCGAAcct                                              >  1:79825/1‑147 (MQ=255)
                                                                                                    aCCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAg                                          >  8:233266/1‑149 (MQ=255)
                                                                                                        tGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGGCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCAACGAACGTGCAGCGCa                                      <  6:157882/149‑1 (MQ=255)
                                                                                                         gATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAg                                     >  4:294460/1‑149 (MQ=255)
                                                                                                               ttggAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTg                                                                        <  7:350239/108‑1 (MQ=255)
                                                                                                               ttggAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTg                                                                        >  8:350239/1‑108 (MQ=255)
                                                                                                                 ggAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTcgt                             >  7:188623/1‑149 (MQ=255)
                                                                                                                  gAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCa                                      >  6:11783/1‑139 (MQ=255)
                                                                                                                  gAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCa                                      <  5:11783/139‑1 (MQ=255)
                                                                                                                      ccTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCGGCGCAGGTATTcgtgcgtg                        <  6:277357/149‑1 (MQ=255)
                                                                                                                            gaACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGa                                                  <  7:21745/117‑1 (MQ=255)
                                                                                                                            gaACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGa                                                  >  8:21745/1‑117 (MQ=255)
                                                                                                                            gaACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAgg                  >  7:118927/1‑149 (MQ=255)
                                                                                                                             aaCCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCa                                      >  7:322103/1‑128 (MQ=255)
                                                                                                                             aaCCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCa                                      <  8:322103/128‑1 (MQ=255)
                                                                                                                             aaCCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAgg                  <  7:233266/148‑1 (MQ=255)
                                                                                                                              aCCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGc                                                       <  7:136147/110‑1 (MQ=255)
                                                                                                                              aCCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGc                                                       >  8:136147/1‑110 (MQ=255)
                                                                                                                               ccTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCGAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAg                                     >  1:266940/1‑127 (MQ=255)
                                                                                                                               ccTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAg                                     <  2:266940/127‑1 (MQ=255)
                                                                                                                               ccTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCt               >  2:186479/1‑149 (MQ=255)
                                                                                                                                  gTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGt            >  7:242601/1‑149 (MQ=255)
                                                                                                                                   ttCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGt            <  8:242601/148‑1 (MQ=255)
                                                                                                                                   ttCAACCTGCACCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCg             <  5:336081/147‑1 (MQ=255)
                                                                                                                                   ttCAACCTGCACCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCg             <  5:375657/147‑1 (MQ=255)
                                                                                                                                    tCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTgg          >  4:246751/1‑149 (MQ=255)
                                                                                                                                     cAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTgg          <  7:159359/148‑1 (MQ=255)
                                                                                                                                      aaCCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCGGGGCTCGTGGtt        >  2:67173/1‑149 (MQ=255)
                                                                                                                                         cTGCGCCGCTCGGCTCCTGCCGCCGCCGCGGAGCCCTAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTtgtg     <  2:79825/149‑1 (MQ=255)
                                                                                                                                         cTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCa                                      >  8:131077/1‑116 (MQ=255)
                                                                                                                                         cTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCa                                      >  8:125640/1‑116 (MQ=255)
                                                                                                                                         cTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTtgtg     >  8:290275/1‑149 (MQ=255)
                                                                                                                                           gcgcCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTgg    >  3:233837/1‑148 (MQ=255)
                                                                                                                                            cgccgcTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTgg    <  3:294460/147‑1 (MQ=255)
                                                                                                                                             gccgcTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGa  >  4:207961/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GCGATGATTCGGTCGATTGTTACCCCGGATGAAAAACACTTTTGATAGGCAAGGCAAAATATCGAGTGCAGGCACCAAGACCGCACCGTGGCTGCCCTTTACCGTGATTGGTTGGAGCCCTGATGAACCTGTTCAACCTGCGCCGCTCGGCTCCTGCCGCCGCCGCCGAGCCCAAGCGTGCGCCGGTGATCGTGCCGCACGCGTCTGAGCAGTCGCCTGAGCGCCTGATGCCCGGCACCGAACGTGCAGCGCAGGTATTCGTGCGTGGCCAGGGCTCGTGGTTGTGGGA  >  NC_002947/451961‑452249

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