Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 259,334 G→C 45.5% P256R (CCA→CGA)  phoE ← outer membrane phosphoporin protein E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913259,3340GC45.5% 58.3 / 200.4 136P256R (CCA→CGA) phoEouter membrane phosphoporin protein E
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base C (31/30);  total (67/69)
Fisher's exact test for biased strand distribution p-value = 8.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

ATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCTCGTTGGTGCGATCTGAGTTGGTATA  >  NC_000913/259189‑259479
                                                                                                                                                 |                                                                                                                                                 
aTTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCg                                                                                                                                                <  4:356204/149‑1 (MQ=255)
 ttGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGt                                                                                                                                               >  6:96286/1‑149 (MQ=255)
   gACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCa                                                                                                                                             >  1:451972/1‑149 (MQ=255)
    aCCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCAt                                                                                                                                            >  1:366930/1‑149 (MQ=255)
     ccAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCAtt                                                                                                                                           >  6:152153/1‑149 (MQ=255)
     ccAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCAtt                                                                                                                                           <  2:458423/149‑1 (MQ=255)
       aGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCAtttt                                                                                                                                         >  8:382183/1‑149 (MQ=255)
            ttCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCgtg                                                                                                                                    >  7:430477/1‑149 (MQ=255)
             tcatcaCCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCgtgt                                                                                                                                   >  1:379039/1‑149 (MQ=255)
             tcatcaCCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGc                                                                                                                                                 <  8:183261/135‑1 (MQ=255)
              catcaCCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGc                                                                                                                                       <  4:325846/144‑1 (MQ=255)
                tccccGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTc                                                                                                                                <  7:15749/146‑1 (MQ=255)
                tcaCCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGttt                                                                                                                                 >  4:224272/1‑148 (MQ=255)
                  aCCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCa                                                                                                                               <  2:366930/148‑1 (MQ=255)
                      aTACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGaa                                                                                                                            >  4:98953/1‑147 (MQ=255)
                      aTACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATa                                                                                                                          <  2:389200/149‑1 (MQ=255)
                          cTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGa                                                                                                                        >  7:326168/1‑147 (MQ=255)
                          cTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGa                                                                                                                        >  3:47259/1‑147 (MQ=255)
                          cTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGa                                                                                                                        >  3:53494/1‑147 (MQ=255)
                           ttCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCTCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCa                                                                                                                               <  2:310522/139‑1 (MQ=255)
                           ttCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAg                                                                                                                     >  8:260811/1‑149 (MQ=255)
                             cAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGt                                                                                                                    <  4:94124/148‑1 (MQ=255)
                                tatCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGcc                                                                                                                <  4:176756/149‑1 (MQ=255)
                                  tCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCa                                                                                                               <  1:4834/148‑1 (MQ=255)
                                   cTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGaaa                                                                                                                      <  2:437389/140‑1 (MQ=255)
                                   cTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGaaa                                                                                                                      >  1:437389/1‑140 (MQ=255)
                                    tttCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGa                                                                                                             >  4:284303/1‑148 (MQ=255)
                                    tttCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGa                                                                                                             >  1:157952/1‑148 (MQ=255)
                                    tttCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGa                                                                                                             >  4:278234/1‑148 (MQ=255)
                                      tCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGAt                                                                                                            <  5:437570/147‑1 (MQ=255)
                                         cTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATaaa                                                                                                         >  1:326250/1‑147 (MQ=255)
                                          tttCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAatat                                                                                                      >  2:309338/1‑149 (MQ=255)
                                              gATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATattatt                                                                                                  <  2:50053/149‑1 (MQ=255)
                                              gATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATattatt                                                                                                  <  2:157952/149‑1 (MQ=255)
                                              gATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATattatt                                                                                                  >  2:286938/1‑149 (MQ=255)
                                              gATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATattatt                                                                                                  <  7:401502/149‑1 (MQ=255)
                                                tAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTgg                                                                                                >  2:167955/1‑149 (MQ=255)
                                                tAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTgg                                                                                                >  3:462876/1‑149 (MQ=255)
                                                tAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGACAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTgg                                                                                                >  3:68596/1‑149 (MQ=255)
                                                  aGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCa                                                                                              <  1:286938/149‑1 (MQ=255)
                                                    aCATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATc                                                                                            <  5:402075/149‑1 (MQ=255)
                                                     cATAACCCAGCGATGGAAGCAGACCAAAGTCAAACTGGTATTGAGCGACCTCTTCAAAGTTCTGTGTCTTATTTGCAAAGCCGCCAGTTATTCGCGTCATTTTGAGTGTTTCAGAATAGAAATTTGCCAGATAAATATAATTGGCATCg                                                                                           >  6:253197/1‑149 (MQ=255)
                                                      aTAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCg                                                                                           >  3:217350/1‑148 (MQ=255)
                                                      aTAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCg                                                                                           >  8:104942/1‑148 (MQ=255)
                                                      aTAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCg                                                                                           >  3:194913/1‑148 (MQ=255)
                                                         aCCCAGCGATGGACGCAGACCAAACTCAAACTGGCATTGAGCGACCGCTTCAAAGTTCAGTGTCTTATTGGCAAAGCCGCCAGTTATTCCCGTCATTTTGCGTGCTTCAGAATTGAAAGTTGCCAGATAAATATTATTGGCATCGTAtt                                                                                       <  4:187291/149‑1 (MQ=255)
                                                           ccAGTGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTAttt                                                                                      >  8:6139/1‑148 (MQ=255)
                                                            cAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCGAGATAAATATTATTGGCATCGTAttt                                                                                      >  5:465610/1‑147 (MQ=255)
                                                            cAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTAttt                                                                                      >  5:230036/1‑147 (MQ=255)
                                                            cAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTAttt                                                                                      >  3:4017/1‑147 (MQ=255)
                                                            cAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCa                                                                                    <  1:309338/149‑1 (MQ=255)
                                                               cGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGGAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGAc                                                                                 <  4:194913/149‑1 (MQ=255)
                                                               cGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAAtattat                                                                                                   <  8:256195/131‑1 (MQ=255)
                                                               cGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAAtattat                                                                                                   >  7:256195/1‑131 (MQ=255)
                                                                    gACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCATATAAATATTATTGGCATCGTATTTCAGACCTGtt                                                                            <  3:98953/149‑1 (MQ=255)
                                                                    gACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGtt                                                                            <  6:488647/149‑1 (MQ=255)
                                                                       gCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGcc                                                                         >  6:9267/1‑149 (MQ=255)
                                                                             cAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGa                                                                                  >  6:97030/1‑134 (MQ=255)
                                                                             cAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGa                                                                                  <  5:97030/134‑1 (MQ=255)
                                                                                aGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCt                                                                <  7:260811/149‑1 (MQ=255)
                                                                                    aaaCTGGTATTGAGCGACCGCTTCATAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTTTTCGCGTCATTTTGCGTGTTTCAGATTAGATAGTTGCCAGATAAATATTATTGTCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCa                                                             <  3:160463/148‑1 (MQ=255)
                                                                                     aaCTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAat                                                                                                        >  3:26163/1‑104 (MQ=255)
                                                                                          gTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCtt                                                        <  8:430477/147‑1 (MQ=255)
                                                                                          gTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGc                                                      <  3:359547/149‑1 (MQ=255)
                                                                                             ttGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCt                                                    <  4:217350/148‑1 (MQ=255)
                                                                                                 gCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGc                                                >  1:90595/1‑148 (MQ=255)
                                                                                                     ccGCTTCTAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACgg                                           <  3:278234/149‑1 (MQ=255)
                                                                                                     ccGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACgg                                           <  3:284303/149‑1 (MQ=255)
                                                                                                     ccGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACgg                                           <  2:379039/149‑1 (MQ=255)
                                                                                                     ccGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACgg                                           >  4:437738/1‑149 (MQ=255)
                                                                                                       gCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCt                                         >  8:175445/1‑149 (MQ=255)
                                                                                                        cTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCtt                                        >  8:144433/1‑149 (MQ=255)
                                                                                                         ttCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGc                                                      >  4:144375/1‑134 (MQ=255)
                                                                                                           cAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCttt                                       <  1:203904/147‑1 (MQ=255)
                                                                                                            aaaGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATAGGCATCGTAATTCAGACCTGTTTCCCATGCTTCTAAACGCTTTCCTGTGCCACGGATTTg                                      >  8:431835/1‑147 (MQ=255)
                                                                                                             aaGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAg                                   <  8:326168/149‑1 (MQ=255)
                                                                                                              aGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCAc                                             <  4:445528/138‑1 (MQ=255)
                                                                                                              aGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCAc                                             >  3:445528/1‑138 (MQ=255)
                                                                                                              aGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAgg                                  <  4:4017/149‑1 (MQ=255)
                                                                                                               gTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGt                                 <  7:175445/149‑1 (MQ=255)
                                                                                                                ttCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAAAATTATTGGCATCGTATTTCAGACCTTTTGCCCATGCTTCTGCACGATAGCCTGTGCCACGGCTTTGCAGGt                                 >  7:469231/1‑148 (MQ=255)
                                                                                                                  cTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTc                               <  2:450133/148‑1 (MQ=255)
                                                                                                                    gtgtCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGc                            >  2:87914/1‑149 (MQ=255)
                                                                                                                     tgtCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGTAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGc                            <  7:431835/148‑1 (MQ=255)
                                                                                                                     tgtCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCt                           >  6:384726/1‑149 (MQ=255)
                                                                                                                       tCTTATTGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCTCg                         >  4:364238/1‑149 (MQ=255)
                                                                                                                           aTTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCTCGt                        <  4:454865/146‑1 (MQ=255)
                                                                                                                            ttGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCTCGTTGGt                    <  6:230036/149‑1 (MQ=255)
                                                                                                                            ttGGCAAAGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCTCGTTGGt                    <  5:164155/149‑1 (MQ=255)
                                                                                                                               gCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCa                                                                                    <  5:11091/82‑1 (MQ=255)
                                                                                                                               gCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCa                                                                                    >  6:11091/1‑82 (MQ=255)
                                                                                                                               gCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCTCGTTGGTGc                  >  3:452914/1‑148 (MQ=255)
                                                                                                                                 aaaGCCGCCAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCTCGTTGGTGCGAt               >  1:301845/1‑149 (MQ=255)
                                                                                                                                       gccAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCTCGTTGGTGCGATCTGAGt         >  1:390360/1‑149 (MQ=255)
                                                                                                                                        ccAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCTCGTTGGTGCGATCTGAGtt        <  1:408292/149‑1 (MQ=255)
                                                                                                                                        ccAGTTATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCTCGTTGGTGCGATCTGAGtt        >  8:486488/1‑149 (MQ=255)
                                                                                                                                         cAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCTCGTTGGTGCGATCTGAGTTg       <  1:167955/149‑1 (MQ=255)
                                                                                                                                            ttATTCGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCTCGTTGGTGCGATCTGAGTTGGta    <  3:364238/149‑1 (MQ=255)
                                                                                                                                              aTTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCTCGTTGGTGCGATCTGAGTTGGtat   <  7:144433/148‑1 (MQ=255)
                                                                                                                                               ttGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCTAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCTCGTTGGTGCGATCTGAGTTGGtata  <  6:465610/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
ATTGACCAGATCTTCATCACCGATACCTTCAATATCTTTCCCTTTCGATAAGACATAACCCAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTTGCAGGTTCTGCTCGTTGGTGCGATCTGAGTTGGTATA  >  NC_000913/259189‑259479

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