Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,983,452 A→T 5.9% V30E (GTG→GAG)  araH ← L‑arabinose ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,983,4520AT 5.9% 210.7 / 3.8 68V30E (GTG→GAG) araHL‑arabinose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (29/35);  new base T (2/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

TCAAAGTCACCGGAAGCGAGGCAGAACAACATGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGACTTTAGGCATCGCAAGGCTCA  >  NC_000913/1983308‑1983595
                                                                                                                                                |                                                                                                                                               
tCAAAGTCACCGGAAGCGAGGCAGAACAACATGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAAcaccac                                                                                                                                             >  7:224708/1‑149 (MQ=255)
tCAAAGTCACCGGAAGCGAGGCAGAACAACATGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAAcaccac                                                                                                                                             <  6:340813/149‑1 (MQ=255)
   aaGTCACCGGAAGCGAGGCAGAACAACATGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAg                                                                                                                                          >  5:21335/1‑149 (MQ=255)
    aGTCACCGGAAGCGAGGCAGAACAACATGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGc                                                                                                                                         <  4:317673/149‑1 (MQ=255)
      tCACCGGAAGCGAGGCAGAACAACATGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAGGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCAt                                                                                                                                       <  1:373686/149‑1 (MQ=255)
          cGGAAGCGAGGCAGAACAACATGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCg                                                                                                                                   >  1:89276/1‑149 (MQ=255)
                  aGGCAGAACAACATGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACTCCACCAGc                                                                                                                                         >  7:182154/1‑135 (MQ=255)
                  aGGCAGAACAACATGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACTCCACCAGc                                                                                                                                         <  8:182154/135‑1 (MQ=255)
                  aGGCAGAACAACATGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACTCCACCAGc                                                                                                                                         <  8:187777/135‑1 (MQ=255)
                  aGGCAGAACAACATGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACTCCACCAGc                                                                                                                                         >  7:187777/1‑135 (MQ=255)
                       gAACAACATGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGa                                                                                                                      <  3:8407/149‑1 (MQ=255)
                       gAACAACATGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGa                                                                                                                      <  8:401498/149‑1 (MQ=255)
                            acaTGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCGAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCTAACACCACCAGCATGCCGTACTGATCCCAGATACGc                                                                                                                 <  8:144338/149‑1 (MQ=255)
                             caTGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGcc                                                                                                                >  5:149948/1‑149 (MQ=255)
                              aTGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGccc                                                                                                               <  6:21335/149‑1 (MQ=255)
                                  cacaAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAg                                                                                                           <  7:23998/149‑1 (MQ=255)
                                  cacaAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAg                                                                                                           >  5:76382/1‑149 (MQ=255)
                                   acaAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAAACCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGc                                                                                                          <  2:25872/149‑1 (MQ=255)
                                       cccaccaTCCCCGACATGGAAATTGCCAGCCCCCACCCTTTCATATTAATGAAGCTCGCAAAAATAGGGACAAAAATGGCACAGCCGATAAATAGAAAAGCAAAAACCACCAGCAACCCGAACTGATCCCAGATCCGACCGAAGCTGa                                                                                                       >  2:272410/2‑148 (MQ=255)
                                         caccaTCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAAt                                                                                                     >  5:378239/1‑148 (MQ=255)
                                         caccaTCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATg                                                                                                    <  8:224708/149‑1 (MQ=255)
                                              tCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATgacgac                                                                                               >  8:220608/1‑149 (MQ=255)
                                                 ccGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACt                                                                                              >  7:57109/1‑147 (MQ=255)
                                                  cGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAg                                                                                           >  3:123149/1‑149 (MQ=255)
                                                       tGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTgcg                                                                                      >  2:114714/1‑149 (MQ=255)
                                                         gAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTgcgc                                                                                     >  7:413609/1‑148 (MQ=255)
                                                            aTTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGa                                                                                 <  2:89276/149‑1 (MQ=255)
                                                                ccAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGAcccc                                                                             >  1:182338/1‑149 (MQ=255)
                                                                ccAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGAccc                                                                              <  8:413609/148‑1 (MQ=255)
                                                                ccAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGAccc                                                                              <  6:378239/148‑1 (MQ=255)
                                                                 cAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCg                                                                            <  6:149948/149‑1 (MQ=255)
                                                                   ggCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGAt                                                                          <  4:123149/149‑1 (MQ=255)
                                                                       cAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATg                                                                         >  2:124649/1‑146 (MQ=255)
                                                                        aaCCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGa                                                                     <  1:390460/149‑1 (MQ=255)
                                                                          cccTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGaaa                                                                   <  7:220608/149‑1 (MQ=255)
                                                                             tttCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGAGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGaaa                                                                   >  3:139186/1‑146 (MQ=255)
                                                                               tCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACag                                                                >  6:219855/1‑147 (MQ=255)
                                                                                cATAATAATGAAGGGGGCAAAATTTGGGACCAAAATGGCACAGGCGATAAAGACCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACagaa                                                              <  2:81080/148‑1 (MQ=255)
                                                                                   atTAATGAAGGTGTCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACagaaga                                                            <  6:197550/147‑1 (MQ=255)
                                                                                   atTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACagaaga                                                            >  5:197550/1‑147 (MQ=255)
                                                                                   atTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGAca                                                          >  3:163496/1‑149 (MQ=255)
                                                                                    ttAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGAcat                                                         >  7:341462/1‑149 (MQ=255)
                                                                                        tGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATa                                                     <  1:124649/149‑1 (MQ=255)
                                                                                          aaGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATAc                                                    <  2:182338/148‑1 (MQ=255)
                                                                                            ggTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATActct                                                 <  4:163496/149‑1 (MQ=255)
                                                                                              tGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACtctc                                                <  4:139186/148‑1 (MQ=255)
                                                                                                  aaaaTTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTAc                                           <  3:50920/149‑1 (MQ=255)
                                                                                                  aaaaTTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTAc                                           >  4:6712/1‑149 (MQ=255)
                                                                                                  aaaaTTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTAc                                           >  1:342577/1‑149 (MQ=255)
                                                                                                         gggACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGc                                     <  1:114714/148‑1 (MQ=255)
                                                                                                                aaaTGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCt                                              >  4:169848/1‑132 (MQ=255)
                                                                                                                aaaTGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCt                                              <  3:169848/132‑1 (MQ=255)
                                                                                                                aaaTGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCt                             <  8:341462/149‑1 (MQ=255)
                                                                                                                  aTGCCACAGGCGATAAAGAGCACCGCAAACACCACCCCCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGCCCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTgg                           <  2:342577/149‑1 (MQ=255)
                                                                                                                      cacaGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGa                       <  3:6712/149‑1 (MQ=255)
                                                                                                                        caGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTg                            <  5:85196/142‑1 (MQ=255)
                                                                                                                        caGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACATAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTg                            >  6:85196/1‑142 (MQ=255)
                                                                                                                         aGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGACtt                    >  4:59014/1‑149 (MQ=255)
                                                                                                                           gCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGCATGACGACTTAGGTGCGCCCGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGACtg                    >  1:363615/1‑146 (MQ=255)
                                                                                                                           gCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGACttt                   >  7:292193/1‑148 (MQ=255)
                                                                                                                           gCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGACttt                   <  8:292193/148‑1 (MQ=255)
                                                                                                                           gCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGACTTTa                  <  1:272410/149‑1 (MQ=255)
                                                                                                                             gATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGACTTTAgg                >  8:159871/1‑149 (MQ=255)
                                                                                                                                     gCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGACTTTAGGCATCGCaa        <  8:57109/149‑1 (MQ=255)
                                                                                                                                         cGCAAACACCACCAGCATGCCGTACTGATACCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGCCCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGACTTTAGGCATAGCAAGGCt    >  1:110467/1‑149 (MQ=255)
                                                                                                                                          gCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGACTTTa                  <  1:286242/134‑1 (MQ=255)
                                                                                                                                          gCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGACTTTa                  >  2:286242/1‑134 (MQ=255)
                                                                                                                                           cAAACACCCCCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGGGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGACTTTAGGCATCGCAAGGCTCa  <  6:76382/149‑1 (MQ=255)
                                                                                                                                           cAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGACTTTAGGCATCGCAAGGCt    >  1:86665/1‑147 (MQ=255)
                                                                                                                                           cAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGACTTTAGGCATCGCAAGGCTCa  <  5:219855/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TCAAAGTCACCGGAAGCGAGGCAGAACAACATGCCACAAGCCACCATCCCCGACATGGAAATTGCCAGGCCCAACCCTTTCATATTAATGAAGGTGGCAAAATTTGGGACAAAAATGGCACAGGCGATAAAGAGCACCGCAAACACCACCAGCATGCCGTACTGATCCCAGATACGCCCGAAGCTGAATGACGACTTAGGTGCGCCAGACCCCGATGTAGAAACAGAAGACATCATACTCTCCTTACTCAGGCAACAGCCTGGCTGACTTTAGGCATCGCAAGGCTCA  >  NC_000913/1983308‑1983595

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: