Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,418,628 G→A 8.0% intergenic (+17/‑414) yhdV → / → yhdW putative outer membrane protein/pseudogene, ABC transporter periplasmic binding protein family;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,418,6280GA 8.0% 230.1 / 9.1 75intergenic (+17/‑414)yhdV/yhdWputative outer membrane protein/pseudogene, ABC transporter periplasmic binding protein family;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (33/36);  new base A (3/3);  total (36/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 = 6.35e-01

GTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTT  >  NC_000913/3418489‑3418773
                                                                                                                                           |                                                                                                                                                 
gTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGct                                                                                                                                          >  3:382589/1‑149 (MQ=255)
 ttCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGctc                                                                                                                                         <  4:669623/149‑1 (MQ=255)
       ccAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTa                                                                                                                                   >  1:421205/1‑149 (MQ=255)
       ccAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTa                                                                                                                                   >  3:584488/1‑149 (MQ=255)
        cAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTaa                                                                                                                                  >  2:508233/1‑149 (MQ=255)
        cAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTaa                                                                                                                                  >  1:79354/1‑149 (MQ=255)
               tGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACa                                                                                                                           >  3:265240/1‑149 (MQ=255)
               tGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACCGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACa                                                                                                                           >  2:69222/1‑149 (MQ=255)
                 tATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGAGCGCGGAACTCCgcgccctg                                                                                                                                       >  1:423653/1‑131 (MQ=255)
                 tATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGAGCGCGGAACTCCGCGccctg                                                                                                                                       <  2:423653/135‑2 (MQ=255)
                         gCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTcc                                                                                                                 <  4:200661/149‑1 (MQ=255)
                            cGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATa                                                                                                              <  2:382360/149‑1 (MQ=255)
                             gACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATaa                                                                                                             <  3:553795/149‑1 (MQ=255)
                             gACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATaa                                                                                                             <  1:508233/149‑1 (MQ=255)
                             gACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATaa                                                                                                             <  4:642268/149‑1 (MQ=255)
                               cTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCttttt                                                                                                                                    >  1:205146/1‑124 (MQ=255)
                               cTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCttttt                                                                                                                                    <  2:205146/124‑1 (MQ=255)
                               cTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGa                                                                                                           <  4:584488/149‑1 (MQ=255)
                                tCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATa                                                                                                                             <  4:422033/130‑1 (MQ=255)
                                tCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATa                                                                                                                             >  3:422033/1‑130 (MQ=255)
                                  aGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaa                                                                                                        >  2:170176/1‑149 (MQ=255)
                                    gATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaaaa                                                                                                      <  1:417997/149‑1 (MQ=255)
                                            ggAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCta                                                                                              >  4:205733/1‑149 (MQ=255)
                                              aaaCGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCtata                                                                                            <  3:205733/149‑1 (MQ=255)
                                              aaaCGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCtata                                                                                            >  4:491441/1‑149 (MQ=255)
                                                    cTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTctgc                                                                                      <  4:265240/149‑1 (MQ=255)
                                                     tGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTctgct                                                                                     <  4:382589/149‑1 (MQ=255)
                                                               gCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTcc                                                                           <  1:170176/149‑1 (MQ=255)
                                                                         cAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGAttt                                                                 <  2:243933/149‑1 (MQ=255)
                                                                             tgctgTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCa                                                             >  2:628378/1‑149 (MQ=255)
                                                                                tgTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGt                                                          >  3:576100/1‑149 (MQ=255)
                                                                                 gTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGtt                                                         <  3:472070/149‑1 (MQ=255)
                                                                                    gCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACa                                                      <  1:194383/149‑1 (MQ=255)
                                                                                          cTGACGATGGGGATCATCTCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATAtt                                                                                                  <  4:411493/99‑1 (MQ=255)
                                                                                          cTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATAtt                                                                                                  >  3:411493/1‑99 (MQ=255)
                                                                                            gACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATAtt                                                                                                  <  3:366978/97‑1 (MQ=255)
                                                                                            gACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATAtt                                                                                                  >  4:366978/1‑97 (MQ=255)
                                                                                            gACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTa                                              <  2:421205/149‑1 (MQ=255)
                                                                                               gATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTAttt                                           <  2:324234/149‑1 (MQ=255)
                                                                                                aTGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTc                                          >  2:339150/1‑149 (MQ=255)
                                                                                                aTGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTc                                          >  4:545631/1‑149 (MQ=255)
                                                                                                aTGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTc                                          <  3:491441/149‑1 (MQ=255)
                                                                                                  ggggATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCgc                                        <  1:69222/149‑1 (MQ=255)
                                                                                                   gggATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCgct                                       >  3:663328/1‑149 (MQ=255)
                                                                                                     gATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCgctgc                                     >  2:106870/1‑149 (MQ=255)
                                                                                                     gATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCgctgc                                     >  3:412721/1‑149 (MQ=255)
                                                                                                     gATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCgctgc                                     >  3:82935/1‑149 (MQ=255)
                                                                                                                tATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATaa                          <  4:576100/149‑1 (MQ=255)
                                                                                                                 aTCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAAc                         >  4:621862/1‑149 (MQ=255)
                                                                                                                  tCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACg                        >  1:619021/1‑149 (MQ=255)
                                                                                                                    gAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGaaa                             <  4:30087/142‑1 (MQ=255)
                                                                                                                    gAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGaaa                             >  3:30087/1‑142 (MQ=255)
                                                                                                                     aaaTAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCtat                                                                                             <  3:456875/77‑1 (MQ=255)
                                                                                                                     aaaTAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCtat                                                                                             >  4:456875/1‑77 (MQ=255)
                                                                                                                     aaaTAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATaa                          <  2:67846/144‑1 (MQ=255)
                                                                                                                     aaaTAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATaa                          >  1:67846/1‑144 (MQ=255)
                                                                                                                        tAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaa                  <  1:339150/149‑1 (MQ=255)
                                                                                                                          aCTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGAt                                                                   <  3:470546/98‑1 (MQ=255)
                                                                                                                          aCTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGAt                                                                   >  4:470546/1‑98 (MQ=255)
                                                                                                                            tGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGa                               <  3:374001/132‑1 (MQ=255)
                                                                                                                            tGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGa                               >  4:374001/1‑132 (MQ=255)
                                                                                                                              aaCAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaatttaat            <  3:621862/149‑1 (MQ=255)
                                                                                                                              aaCAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaatttaat            >  4:520238/1‑149 (MQ=255)
                                                                                                                               aCAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTctg                                                                                       >  2:654531/1‑73 (MQ=255)
                                                                                                                               aCAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTctg                                                                                       <  1:654531/73‑1 (MQ=255)
                                                                                                                               aCAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCaatttaatt           <  3:545631/149‑1 (MQ=255)
                                                                                                                                aagagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCgct                                       <  3:541358/116‑1 (MQ=255)
                                                                                                                                aagagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCgct                                       >  4:541358/5‑120 (MQ=255)
                                                                                                                                 aGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAAtttt         >  3:79303/1‑149 (MQ=255)
                                                                                                                                  gagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTc        >  4:314785/3‑149 (MQ=255)
                                                                                                                                   ggCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTca       <  4:79303/149‑1 (MQ=255)
                                                                                                                                   agcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTca       <  3:314785/148‑1 (MQ=255)
                                                                                                                                       cgGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACt   <  4:412721/149‑1 (MQ=255)
                                                                                                                                       cgGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACt   <  4:82935/149‑1 (MQ=255)
                                                                                                                                        ggAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACtt  >  2:473538/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
GTTCACACCAATACCTGTATGAAAGGCACGACTCAGGATAACTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTT  >  NC_000913/3418489‑3418773

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