Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,640,864 A→G intergenic (+93/‑78) ytjE → / → yjtD protein YtjE/putative tRNA/rRNA methyltransferase YjtD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,640,8640AG92.5% 100.3 / 0.1 40intergenic (+93/‑78)ytjE/yjtDprotein YtjE/putative tRNA/rRNA methyltransferase YjtD
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (24/13);  total (26/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-02

TCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAACTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATAT  >  NC_000913/4640639‑4641102
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
tCTAAATGCTAGCTTTCGTCACATTCTTTTAATAGTCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTAGTAAATTCGCACCGCTtaat                                                                                                                                                                                                                         >  1:135974/1‑250 (MQ=255)
         tAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAACTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGAc                                                                                                                                                                                                                 >  1:202804/1‑249 (MQ=255)
          aGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAg                                                                                                                                                                                                               <  2:238316/250‑1 (MQ=255)
                   cacaTTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCAccct                                                                                                                                                                                                     >  1:289274/1‑250 (MQ=255)
                   cacaTTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAACTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCAc                                                                                                                                                                                                                                 >  2:500147/1‑223 (MQ=255)
                   cacaTTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAACTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCAc                                                                                                                                                                                                                                 <  1:500147/223‑1 (MQ=255)
                   cacaTTATTTTAACAATCCAACAAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTAGCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAAGTCAACCGTTAGTACAGTCAGGAAATAGGTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGACCAccct                                                                                                                                                                                                     <  2:148705/251‑2 (MQ=255)
                     caTTATTTTAATAATCCAAATAGTTGCATCATACAACTAATAAACGTGGTGAATACAATTGTCGAGATTTATTTTTTATAAAATTATACTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAATTAGTTTAGCCTTTTTTAAACTAGGTAAAGGGCTTTTTCTGCGACTTACGTTAGGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCAcactgt                                                                                                                                                                                                   >  2:456719/1‑246 (MQ=255)
                           tttAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTcgcg                                                                                                                                                                                              <  2:408834/250‑1 (MQ=255)
                              aataatCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGtt                                                                                                                                                                                           >  1:70940/1‑250 (MQ=255)
                                            ttGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTACTAAGTTGACAATGATCACCCTGTTCGCGGTTATTTGGTCaaga                                                                                                                                                                               >  1:355040/1‑248 (MQ=255)
                                             tGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCacgaag                                                                                                                                                                             >  2:31352/1‑249 (MQ=255)
                                                           aataaACGTGGTTAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGTTCAAGAAGAGTGGCAATATGCGta                                                                                                                                                             <  2:254411/251‑1 (MQ=255)
                                                                    ggtgAATCCAATTGGCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGTATATGTAGCATTTTTTAACAATACAACCGTTAGTACATTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAACTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGatt                                                                                                                                                     <  2:426707/250‑1 (MQ=255)
                                                                                    gagaTTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGTAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTATGGTCGCACCCGc                                                                                                                                    >  2:110734/1‑251 (MQ=255)
                                                                                                    aaaaTTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATAtt                                                                                                                    >  1:303334/1‑251 (MQ=255)
                                                                                                                aGTAAACAGAATAATATGTAGTTTTTTTTAACAATTCAACCGTTAGTACAGTCAGGACATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCCTTTTATGCCACTTCCGTTAAGGATTTGTAAACTCGCACCGCCTAATAAGTTGACAGTGCCCTCCCTCTTTGCGTTTATTTGGTCAACAAGAGTGGCAATATGCTTAGAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGTCAgcggcg                                                                                                        <  2:290438/251‑1 (MQ=255)
                                                                                                                  gccgatctaGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGAAAAGGGCTTTTTCTGCGACTTACGTTAAGGAGTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCTCCCTGTTCGCGGTTATTTGGTCAAGAAGAGGGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGgcgc                                                                                                       <  2:477048/242‑1 (MQ=255)
                                                                                                                    aaCAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCggg                                                                                                    >  1:141864/1‑251 (MQ=255)
                                                                                                                         aaGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCCCCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATg                                                                                               >  1:470238/1‑250 (MQ=255)
                                                                                                                          aGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGCAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGgcgc                                                                                                       >  1:477048/1‑242 (MQ=255)
                                                                                                                               atGTAGAATTTTTTAACAATTCAACCGTTAGTACAGTCAGGCAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAACACg                                                                                         >  2:321744/1‑251 (MQ=255)
                                                                                                                                   ggCATTTTTTAACACTTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGCAGAGTGGCAATATGCGTATAATGTTTATTCTGGTCGCTCCCGCCAGAGCAGAACATATTGGGGCGGCGGCGCGGGCAATGAAAACGTTgg                                                                                     >  2:3880/2‑251 (MQ=255)
                                                                                                                                           ttAACAATTCAACCGTTAGTACAGTCAGGCAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGc                                                                             >  1:162694/1‑251 (MQ=255)
                                                                                                                                            tAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCg                                                                            >  1:94522/1‑251 (MQ=255)
                                                                                                                                             aacaatTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGTCTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCg                                                                            <  2:165293/250‑1 (MQ=255)
                                                                                                                                                aattCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGaact                                                                        >  1:238430/4‑248 (MQ=255)
                                                                                                                                                      aCCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTCCGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGAt                                                                  <  1:323122/251‑1 (MQ=255)
                                                                                                                                                       ccGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGAtt                                                                 >  1:491701/1‑251 (MQ=255)
                                                                                                                                                           tAGTACAGTCAGGAAATAGTTTAGCCTTTTTTACGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTCAATTCGCACCGCGTCATAAGTTGACAGTGATCAc                                                                                                                                                                                                        >  1:110809/1‑112 (MQ=255)
                                                                                                                                                           tAGTACAGTCAGGAAATAGTTTAGCCGTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGAGCAc                                                                                                                                                                                                        <  2:110809/112‑1 (MQ=255)
                                                                                                                                                               aCAGTCATGAAATAGTTTAGCCTTTTTTAAGCTAAGGAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAATTTGACAGTGATCACCCTGTTCGCGTTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGTGGATTGTCGATAg                                                         <  2:88495/251‑1 (MQ=255)
                                                                                                                                                                   tCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGCTTTGTAAATTCGCACCGCGTACTAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCAAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGta                                                       >  2:168249/1‑248 (MQ=255)
                                                                                                                                                                   tCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGCTAGTc                                                       >  1:502048/1‑249 (MQ=255)
                                                                                                                                                                            aGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGAGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGAGATCACCCTGTTAGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGACAGCGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGTGTTTAAGCGAACTGCGGGTTGACGCAAGACAGGCCCACCTg                                            >  1:502206/1‑251 (MQ=255)
                                                                                                                                                                                  gCCTTTTTTAAGCTAAGTAAAGGGCTTTTGCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGagcca                                      <  1:110734/251‑1 (MQ=255)
                                                                                                                                                                                                  gTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCg                     >  1:87157/1‑251 (MQ=255)
                                                                                                                                                                                                      aGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCAAACCTGGAGCCAGCCACCCGCTGGGTCGcac                  >  2:447569/1‑250 (MQ=255)
                                                                                                                                                                                                              atctacGACTTAAGTTAAGGATTTGTAAATACGAACCGCGTAATAAGATGCCAGTGATCACCCTGTTCGCGGTAATTTGGTCACGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCTGAAAATATTGGGGCAGCGGCGCGGGCACTGAAAACGATGGGTTTTAGAGAACTGAGGTTTGTCGATAGTAAGGCCCACCTGGAGCCAGCACCACGCTGGGCTGCACATGGATCTg         >  2:360832/6‑251 (MQ=255)
                                                                                                                                                                                                                     aCTTACGTTAAAGATTGTTAAATTCGAACAGGGTATTAAGGTGATAGTGATCACCCTGTTCGCGGTTATTTGGACAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCTCACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGatat  <  2:409646/251‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
TCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAACTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATAT  >  NC_000913/4640639‑4641102

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