Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 156,157 G→A 5.4% T15T (ACC→ACT)  yadV ← putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913156,1570GA 5.4% 252.6 / 3.4 74T15T (ACC→ACT) yadVputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base G (29/41);  new base A (2/2);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

CGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAATACTTA  >  NC_000913/156009‑156301
                                                                                                                                                    |                                                                                                                                                
cGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAg                                                                                                                                                  <  1:702638/149‑1 (MQ=255)
 gCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAgg                                                                                                                                                 >  4:243429/1‑149 (MQ=255)
    aGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAAGTTa                                                                                                                                              >  1:116459/1‑149 (MQ=255)
    aGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAAGTTa                                                                                                                                              >  1:767517/1‑149 (MQ=255)
       aTCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTT‑ATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGaa                                                                                                                                          <  4:91849/149‑1 (MQ=255)
          tAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGc                                                                                                                                        >  1:448439/1‑149 (MQ=255)
            aCCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCAt                                                                                                                                      >  2:12483/1‑149 (MQ=255)
            aCCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCAt                                                                                                                                      >  2:384996/1‑149 (MQ=255)
               aaCTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATaaa                                                                                                                                   >  4:801450/1‑149 (MQ=255)
                aCTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAg                                                                                                                                  <  4:393705/149‑1 (MQ=255)
                  tctGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTa                                                                                                                                              <  3:745686/135‑1 (MQ=255)
                  tctGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTa                                                                                                                                              >  4:745686/1‑135 (MQ=255)
                     ggACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCtgtt                                                                                                                             <  4:264784/149‑1 (MQ=255)
                     ggACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCtgtt                                                                                                                             <  1:384996/149‑1 (MQ=255)
                      gACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCtgtt                                                                                                                             <  1:48261/148‑1 (MQ=255)
                      gACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCtgtt                                                                                                                             >  2:48261/1‑148 (MQ=255)
                           gCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGt                                                                                                                       <  2:436214/149‑1 (MQ=255)
                                gggTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGgt                                                                                                                  <  3:243429/149‑1 (MQ=255)
                                 ggTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGgtg                                                                                                                 <  1:12483/149‑1 (MQ=255)
                                  gTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTAc                                                                                                                                             <  3:105481/120‑1 (MQ=255)
                                  gTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTAc                                                                                                                                             >  4:105481/1‑120 (MQ=255)
                                      ttCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTaaa                                                                                                            >  2:54455/1‑149 (MQ=255)
                                          cTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGttt                                                                                                                     >  4:665211/1‑136 (MQ=255)
                                          cTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGttt                                                                                                                     <  3:665211/136‑1 (MQ=255)
                                           tttATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAg                                                                                                                                         <  4:703950/115‑1 (MQ=255)
                                           tttATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAg                                                                                                                                         >  3:703950/1‑115 (MQ=255)
                                           tttATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACa                                                                                                       <  1:281621/149‑1 (MQ=255)
                                              aTTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATaa                                                                                                    >  1:631786/1‑149 (MQ=255)
                                               ttttCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATaaa                                                                                                   <  3:515834/149‑1 (MQ=255)
                                                     aGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATcc                                                                                             <  3:801450/149‑1 (MQ=255)
                                                     aGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATcc                                                                                             >  3:681897/1‑149 (MQ=255)
                                                      gACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATccc                                                                                            <  1:735302/149‑1 (MQ=255)
                                                         gtacgtTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGa                                                                                         <  3:661725/149‑1 (MQ=255)
                                                         gtacgtTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGa                                                                                         <  3:93618/149‑1 (MQ=255)
                                                             gtTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATaga                                                                                     <  2:631786/149‑1 (MQ=255)
                                                               tGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGt                                                                                   <  4:413422/149‑1 (MQ=255)
                                                                  cacTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGttaa                                                                                <  2:113984/149‑1 (MQ=255)
                                                                    cTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGttaatt                                                                              <  1:54455/149‑1 (MQ=255)
                                                                       ttttGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTaaa                                                                           >  2:702591/1‑149 (MQ=255)
                                                                       ttttGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTaaa                                                                           >  2:370546/1‑149 (MQ=255)
                                                                       ttttGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTaaa                                                                           <  2:10186/149‑1 (MQ=255)
                                                                        tttGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCCTCCCTGATAGAGTTAATTaaaa                                                                          <  3:397032/149‑1 (MQ=255)
                                                                        tttGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTaaaa                                                                          <  3:57982/149‑1 (MQ=255)
                                                                           gATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACac                                                                       <  2:690131/149‑2 (MQ=255)
                                                                            aTCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAGCAgg                                                                      <  4:560396/149‑1 (MQ=255)
                                                                            aTCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAgg                                                                      >  4:343532/1‑149 (MQ=255)
                                                                            aTCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAgg                                                                      >  3:560396/1‑149 (MQ=255)
                                                                               gCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGt                                                                                                                                                >  2:97752/1‑72 (MQ=255)
                                                                               gCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGt                                                                                                                                                <  1:97752/72‑1 (MQ=255)
                                                                                   ttATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAAGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGaataat                                                               <  2:767517/149‑1 (MQ=255)
                                                                                   ttATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAAGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGaataat                                                               <  2:116459/149‑1 (MQ=255)
                                                                                     atatatTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAatat                                                             <  3:685349/149‑1 (MQ=255)
                                                                                        tatTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATaaa                                                          <  2:448439/149‑1 (MQ=255)
                                                                                                    aCCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAattatt                                              <  2:411491/149‑1 (MQ=255)
                                                                                                      ccGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAattattat                                            >  3:89566/1‑149 (MQ=255)
                                                                                                      ccGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAattattat                                            >  1:148041/1‑149 (MQ=255)
                                                                                                            tGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTg                                      <  1:370546/149‑1 (MQ=255)
                                                                                                            tGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTg                                      <  1:702591/149‑1 (MQ=255)
                                                                                                                       gCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGaataat                                                               >  4:634854/1‑113 (MQ=255)
                                                                                                                       gCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGaataat                                                               <  3:634854/113‑1 (MQ=255)
                                                                                                                         aaTAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGt                                     <  1:746874/137‑1 (MQ=255)
                                                                                                                         aaTAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGt                                     >  2:746874/1‑137 (MQ=255)
                                                                                                                             gACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATtcaatcaa                     <  3:422361/149‑1 (MQ=255)
                                                                                                                                 aTGAACTAAAAGCCATACAGGTTACCAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGTATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTa                 >  3:535657/1‑149 (MQ=255)
                                                                                                                                  tGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTaa                <  4:89566/149‑1 (MQ=255)
                                                                                                                                   gAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAAc               >  1:379685/1‑149 (MQ=255)
                                                                                                                                    aaCTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACt              <  2:466867/149‑1 (MQ=255)
                                                                                                                                       tAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATattc                               <  1:407323/129‑1 (MQ=255)
                                                                                                                                       tAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATattc                               >  2:407323/1‑129 (MQ=255)
                                                                                                                                       tAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGt           >  4:160283/1‑149 (MQ=255)
                                                                                                                                       tAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGt           >  3:801911/1‑149 (MQ=255)
                                                                                                                                         aaaGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGCAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGa         >  3:509332/1‑149 (MQ=255)
                                                                                                                                               aTACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAAtactt   >  2:553349/1‑149 (MQ=255)
                                                                                                                                                tACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAAttt                  <  2:76628/133‑1 (MQ=255)
                                                                                                                                                tACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAAttt                  >  1:76628/1‑133 (MQ=255)
                                                                                                                                                tACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAAtactta  <  2:148041/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAATACTTA  >  NC_000913/156009‑156301

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