Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,642,778 G→A 5.5% R23C (CGC→TGC)  rodZ ← MreB assembly cytoskeletal protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,642,7780GA 5.5% 244.0 / 2.9 73R23C (CGC→TGC) rodZMreB assembly cytoskeletal protein
Reads supporting (aligned to +/- strand):  ref base G (30/39);  new base A (2/2);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

CCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCG  >  NC_000913/2642631‑2642918
                                                                                                                                                   |                                                                                                                                            
ccAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTcgcg                                                                                                                                              <  1:42849/148‑1 (MQ=255)
   gACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGc                                                                                                                                           <  8:16489/148‑1 (MQ=255)
   gACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGc                                                                                                                                           <  6:122225/148‑1 (MQ=255)
    aCGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCa                                                                                                                                          >  3:219287/1‑148 (MQ=255)
     cgcgCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCAt                                                                                                                                         >  7:298383/1‑148 (MQ=255)
       cgcATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTTCCTTATCTTCTTCATTGTCGCGTACCGTGGAAACCTTCAGGCATAGTCGCTCGGGAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTa                                                                                                                                       <  2:38295/148‑1 (MQ=255)
        gcATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGc                                                                                                                                           <  2:43659/143‑1 (MQ=255)
        gcATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGc                                                                                                                                           >  1:43659/1‑143 (MQ=255)
        gcATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTAc                                                                                                                                      <  2:266316/148‑1 (MQ=255)
         cATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACg                                                                                                                                     >  1:184662/1‑148 (MQ=255)
            aaGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAg                                                                                                                                  >  2:219435/1‑148 (MQ=255)
              gagCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGc                                                                                                                                >  7:208290/1‑148 (MQ=255)
              gagCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGc                                                                                                                                >  8:162011/1‑148 (MQ=255)
               agCGGATATATCCGCGCAGGAATGTTGAAGCAAGAGCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCGCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCg                                                                                                                               <  5:44237/148‑1 (MQ=255)
                         tCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAg                                                                                                                             >  5:224624/1‑140 (MQ=255)
                         tCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAg                                                                                                                             <  6:224624/140‑1 (MQ=255)
                            gcgcAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGgtagta                                                                                                                  <  4:219287/148‑1 (MQ=255)
                             cgcAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTaa                                                                                                                 <  8:298383/148‑1 (MQ=255)
                             cgcAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTaa                                                                                                                 <  5:56275/148‑1 (MQ=255)
                                aGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTg                                                                                                              <  7:191847/148‑1 (MQ=255)
                                      gTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCAt                                                                                                        >  3:296755/1‑148 (MQ=255)
                                       ttGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGgtagta                                                                                                                  <  4:130399/137‑1 (MQ=255)
                                       ttGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGgtagta                                                                                                                  >  3:130399/1‑137 (MQ=255)
                                       ttGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCAAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCAtt                                                                                                       >  2:221235/1‑148 (MQ=255)
                                          aaGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCAtttt                                                                                                     <  3:219254/147‑1 (MQ=255)
                                                 aTCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTcgtgc                                                                                             >  2:326674/1‑148 (MQ=255)
                                                 aTCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTcgtgc                                                                                             <  1:219435/148‑1 (MQ=255)
                                                    ggcggGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTcgtgcgtg                                                                                          <  3:73978/148‑1 (MQ=255)
                                                    ggcggGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTcgtgcgt                                                                                           <  8:94647/147‑1 (MQ=255)
                                                         gTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTc                                                                                     <  5:199316/148‑1 (MQ=255)
                                                            ccTTAGCTGCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGt                                                                                  <  4:20409/148‑1 (MQ=255)
                                                                aTCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTTAGTCCTAGTTGTTCGCAAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTatt                                                                               >  1:53810/1‑147 (MQ=255)
                                                                aTCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTcgtt                                                                                              <  1:25022/132‑2 (MQ=255)
                                                                aTCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTcgtt                                                                                              >  2:25022/1‑131 (MQ=255)
                                                                 tcttctTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTattca                                                                             >  2:302012/1‑148 (MQ=255)
                                                                   ttcttcAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTattcatt                                                                           >  1:247913/1‑148 (MQ=255)
                                                                     cttcAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCg                                                                         >  3:209995/1‑148 (MQ=255)
                                                                         aaTGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCTGTTGTAAGTGCTTCTTTTTGGTCGTGGGTGGCTTCAGTATTCATTCGCTGc                                                                     <  8:208290/148‑1 (MQ=255)
                                                                         aaTGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTg                                                                      <  5:58286/147‑1 (MQ=255)
                                                                         aaTGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTg                                                                      <  6:277885/147‑1 (MQ=255)
                                                                         aaTGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTg                                                                      >  8:302651/1‑147 (MQ=255)
                                                                              cgcgTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCa                                                                    >  3:18472/1‑144 (MQ=255)
                                                                              cgcgTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCa                                                                    <  4:18472/144‑1 (MQ=255)
                                                                                     cGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTaa                                                         <  1:326674/148‑1 (MQ=255)
                                                                                     cGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTaa                                                         <  4:296755/148‑1 (MQ=255)
                                                                                     cGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTaa                                                         <  7:302651/148‑1 (MQ=255)
                                                                                     cGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCAAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTaa                                                         <  2:53810/148‑1 (MQ=255)
                                                                                      gTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTaaa                                                        <  1:302012/148‑1 (MQ=255)
                                                                                           aaCCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATgaga                                                   <  7:72351/148‑1 (MQ=255)
                                                                                           aaCCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATgaga                                                   <  4:209995/148‑1 (MQ=255)
                                                                                               ttCTGGCAAAGTGGCTCGGCAACGGCCTGCTGACTAAGTCCTTGTTGTTCGCGAGCATTACGCAGGCGAGGGCCGGTTGTAAGTGCTTCTTTTTGGTCGTGGGTGGCTTCAGTATTCATTCGCTGCAGGTAGGTTTAAATGAGAATTg                                               <  1:321443/148‑1 (MQ=255)
                                                                                                 cAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTga                                              <  8:314833/147‑1 (MQ=255)
                                                                                                 cAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTga                                              <  1:157888/147‑1 (MQ=255)
                                                                                                 cAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCACTCGCTGCAGGTACGTTTAAATGAGAATTg                                               <  7:162011/146‑1 (MQ=255)
                                                                                                  aGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTgag                                             <  2:247913/147‑1 (MQ=255)
                                                                                                       aaGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCg                                       >  5:219646/1‑148 (MQ=255)
                                                                                                        aGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCgg                                      >  2:72116/1‑148 (MQ=255)
                                                                                                               cGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGt                                                                                                  >  7:235311/1‑81 (MQ=255)
                                                                                                               cGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGt                                                                                                  <  8:235311/81‑1 (MQ=255)
                                                                                                                 gCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACcatc                             <  6:219646/148‑1 (MQ=255)
                                                                                                                 gCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCAAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACcatc                             <  1:221235/148‑1 (MQ=255)
                                                                                                                        ccTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCa                      >  8:54756/1‑148 (MQ=255)
                                                                                                                           gctgACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGc                                                                     >  7:25700/1‑98 (MQ=255)
                                                                                                                           gctgACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGc                                                                     <  8:25700/98‑1 (MQ=255)
                                                                                                                           gctgACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTg                                               <  8:91685/120‑1 (MQ=255)
                                                                                                                           gctgACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTg                                               >  7:91685/1‑120 (MQ=255)
                                                                                                                           gctgACTAAGTCCTAGTTGTTCGCGAGCATCACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTgag                                             <  8:120748/122‑1 (MQ=255)
                                                                                                                           gctgACTAAGTCCTAGTTGTTCGCGAGCATCACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTgag                                             >  7:120748/1‑122 (MQ=255)
                                                                                                                                cTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACAcc              >  1:255043/1‑148 (MQ=255)
                                                                                                                                   aGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACcgc            >  8:198232/1‑147 (MQ=255)
                                                                                                                                     tCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCAATCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGa          >  2:97519/1‑147 (MQ=255)
                                                                                                                                     tCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCAGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGa          >  2:45798/1‑147 (MQ=255)
                                                                                                                                        tAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACa      >  8:181496/1‑148 (MQ=255)
                                                                                                                                            tgttCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCg  >  4:271977/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
CCAGACGCGCATAAGAGCGGATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCG  >  NC_000913/2642631‑2642918

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: