Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,127,800 T→A 5.9% intergenic (+71/‑39) yceM → / → murJ putative oxidoreductase/putative lipid II flippase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,127,8000TA 5.9% 292.6 / 7.2 84intergenic (+71/‑39)yceM/murJputative oxidoreductase/putative lipid II flippase
Reads supporting (aligned to +/- strand):  ref base T (41/38);  new base A (2/3);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

CAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTGGCGCAGGGATG  >  NC_000913/1127655‑1127946
                                                                                                                                                 |                                                                                                                                                  
cAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTg                                                                                                                                                  >  1:747640/1‑148 (MQ=255)
  aaCCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCa                                                                                                                                                >  2:448930/1‑148 (MQ=255)
  aaCCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAg                                                                                                                                               <  3:355305/149‑1 (MQ=255)
  aaCCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAg                                                                                                                                               >  4:608708/1‑149 (MQ=255)
     cgccgGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGc                                                                                                                                             <  2:396059/148‑1 (MQ=255)
      gccgGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCAAGGCAGGCgc                                                                                                                                           >  3:128635/1‑148 (MQ=255)
      gccgGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCAAGGCAGGCgc                                                                                                                                           <  4:128635/149‑2 (MQ=255)
      gccgGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCAAGGCAGGCg                                                                                                                                            <  2:137077/148‑1 (MQ=255)
          gCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTg                                                                                                                                       >  4:641412/1‑149 (MQ=255)
            gAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGc                                                                                                                                      <  2:522854/148‑1 (MQ=255)
             aaCAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGcc                                                                                                                                     >  1:358566/1‑148 (MQ=255)
              aCAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGccc                                                                                                                                    >  1:151779/1‑148 (MQ=255)
               cAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCg                                                                                                                                   >  2:891504/1‑148 (MQ=255)
                aGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTg                                                                                                                                 <  3:608708/149‑1 (MQ=255)
                 ggCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTg                                                                                                                                 >  1:616354/1‑148 (MQ=255)
                  gCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTggg                                                                                                                               >  4:108417/1‑149 (MQ=255)
                    cGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGt                                                                                                                              <  2:925022/148‑1 (MQ=255)
                      tGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCt                                                                                                                            >  1:28624/1‑148 (MQ=255)
                        cTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTgg                                                                                                                          <  1:614888/148‑1 (MQ=255)
                        cTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTgg                                                                                                                          <  1:448930/148‑1 (MQ=255)
                        cTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGa                                                                                                                         <  3:454576/149‑1 (MQ=255)
                         tGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGa                                                                                                                         >  2:18751/1‑148 (MQ=255)
                         tGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGAt                                                                                                                        <  4:612841/149‑1 (MQ=255)
                           gcgcAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGAtt                                                                                                                       <  2:726869/148‑1 (MQ=255)
                           gcgcAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCAAGGCAGGCgctga                                                                                                                                        <  3:882350/131‑5 (MQ=255)
                           gcgcAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCAAGGCAGGCgctga                                                                                                                                        >  4:882350/1‑127 (MQ=255)
                                aCGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAac                                                                                                                  >  2:721642/1‑148 (MQ=255)
                                 cGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAaca                                                                                                                 >  2:928682/1‑148 (MQ=255)
                                   tATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCg                                                                                                              <  3:160036/149‑1 (MQ=255)
                                       gTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGa                                                                                                          <  3:290437/149‑1 (MQ=255)
                                        ttGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGa                                                                                                          >  1:711490/1‑148 (MQ=255)
                                        ttGACAAGATCTGGCGCGATGCGATGAGGGAATAAACCGTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGGTTTGCCCGTGGGTCTGGATTAGAACACCGATGaa                                                                                                         >  4:435426/1‑149 (MQ=255)
                                         tGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGaa                                                                                                         <  1:589595/148‑1 (MQ=255)
                                         tGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAAt                                                                                                        <  4:462343/149‑1 (MQ=255)
                                           aCAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTa                                                                                                                      >  1:130014/1‑133 (MQ=255)
                                           aCAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTa                                                                                                                      <  2:130014/133‑1 (MQ=255)
                                            cAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAAttt                                                                                                      >  1:573811/1‑148 (MQ=255)
                                            cAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAAttt                                                                                                      >  2:820/1‑148 (MQ=255)
                                              aGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATttat                                                                                                    >  2:644792/1‑148 (MQ=255)
                                                aTCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTaa                                                                                                 <  4:375853/149‑1 (MQ=255)
                                                      cgcgATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCg                                                                                            <  2:196211/148‑1 (MQ=255)
                                                       gcgatgcgatGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCt                                                                                          <  4:748245/149‑1 (MQ=255)
                                                             cgatgAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGc                                                                                           <  1:9218/142‑1 (MQ=255)
                                                             cgatgAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGc                                                                                           >  2:9218/1‑142 (MQ=255)
                                                              gatgAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGgccgcc                                                                                   <  4:598239/149‑1 (MQ=255)
                                                                 gAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTc                                                                                >  4:403960/1‑149 (MQ=255)
                                                                  aGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTc                                                                                >  1:385931/1‑148 (MQ=255)
                                                                      aaTAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGAt                                                                                                                        <  3:821552/104‑1 (MQ=255)
                                                                      aaTAAACCCAGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGAt                                                                                                                        >  4:821552/1‑104 (MQ=255)
                                                                        tAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCg                                                                          <  1:37583/148‑1 (MQ=255)
                                                                           aCCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATg                                                                                                           <  1:312705/112‑1 (MQ=255)
                                                                           aCCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATg                                                                                                           >  2:312705/1‑112 (MQ=255)
                                                                             ccTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTg                                                                                                                                 <  2:611496/88‑1 (MQ=255)
                                                                             ccTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTg                                                                                                                                 >  1:611496/1‑88 (MQ=255)
                                                                               tGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCa                                                                   <  1:928682/148‑1 (MQ=255)
                                                                                  aaCATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCa                                                                               >  4:203585/1‑133 (MQ=255)
                                                                                  aaCATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCa                                                                               <  3:203585/133‑1 (MQ=255)
                                                                                           cGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCgtgt                                                      <  3:794704/149‑1 (MQ=255)
                                                                                           cGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCgtg                                                       <  1:320174/148‑1 (MQ=255)
                                                                                              tAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGc                                                    <  2:577891/148‑1 (MQ=255)
                                                                                                  aaTTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGc                                               >  3:667372/1‑149 (MQ=255)
                                                                                                   aTTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCAGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTCGc                                               >  2:326202/1‑148 (MQ=255)
                                                                                                    ttCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCt                                              >  1:385255/1‑148 (MQ=255)
                                                                                                     tCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCtt                                             <  2:492022/148‑1 (MQ=255)
                                                                                                        ccATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGc                                          >  2:615072/1‑148 (MQ=255)
                                                                                                        ccATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGc                                          >  2:688935/1‑148 (MQ=255)
                                                                                                         cATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCa                                         >  2:953790/1‑148 (MQ=255)
                                                                                                         cATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCAc                                        >  3:911922/1‑149 (MQ=255)
                                                                                                         cATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCAc                                        >  4:407987/1‑149 (MQ=255)
                                                                                                          aTAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACg                                       >  4:789686/1‑149 (MQ=255)
                                                                                                             aTCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACgag                                     <  1:644792/148‑1 (MQ=255)
                                                                                                              tCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACgaga                                    >  2:410629/1‑148 (MQ=255)
                                                                                                                  ggTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCa                                >  1:113023/1‑148 (MQ=255)
                                                                                                                    tAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAAtt                             >  3:39803/1‑149 (MQ=255)
                                                                                                                      gACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGt                           >  3:639128/1‑149 (MQ=255)
                                                                                                                       aCTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTc                          <  3:789686/149‑1 (MQ=255)
                                                                                                                        cTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTc                          <  1:721642/148‑1 (MQ=255)
                                                                                                                                ctctTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATc                 <  3:61509/149‑1 (MQ=255)
                                                                                                                                  ctTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCt                <  2:823387/148‑1 (MQ=255)
                                                                                                                                   tttCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCttt              >  4:370169/1‑149 (MQ=255)
                                                                                                                                   tttCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCttt              <  3:888999/149‑1 (MQ=255)
                                                                                                                                    ttCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTg             <  3:76101/149‑1 (MQ=255)
                                                                                                                                    ttCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTg             <  4:639128/149‑1 (MQ=255)
                                                                                                                                    ttCAGCGCCTGCCTTGCAGGCGTTTTGCCAGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATAACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCttt              >  1:39131/1‑148 (MQ=255)
                                                                                                                                          gcctgcctTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTGGCGCa        <  1:532352/148‑1 (MQ=255)
                                                                                                                                          gcctgcctTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTGGCGCAg       >  4:186939/1‑149 (MQ=255)
                                                                                                                                           cctgcctTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTGGCGCAg       <  1:88198/148‑1 (MQ=255)
                                                                                                                                           cctgcctTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTGGCGCAg       <  1:640380/148‑1 (MQ=255)
                                                                                                                                             tgcctTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCt                                                                            <  3:154201/77‑1 (MQ=255)
                                                                                                                                             tgcctTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCt                                                                            >  3:869279/1‑77 (MQ=255)
                                                                                                                                             tgcctTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCt                                                                            >  4:154201/1‑77 (MQ=255)
                                                                                                                                             tgcctTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCt                                                                            <  4:869279/77‑1 (MQ=255)
                                                                                                                                              gcctTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTGGCGCAGGGAt   <  3:641412/149‑1 (MQ=255)
                                                                                                                                               cctTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTc                          <  2:163224/125‑1 (MQ=255)
                                                                                                                                               cctTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTc                          >  1:163224/1‑125 (MQ=255)
                                                                                                                                               cctTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTGGCGCAGGGATg  <  4:24024/149‑1 (MQ=255)
                                                                                                                                                ctTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTGGCGCAGGGATg  <  2:593046/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTGGCGCAGGGATG  >  NC_000913/1127655‑1127946

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