Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 611,895 G→T 5.4% S200R (AGC→AGA)  fepA ← ferrienterobactin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913611,8950GT 5.4% 299.9 / 3.6 87S200R (AGC→AGA) fepAferrienterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (47/35);  new base T (3/2);  total (50/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCCCTTCGCGCCCGGCTGGTAACGTCGTGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGCCGCCCGCCGCGCCGTTGCCATAACGCGCAGCTG  >  NC_000913/611747‑612043
                                                                                                                                                    |                                                                                                                                                    
ccccTTCGCGCCCGGCTGGTAACGTCGTGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAgg                                                                                                                                                      <  1:78529/149‑1 (MQ=255)
   cTTCGCGCCCGGCTGGTAACGTCGTGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGATGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCt                                                                                                                                                    >  5:367080/1‑148 (MQ=255)
    ttCGCGCCCGGCTGGTAACGTCGTGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTa                                                                                                                                                   >  5:242961/1‑148 (MQ=255)
      cgcgCCCGGCTGGTAACGTCGTGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTGGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTaa                                                                                                                                                  <  5:311441/147‑1 (MQ=255)
       gcgcCCGGCTGGTAACGTCGTGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAg                                                                                                                                                <  7:314482/148‑1 (MQ=255)
       gcgcCCGGCTGGTAACGTCGTGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGt                                                                                                                                               <  5:50689/149‑1 (MQ=255)
         gcCGGGCTGGTAACGTCGTGGGATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATAGAAACGGAAGCTGATTTCGTCGCCGAGCGGACCGGTCAGGGGAAAGt                                                                                                                                               >  2:139437/1‑147 (MQ=255)
           ccGGCTGGTAACGTCGTGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGt                                                                                                                                           <  4:60547/149‑1 (MQ=255)
             ggCTGGTAACGTCGTGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGc                                                                                                                                         <  1:327311/149‑1 (MQ=255)
                   tAACGTCGTGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGTCTAAAGTTAGTGCGTTtgg                                                                                                                                   >  8:67263/1‑149 (MQ=255)
                   tAACGTCGTGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGTCTAAAGTTAGTGCGTTtgg                                                                                                                                   >  8:60264/1‑149 (MQ=255)
                   tAACGTCGTGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTaaa                                                                                                                                                 <  4:1004/135‑1 (MQ=255)
                   tAACGTCGTGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTaaa                                                                                                                                                 >  3:1004/1‑135 (MQ=255)
                     aCGTCGTGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTtggt                                                                                                                                  <  4:207161/148‑1 (MQ=255)
                           tGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCt                                                                                                                           >  3:196953/1‑149 (MQ=255)
                               aTACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGTCTAAAGTTAGTGCGTTTGGTGGCACCTTCct                                                                                                                       <  7:67263/149‑1 (MQ=255)
                               aTACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGTCTAAAGTTAGTGCGTTTGGTGGCACCTTCct                                                                                                                       <  7:60264/149‑1 (MQ=255)
                                   gTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCttt                                                                                                                   <  4:196953/149‑1 (MQ=255)
                                   gTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCttt                                                                                                                   <  3:358342/149‑1 (MQ=255)
                                      ccGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTAt                                                                                                                 <  7:287753/148‑1 (MQ=255)
                                                ctgACTGATTGCCCTGGTTGATATCCCACGCTTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGTTTTCGTCTCCCATCTGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACGTTCCTCTTTATGTTCTGGCGCa                                                                                                      <  4:35784/147‑1 (MQ=255)
                                                    cTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTg                                                                                                   >  4:204086/1‑148 (MQ=255)
                                                     tGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGa                                                                                                  >  5:189629/1‑148 (MQ=255)
                                                       aTGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGaaa                                                                                                >  4:63130/1‑148 (MQ=255)
                                                           cccTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGc                                                                                           >  1:404669/1‑149 (MQ=255)
                                                            ccTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATTTTCTGGCGCATTGAAATATGCg                                                                                          >  1:51050/1‑149 (MQ=255)
                                                            ccTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCg                                                                                          >  3:161530/1‑149 (MQ=255)
                                                            ccTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCACCCAGCGGACCGGTCAGGCAAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCg                                                                                          >  1:42776/1‑149 (MQ=255)
                                                               ggTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTcc                                                                                       >  8:197842/1‑149 (MQ=255)
                                                                gTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTccc                                                                                      <  3:63130/149‑1 (MQ=255)
                                                                     tatCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAg                                                                                 >  8:118770/1‑149 (MQ=255)
                                                                      atCCCACGCGTCAGCCTGGGTTTTGTCTATGTTGCCATACAAAAGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGGGTTTGGTGTCACCTTGCTCTTTATGTTCTGGCGCATTGAAatat                                                                                             <  6:254097/136‑1 (MQ=255)
                                                                      atCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGc                                                                                >  8:283022/1‑149 (MQ=255)
                                                                      atCACACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGTAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGTCTAAAGTTAGTGGGTTTGGTGGCACCTTGCTCTTTATTTTCTGGCGCATTGAAatat                                                                                             >  5:254097/1‑136 (MQ=255)
                                                                        cccACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCg                                                                              >  6:117591/1‑149 (MQ=255)
                                                                        cccACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCg                                                                              >  5:116542/1‑149 (MQ=255)
                                                                         ccccgcgTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGt                                                                             >  2:181248/4‑149 (MQ=255)
                                                                          cACGCGTCAGCCTGGGTTTTGTCGATGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTg                                                                            >  7:332291/1‑149 (MQ=255)
                                                                           aCGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGc                                                                           <  2:265927/149‑1 (MQ=255)
                                                                                tCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTc                                                                      >  8:252070/1‑149 (MQ=255)
                                                                                    ccTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTcgccg                                                                  >  2:231060/1‑149 (MQ=255)
                                                                                     cTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCCCTcgccac                                                                 >  1:383004/1‑147 (MQ=255)
                                                                                        ggTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCt                                                                <  7:118770/147‑1 (MQ=255)
                                                                                         gTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGgc                                                                                                         >  1:349396/1‑105 (MQ=255)
                                                                                         gTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGgc                                                                                                         <  2:349396/105‑1 (MQ=255)
                                                                                            ttGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCTGTGCCACTCGCCGCTGCCttt                                                          >  8:325907/1‑149 (MQ=255)
                                                                                               tCGAGGATGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCtttttt                                                       <  6:191252/149‑1 (MQ=255)
                                                                                                 gAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTgg                                                     >  2:421909/1‑149 (MQ=255)
                                                                                                 gAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTgg                                                     <  6:242961/149‑1 (MQ=255)
                                                                                                  aGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGt                                                    <  1:231060/149‑1 (MQ=255)
                                                                                                      ggCCATACAAACCGAAGCTGAATTCGGCGCCCAGCGGACCGGTCAGGCGAAAGTTAGTGCGTTTGGGGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAatg                                                <  6:92262/148‑1 (MQ=255)
                                                                                                        ccATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAatgat                                              <  7:283022/149‑1 (MQ=255)
                                                                                                         cATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGt                                                                                         >  1:378899/1‑105 (MQ=255)
                                                                                                         cATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGt                                                                                         <  2:378899/105‑1 (MQ=255)
                                                                                                         cATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAatgatg                                             >  4:106727/1‑149 (MQ=255)
                                                                                                         cATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCAACTTCCTCTATATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCCCTCGCCGCTGCCTTTTTAGGGAatgatc                                             >  8:54840/1‑148 (MQ=255)
                                                                                                          aTACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGATTGGTGGAACCTTCCTCTTTATGTTCTGGCGCATTGTAATATGCGTCCCAGCAGCCGTGCCACACGCCGCTGCATTATTTGGTAatgatga                                            >  6:321847/1‑148 (MQ=255)
                                                                                                          aTACAAACGGAAGCTGAATTCGTCGCCCAGCGGAACGGTCAGGATAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAatgatg                                             >  6:406403/1‑148 (MQ=255)
                                                                                                           tACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCTGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGtt                                           <  7:325907/149‑1 (MQ=255)
                                                                                                            aCAAACTGAAGCTGAATTCGTCGCCCAGAGGACCGGTCGGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCATCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGtt                                           <  5:321847/148‑1 (MQ=255)
                                                                                                            aCAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTCAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGtt                                           >  1:31784/1‑148 (MQ=255)
                                                                                                             cAAACGGAAGCGGAATTCGTCGCCCAGCGGACCGGGCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTaa                                         <  4:337365/149‑1 (MQ=255)
                                                                                                              aaaCGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAAc                                        <  4:15773/149‑1 (MQ=255)
                                                                                                               aaCGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAAc                                        >  8:112245/1‑148 (MQ=255)
                                                                                                                aCGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCa                                      >  7:288779/1‑149 (MQ=255)
                                                                                                                  ggAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGt                                            >  4:303231/1‑141 (MQ=255)
                                                                                                                  ggAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGt                                            <  3:303231/141‑1 (MQ=255)
                                                                                                                   gAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTGCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCAcgc                                   <  7:374011/149‑1 (MQ=255)
                                                                                                                      gCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCAcgccgc                                >  8:104111/1‑149 (MQ=255)
                                                                                                                       cTGAATTCGTCGCCCAGCGGCCCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACgccgcc                               >  6:352933/1‑149 (MQ=255)
                                                                                                                        tGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGccgccc                              >  4:133827/1‑149 (MQ=255)
                                                                                                                        tGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGccgccc                              >  3:211284/1‑149 (MQ=255)
                                                                                                                         gAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGccgcccg                             >  7:210853/1‑149 (MQ=255)
                                                                                                                              cgtcgCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGCCGCCcgccgcg                        >  6:261004/1‑149 (MQ=255)
                                                                                                                               gtcgCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGAAGTTAACCACGCCGCCcgccgcgc                       >  7:59311/1‑149 (MQ=255)
                                                                                                                                 cgCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGtt                                           <  8:418992/127‑1 (MQ=255)
                                                                                                                                 cgCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGtt                                           >  7:418992/1‑127 (MQ=255)
                                                                                                                                 cgCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCa                                      >  1:364720/1‑132 (MQ=255)
                                                                                                                                 cgCCCACCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCa                                      <  2:364720/132‑1 (MQ=255)
                                                                                                                                  gCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCAACTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGCCGCCCGCCGCGCCGt                    >  6:87295/1‑149 (MQ=255)
                                                                                                                                   cccAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGCCGCCCGCCGCGCCGtt                   <  3:155669/149‑1 (MQ=255)
                                                                                                                                    ccAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGCCGCCCGCCGCGCCGTTg                  <  7:197842/149‑1 (MQ=255)
                                                                                                                                    ccAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGCCGCCCGCCGCGCCGTTg                  >  8:279238/1‑149 (MQ=255)
                                                                                                                                     caccGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGCCACCTTCCTCATTAAGTTCTCGCGCAGTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGCCGCCCGCCGCGCCGTTGc                 <  1:139437/146‑1 (MQ=255)
                                                                                                                                     cAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGCCGCCCGCCGCGCCGTTGc                 >  8:196808/1‑149 (MQ=255)
                                                                                                                                            ccGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTaa                                                  >  6:324309/1‑109 (MQ=255)
                                                                                                                                            ccGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTaa                                                  <  5:324309/109‑1 (MQ=255)
                                                                                                                                               gTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGCCGCCCGCCGCGCCGTTGCCATAACgcgc       <  2:18612/149‑1 (MQ=255)
                                                                                                                                                    gCTAAAGTTAGTGCGTTTGGTGGCGCCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGCCGCCCGCCGCGCCGTTGCCATAACGCGCAGCTg  >  3:209738/1‑149 (MQ=255)
                                                                                                                                                    gCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGCGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTTACCCCGCCGCCCGCCGCGCCGGTGCCATAACGcgccgcag  >  1:413296/1‑147 (MQ=255)
                                                                                                                                                    gCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGCCGCCCGCCGCGCCGTTGCCATAACGCGCAGCTg  >  1:429266/1‑149 (MQ=255)
                                                                                                                                                    gCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGCCGCCCGCCGCGCCGTTGCCATAACGCGCAGCTg  <  8:288779/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCCTTCGCGCCCGGCTGGTAACGTCGTGGCATACGTTCCGGCACGCGCGGACTGATGGCCCTGGTTGATATCCCACGCGTCAGCCTGGGTTTTGTCGAGGTTGCCATACAAACGGAAGCTGAATTCGTCGCCCAGCGGACCGGTCAGGCTAAAGTTAGTGCGTTTGGTGGCACCTTCCTCTTTATGTTCTGGCGCATTGAAATATGCGTCCCAGGAGCCGTGCCACTCGCCGCTGCCTTTTTTGGTAATGATGTTAACCACGCCGCCCGCCGCGCCGTTGCCATAACGCGCAGCTG  >  NC_000913/611747‑612043

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