Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,103,106 T→A 7.5% pseudogene (284/446 nt) wbbL ← pseudogene, lipopolysaccharide biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,103,1060TA 7.5% 376.0 / 15.0 107pseudogene (284/446 nt)wbbLpseudogene, lipopolysaccharide biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base T (46/53);  new base A (3/5);  total (49/58)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 110 total aligned reads displayed.

TGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTGGCA  >  NC_000913/2102961‑2103254
                                                                                                                                                 |                                                                                                                                                    
tGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCa                                                                                                                                                   >  4:894209/1‑149 (MQ=255)
   ctcaATAAAATCAGAGAGCACAGGAAATTTTCTTACGGAATTATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGAAAATCATCATGCTTCatgg                                                                                                                                                >  2:918446/3‑148 (MQ=255)
   cacaATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatga                                                                                                                                                <  2:834264/149‑1 (MQ=255)
      aaTAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCAtgatga                                                                                                                                             >  4:421798/1‑149 (MQ=255)
      aaTAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCAtgatga                                                                                                                                             >  2:548977/1‑149 (MQ=255)
        tAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGata                                                                                                                                           <  1:548977/149‑1 (MQ=255)
         aaaaTCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgattt                                                                                                                                          <  4:182123/149‑3 (MQ=255)
         aaaaTCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgattt                                                                                                                                          <  2:924204/149‑3 (MQ=255)
           aaTCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCg                                                                                                                                        <  2:375031/149‑1 (MQ=255)
           aaTCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCg                                                                                                                                        <  3:361751/149‑1 (MQ=255)
           aaTCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCg                                                                                                                                        <  2:412295/149‑1 (MQ=255)
           aaTCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgatttgc                                                                                                                                        <  2:502354/149‑5 (MQ=255)
              cAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGa                                                                                                                                     <  2:905548/149‑1 (MQ=255)
              cAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGa                                                                                                                                     <  4:113254/149‑1 (MQ=255)
                      acaGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTTGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCaaaaa                                                                                                                             <  2:204066/149‑1 (MQ=255)
                      acaGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCaaaaa                                                                                                                             >  1:354601/1‑149 (MQ=255)
                        aGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACa                                                                                                                           >  2:271535/1‑149 (MQ=255)
                         ggAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGa                                                                                                                                     >  1:16412/1‑138 (MQ=255)
                         ggAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGa                                                                                                                                     <  2:16412/138‑1 (MQ=255)
                         ggAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACaa                                                                                                                          >  3:31610/1‑149 (MQ=255)
                          gAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACaaa                                                                                                                         <  3:840884/149‑1 (MQ=255)
                           aaaTTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACaaaa                                                                                                                        <  3:710262/149‑1 (MQ=255)
                              ttttCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGt                                                                                                                     >  2:101123/1‑149 (MQ=255)
                                ttCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa                                                                                                                   <  1:491576/149‑1 (MQ=255)
                                 tCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAat                                                                                                                  <  2:354601/149‑1 (MQ=255)
                                  cTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatc                                                                                                                 <  4:31805/149‑1 (MQ=255)
                                  cTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatc                                                                                                                 <  3:487676/149‑1 (MQ=255)
                                     aCGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatcatc                                                                                                              >  1:430762/1‑149 (MQ=255)
                                      cGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatcatca                                                                                                             <  4:466841/149‑1 (MQ=255)
                                           taatCATGTAAGGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGc                                                                                                        <  3:786165/149‑1 (MQ=255)
                                             atCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCgg                                                                                                      <  1:809575/149‑1 (MQ=255)
                                              tCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACaaa                                                                                                                         >  4:54450/1‑129 (MQ=255)
                                              tCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCggg                                                                                                     >  2:889312/1‑149 (MQ=255)
                                                   tAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCgg                                                                                                                                       <  1:30050/110‑1 (MQ=255)
                                                   tAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCgg                                                                                                                                       >  2:30050/1‑110 (MQ=255)
                                                   tAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCtat                                                                                                >  4:781953/1‑149 (MQ=255)
                                                   tAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCtat                                                                                                >  2:414161/1‑149 (MQ=255)
                                                   tAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCtat                                                                                                >  1:451468/1‑149 (MQ=255)
                                                    aaaGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCtata                                                                                               >  2:612921/1‑149 (MQ=255)
                                                     aaGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTatat                                                                                              <  3:407835/149‑1 (MQ=255)
                                                     aaGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTatat                                                                                              >  2:175477/1‑149 (MQ=255)
                                                     aaGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTatat                                                                                              <  1:889312/149‑1 (MQ=255)
                                                         tttCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATAttttt                                                                                          <  1:140683/149‑1 (MQ=255)
                                                          ttCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc                                                                                         <  4:31610/149‑1 (MQ=255)
                                                          ttCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc                                                                                         <  4:852485/149‑1 (MQ=255)
                                                          ttCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc                                                                                         <  2:530636/149‑1 (MQ=255)
                                                              cTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCttt                                                                                     <  1:735769/149‑1 (MQ=255)
                                                              cTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCttt                                                                                     <  2:912393/149‑1 (MQ=255)
                                                               ttCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTa                                                                                    <  3:706730/149‑1 (MQ=255)
                                                               ttCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTa                                                                                    <  1:612921/149‑1 (MQ=255)
                                                                    ctcGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATAc                                                                               <  4:460644/149‑1 (MQ=255)
                                                                     tcGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAGTGTAATCATCATCTGCGGGTCTATATTTTTCCTTTa                                                                                    >  1:684472/1‑143 (MQ=255)
                                                                     tcGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAGTGTAATCATCATCTGCGGGTCTATATTTTTCCTTTa                                                                                    <  2:684472/143‑1 (MQ=255)
                                                                     tcGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACg                                                                              >  2:898519/1‑149 (MQ=255)
                                                                              gCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCaa                                                                     >  1:22233/1‑149 (MQ=255)
                                                                              gCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCaa                                                                     <  2:821435/149‑1 (MQ=255)
                                                                                 tAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAatat                                                                  >  3:411050/1‑149 (MQ=255)
                                                                                    tGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattat                                                               <  2:22233/149‑1 (MQ=255)
                                                                                     gTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattatt                                                              <  4:685032/149‑1 (MQ=255)
                                                                                        cTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattattatt                                                           >  4:438042/1‑149 (MQ=255)
                                                                                         tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAat                                                                                                                  <  2:499474/93‑1 (MQ=255)
                                                                                         tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAat                                                                                                                  >  1:499474/1‑93 (MQ=255)
                                                                                         tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattattatta                                                          >  2:267721/1‑149 (MQ=255)
                                                                                         tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattattatta                                                          >  2:316914/1‑149 (MQ=255)
                                                                                           aaaGCATAACGTTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATg                                                        >  3:598253/1‑149 (MQ=255)
                                                                                           aaaGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATg                                                        >  1:391454/1‑149 (MQ=255)
                                                                                               cATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCa                                                    >  2:41172/1‑149 (MQ=255)
                                                                                                 tAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattattatt                                                           >  3:872217/1‑140 (MQ=255)
                                                                                                 tAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattattatt                                                           <  4:872217/140‑1 (MQ=255)
                                                                                                 tAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCaaa                                                  >  3:395291/1‑149 (MQ=255)
                                                                                                  aaCGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgatt                                                                                                                                           <  3:820810/59‑2 (MQ=37)
                                                                                                  aaCGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgatt                                                                                                                                           >  4:820810/1‑58 (MQ=37)
                                                                                                       ttACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTatat                                                                                              <  2:328905/99‑1 (MQ=255)
                                                                                                       ttACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTatat                                                                                              >  1:328905/1‑99 (MQ=255)
                                                                                                       ttACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGta                                            <  1:271535/149‑1 (MQ=255)
                                                                                                           tttCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATaca                                        >  3:832955/1‑149 (MQ=255)
                                                                                                           tttCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATaca                                        >  2:846560/1‑149 (MQ=255)
                                                                                                            ttCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATacac                                       <  3:611784/149‑1 (MQ=255)
                                                                                                             tCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACAcc                                      >  3:19315/1‑149 (MQ=255)
                                                                                                                aCATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTcc                                   <  2:572522/149‑1 (MQ=255)
                                                                                                                    atTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTa                               >  2:153601/1‑149 (MQ=255)
                                                                                                                         atatatGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAAtata                          >  2:147950/1‑149 (MQ=255)
                                                                                                                               gTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc                                                                                         >  3:603511/1‑80 (MQ=255)
                                                                                                                               gTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc                                                                                         <  4:603511/80‑1 (MQ=255)
                                                                                                                                  ccgatatcatcATGAAGCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGcc                 >  2:52557/6‑149 (MQ=255)
                                                                                                                                    gatatcatcATGAAGCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTg               >  1:863300/4‑149 (MQ=255)
                                                                                                                                     aaaTCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGc              >  3:21101/1‑149 (MQ=255)
                                                                                                                                          atcatGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAAt         >  2:503558/1‑149 (MQ=255)
                                                                                                                                            catGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattatt                                                              >  4:198035/1‑94 (MQ=255)
                                                                                                                                            catGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattatt                                                              <  3:198035/94‑1 (MQ=255)
                                                                                                                                            catGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGc       >  4:656210/1‑149 (MQ=255)
                                                                                                                                             atGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCt      <  2:430762/149‑1 (MQ=255)
                                                                                                                                              tGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTg     <  1:846560/149‑1 (MQ=255)
                                                                                                                                              tGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTg     <  1:925072/149‑1 (MQ=255)
                                                                                                                                              tGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTg     <  1:574870/149‑1 (MQ=255)
                                                                                                                                               gCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGTAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTgg    <  2:530289/149‑1 (MQ=255)
                                                                                                                                               gCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTgg    <  4:75637/149‑1 (MQ=255)
                                                                                                                                                 ttCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCa                     >  4:19145/1‑130 (MQ=255)
                                                                                                                                                 ttCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCa                     <  3:19145/130‑1 (MQ=255)
                                                                                                                                                 ttCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCATTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTGGCa  <  1:918446/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTGGCA  >  NC_000913/2102961‑2103254

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: