Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,521,083 A→G 5.2% L313P (CTG→CCG)  yjhF ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,521,0830AG 5.2% 245.3 / 2.9 77L313P (CTG→CCG) yjhFputative transporter
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base G (2/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

ACCCAGACAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATC  >  NC_000913/4520935‑4521226
                                                                                                                                                    |                                                                                                                                               
aCCCAGACAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATcagca                                                                                                                                                 <  3:142900/149‑1 (MQ=255)
  ccAGACAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATcagcagc                                                                                                                                               >  7:341313/1‑149 (MQ=255)
   cAGACACATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATcagcagca                                                                                                                                              <  2:102619/149‑1 (MQ=255)
    aGACAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATcagcagc                                                                                                                                               <  6:370453/147‑1 (MQ=255)
    aGACAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATcagcagc                                                                                                                                               >  5:370453/1‑147 (MQ=255)
        aaaTACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCg                                                                                                                                         >  7:358389/1‑149 (MQ=255)
                  gAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCcg                                                                                                                               <  7:138966/149‑1 (MQ=255)
                     tGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGa                                                                                                                            >  3:345443/1‑149 (MQ=255)
                         gCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGc                                                                                                                        <  8:358389/149‑1 (MQ=255)
                               ccccAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTcc                                                                                                                  <  5:331519/149‑1 (MQ=255)
                                 ccAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAAATTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCg                                                                                                                 >  5:388128/1‑148 (MQ=255)
                                    ggCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCCCCACCACCACCGATAATCAGCAGCAGACCGGCAAATCCCCCGATAGCTTTTCCGc                                                                                                                >  5:247475/1‑146 (MQ=255)
                                     gCCATCAGGATCTGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGa                                                                                                            <  1:244568/149‑1 (MQ=255)
                                     gCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCGGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAg                                                                                                             >  3:154996/1‑148 (MQ=255)
                                     gCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCGGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAg                                                                                                             >  3:160676/1‑148 (MQ=255)
                                             gATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCCCCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCAt                                                                                                      >  2:214448/1‑147 (MQ=255)
                                              aTCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGc                                                                                                                                       >  7:122889/1‑113 (MQ=255)
                                              aTCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGc                                                                                                                                       <  8:122889/113‑1 (MQ=255)
                                                cGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCAtcaggt                                                                                                 >  7:14803/1‑149 (MQ=255)
                                                 gggTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAAACCAtcaggt                                                                                                 >  8:101096/1‑148 (MQ=255)
                                                   gTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCCCCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCAtcaggtcag                                                                                              >  2:293409/1‑149 (MQ=255)
                                                   gTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCAtcaggtcag                                                                                              >  6:397495/1‑149 (MQ=255)
                                                        tatCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCAtcaggtcagg                                                                                             >  3:177536/1‑145 (MQ=255)
                                                        tatCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCAtcaggtcagg                                                                                             <  4:177536/145‑1 (MQ=255)
                                                        tatCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGAtt                                                                                         <  5:228270/149‑1 (MQ=255)
                                                         atCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACAACAACCGATAATCAGCAGCAGACCGGCAATACCCCCGATAGCTTTTCAGCAAGAATACATCAGGTCAGGGAttc                                                                                        >  2:99268/1‑148 (MQ=255)
                                                          tCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGt                                                                                       >  7:335370/1‑149 (MQ=255)
                                                             aTACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGt                                                                                       >  5:359549/1‑146 (MQ=255)
                                                             aTACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGt                                                                                       <  6:359549/146‑1 (MQ=255)
                                                                cccGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCGGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTccc                                                                                 <  4:154996/149‑1 (MQ=255)
                                                                cccGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCGGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTccc                                                                                 <  4:160676/149‑1 (MQ=255)
                                                                cccGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTccc                                                                                 <  2:5377/149‑1 (MQ=255)
                                                                 ccGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCcg                                                                                >  2:387617/1‑149 (MQ=255)
                                                                 ccGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCcg                                                                                >  1:400987/1‑149 (MQ=255)
                                                                   gAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCcgcg                                                                              >  8:96294/1‑149 (MQ=255)
                                                                          gggTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAgg                                                                       <  6:63372/149‑1 (MQ=255)
                                                                            gTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCGGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGcc                                                                     >  4:135128/1‑149 (MQ=255)
                                                                            gTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGcc                                                                     >  8:196596/1‑149 (MQ=255)
                                                                                aaaTATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATg                                                                 <  8:341313/149‑1 (MQ=255)
                                                                                   tataCTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCATCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGta                                                              <  1:99268/149‑1 (MQ=255)
                                                                                   tataCTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGt                                                               >  4:217323/1‑148 (MQ=255)
                                                                                     taCTGACCAACGCCAGAGTCGATGAGAACTTGCTTGACTGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGtata                                                            <  5:187820/149‑1 (MQ=255)
                                                                                     taCTGACCAACGACAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCTCCCCAGGCCCATGGtat                                                             >  2:352072/1‑148 (MQ=255)
                                                                                         gACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAAGCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATc                                                        >  8:362295/1‑149 (MQ=255)
                                                                                             aaCGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATcgcg                                                     >  3:69219/1‑148 (MQ=255)
                                                                                              aCGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCAGCAAGAATCCATCAGGTCAGGGATTGTTTTCACCCGCCCCAGGCCCATGGTATAGATAGCGaa                                                   >  1:209942/1‑149 (MQ=255)
                                                                                                   agagTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCa                                              <  7:101096/149‑1 (MQ=255)
                                                                                                           tGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATTGTATAGATCGCGAATAGCAGTGAAATc                                      <  7:196596/149‑1 (MQ=255)
                                                                                                            gagaACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCCGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATcgc                                                      <  1:282983/132‑1 (MQ=255)
                                                                                                            gagaACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCATGGATTGTTTTCCCGCGCCCACGGCCCATGGTATAGATcgc                                                      >  2:282983/1‑132 (MQ=255)
                                                                                                             agaACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCa                                     <  7:96294/148‑1 (MQ=255)
                                                                                                             agaACTTGCTTGAATGCAACACCACCACAGATAATCATCAGCAGACCGGCAATACCCGAGATAGCTTTTCCGCAAGAAt                                                                                                          <  6:174299/79‑1 (MQ=255)
                                                                                                             agaACTTGCTTGAATGAACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAAt                                                                                                          >  5:174299/1‑79 (MQ=255)
                                                                                                              gaACTTGCTTGAATTCAACACAACCACCGATAATCAGCAGAAGACCGGATATACGCGCGATAGCTTTTCAGCAAGAATCCATCAGGTTAGGGATTGTTTTCCCGCGCCCAAGGCCCATTGTATAGATCGCTAATAGCAGTGAAATCAGc                                   <  4:345443/149‑1 (MQ=255)
                                                                                                               aaCTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCa                                  <  5:397495/149‑1 (MQ=255)
                                                                                                                 cttgcttgAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATg                                >  6:235685/1‑149 (MQ=255)
                                                                                                                  ttgcttgAATGCACCACCCCCAGAGATAATCAGCAGCAGACCGGCAATACCGGCGATAGCTTTTCCGCAAGAATCCATCTGGTCAGGGTTTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCa                                  <  6:167986/146‑1 (MQ=255)
                                                                                                                  ttgcttgAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCa                                              >  2:66127/1‑134 (MQ=255)
                                                                                                                  ttgcttgAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCa                                              <  1:66127/134‑1 (MQ=255)
                                                                                                                  ttgcttgAATGCACCACAACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCAACAGGTCAGGCATTGTTTTACCCCGCCCCCCGCCCATGGAATAGATCCCCAAAAGCAGTGAAAACAGCa                                  >  5:167986/1‑146 (MQ=255)
                                                                                                                      ttgAATGCACCACCACCACCGATAATCAGCAGCGGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGAt                           >  4:346025/1‑149 (MQ=255)
                                                                                                                         aaTGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATg                                                                 <  1:242973/108‑1 (MQ=255)
                                                                                                                         aaTGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATg                                                                 >  2:242973/1‑108 (MQ=255)
                                                                                                                          aTGCACCACCACCACCGATAATCAGCAGCGGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGtt                       <  3:135128/149‑1 (MQ=255)
                                                                                                                          aTGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGtt                       <  4:69219/149‑1 (MQ=255)
                                                                                                                          aTGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGtt                       <  6:388128/149‑1 (MQ=255)
                                                                                                                          aTGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGtt                       <  1:387617/149‑1 (MQ=255)
                                                                                                                            gCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTg                      >  2:250343/1‑148 (MQ=255)
                                                                                                                             caccagcagcacCGATAATCAGCAGCAGACCGGCAATACGCGCGATAGCTTTTCCTCAAGAATCCATCAGGTCAGTGATTGTTTTCCCGCGCCCCAGGGCCATGGTATAGATCGCGAATAGGAGTGAAATCAGCATGGCGATCGTTGaa                    <  8:14803/149‑1 (MQ=255)
                                                                                                                             caccaccaccacCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGa                     <  2:400987/148‑1 (MQ=255)
                                                                                                                              acctccaccaccGATAATCAGCAGCTTTCGGGCAATACCGGCGATAGCTTTTCCGCTTGTATCCATCAGGTGAGGGTTTGTTTTCCCTCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAAt                   <  5:238340/149‑1 (MQ=255)
                                                                                                                              accaccaccaccGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAg                                                                        >  7:163253/1‑96 (MQ=255)
                                                                                                                              accaccaccaccGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAg                                                                        <  8:163253/96‑1 (MQ=255)
                                                                                                                              accaccaccaccGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATTGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAAt                   <  1:352072/149‑1 (MQ=255)
                                                                                                                                  ccaccaccGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTAcc               >  3:376770/1‑149 (MQ=255)
                                                                                                                                   aaccaacGATAATCAGCAGCAGCCAGGCAATACCCGAGATAGCATATACGCAAGAATCCAACAGGTCAGCGATAGATTTCCCGCGCACAAGCCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGATGAAATAcc               >  2:337736/2‑148 (MQ=255)
                                                                                                                                    accaccGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTa             <  7:362295/149‑1 (MQ=255)
                                                                                                                                     ccaccGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTaa            <  2:242217/149‑1 (MQ=255)
                                                                                                                                      caccGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAg           <  1:214448/149‑1 (MQ=255)
                                                                                                                                       cccGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAg           >  5:272477/2‑148 (MQ=255)
                                                                                                                                         cGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAAc        <  1:293409/149‑1 (MQ=255)
                                                                                                                                               tCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATc  <  1:250343/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ACCCAGACAAATACGCAGGAATGCGGCCACGCCCCAGGCCATCAGGATCGGGTTAATATCCATACCCGAAACCAGGGTGGAAATATACTGACCAACGCCAGAGTCGATGAGAACTTGCTTGAATGCACCACCACCACCGATAATCAGCAGCAGACCGGCAATACCCGCGATAGCTTTTCCGCAAGAATCCATCAGGTCAGGGATTGTTTTCCCGCGCCCCAGGCCCATGGTATAGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATC  >  NC_000913/4520935‑4521226

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