Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,161,358 G→C 5.1% Q438E (CAG→GAG)  ftsP ← septal ring component that protects the divisome from stress; multicopy suppressor of ftsI(Ts)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,161,3580GC 5.1% 267.6 / 4.6 78Q438E (CAG→GAG) ftsPseptal ring component that protects the divisome from stress; multicopy suppressor of ftsI(Ts)
Reads supporting (aligned to +/- strand):  ref base G (37/37);  new base C (2/2);  total (39/39)
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.51e-01

ATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCG  >  NC_000913/3161212‑3161505
                                                                                                                                                  |                                                                                                                                                   
aTCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGAc                                                                                                                                                   <  8:323850/149‑1 (MQ=255)
 tCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCACTCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGAcc                                                                                                                                                  <  6:19591/149‑1 (MQ=255)
    aaaGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGGCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGaa                                                                                                                                               <  6:369327/149‑1 (MQ=255)
      aGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAAt                                                                                                                                             <  6:142887/149‑1 (MQ=255)
        gTAATCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATaa                                                                                                                                           <  7:89105/149‑1 (MQ=255)
           ttCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCTCGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTCACCGAAATAAAc                                                                                                                                         >  3:272688/1‑148 (MQ=255)
           ttCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATGTCCAGCGTCGGACTGTTGAAGTAGCACGGGAAGTGCGCCCCGGAAGGCTCACCGAAATAAAc                                                                                                                                         >  3:267620/1‑148 (MQ=255)
            tCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGTCCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAAc                                                                                                                                         >  7:148118/1‑147 (MQ=255)
              cATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGTATTGAACAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAagc                                                                                                                                     <  8:380802/149‑1 (MQ=255)
              cATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGTCCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAagc                                                                                                                                     <  8:148118/149‑1 (MQ=255)
               aTACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCTTATTGACCAACAGTTGCCCAATCGAGACACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCATGAAGGCTGACCGAAATAAACAagca                                                                                                                                    <  1:207586/149‑1 (MQ=255)
                tACTCGCCCTCCGCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGCGCCCCGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAagcag                                                                                                                                   >  3:228979/1‑149 (MQ=255)
                   tCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTc                                                                                                                                >  7:53033/1‑149 (MQ=255)
                       cctcctCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACt                                                                                                                            >  5:326399/1‑149 (MQ=255)
                        ctcctcGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACtt                                                                                                                           <  2:79463/149‑1 (MQ=255)
                         tcctcGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTg                                                                                                                          >  7:77803/1‑149 (MQ=255)
                           ctcGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCATTCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGt                                                                                                                         <  2:166040/148‑1 (MQ=255)
                                 cGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCaa                                                                                                                  <  8:137989/149‑1 (MQ=255)
                                   aGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACaa                                                                                                                                       >  3:45049/1‑126 (MQ=255)
                                   aGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACaa                                                                                                                                       <  4:45049/126‑1 (MQ=255)
                                   aGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAAcc                                                                                                                >  7:298936/1‑149 (MQ=255)
                                        gAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAAc                                                                                                           >  1:228785/1‑149 (MQ=255)
                                            aTTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTa                                                                                                       >  6:130097/1‑149 (MQ=255)
                                            aTTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTa                                                                                                       <  6:326399/149‑1 (MQ=255)
                                             ttACGGTACCGGATTGACCAAGAGTTGCCCAATCGAGCCACGTGACGCCATTTCCATCGTCTGACTGTTGAAGTAGAACTGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTAt                                                                                                      <  5:218958/149‑1 (MQ=255)
                                              tACTGTACCGGATTGACCAACAGTTGCGCAATCGAGCCACGGTCAGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCTCCCAGGAAGGCTGACCGATATACACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTag                                                                                                      <  6:31385/148‑2 (MQ=255)
                                                cGGTACCGGATGGACCACCAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCGGACCGAAATAAACCAGCAGCTCCACTTGTCCGGCAACCCAAACGGTATCtt                                                                                                   <  4:228979/149‑1 (MQ=255)
                                                 ggTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCttt                                                                                                  <  8:228531/149‑1 (MQ=255)
                                                        gATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCct                                                                                           >  2:314442/1‑149 (MQ=255)
                                                            gACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTg                                                                                        <  8:53033/148‑1 (MQ=255)
                                                             aCCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTc                                                                                      <  1:372515/149‑1 (MQ=255)
                                                              ccAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCt                                                                                     <  8:77803/149‑1 (MQ=255)
                                                               cAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCACCCCAAACGGTATCTTTCCAGCCTCTGTCtt                                                                                    >  6:134361/1‑149 (MQ=255)
                                                                aaCAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTc                                                                                   >  8:276411/1‑149 (MQ=255)
                                                                  cAGTTGCCCAATCGATCGACGGTCCTCCATTTCAAGCGTATGACTGTTGAATTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCACCCCAAACGGTATCTTTCCAGCCTCTGTCTTCCg                                                                                 <  8:337232/149‑1 (MQ=255)
                                                                     ttGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCggg                                                                               >  2:152534/1‑148 (MQ=255)
                                                                        cccAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACg                                                                           <  1:152534/149‑1 (MQ=255)
                                                                         ccAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAACTAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACgg                                                                          >  1:190770/1‑149 (MQ=255)
                                                                           aaTCGCGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAACCGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCa                                                                        >  4:83042/1‑149 (MQ=255)
                                                                            aTCGAGCCACGGTCCGCGATTTCCAGCGTCTGACTGTTGAAGTAGAACGTGAAGTGCGCCCAGGAAGGCTGACGGAAATAAACTAGCAGCTCCACTTGTCCGTCAACCAAAACGGTATCTTTCCAGCCTCTGTGTTCCGGGAACg                                                                           <  1:365882/145‑1 (MQ=255)
                                                                            aTCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCAt                                                                       >  1:326227/1‑149 (MQ=255)
                                                                            aTCGAGCCACGGTCCGCCATTTCCAGCGACTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGAAAGGCTGACCGAAATAAACAAGCAGCTCCACTTTTCCGTCAACCCAAACGGTAACTTTCCAGCCTCTCTCTTCCGCGAACg                                                                           >  2:365882/1‑145 (MQ=255)
                                                                                    aCGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCATGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCAt                                                               <  5:134361/149‑1 (MQ=255)
                                                                                            ccATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTa                                                       >  1:136651/1‑149 (MQ=255)
                                                                                             cATTTCCAGCGTCTGACTGTTTAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTAc                                                      >  5:197887/1‑149 (MQ=255)
                                                                                              aTTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACg                                                     >  8:23570/1‑149 (MQ=255)
                                                                                                   cAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCt                                                <  3:83042/149‑1 (MQ=255)
                                                                                                     gCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTg                                               >  8:145808/1‑148 (MQ=255)
                                                                                                         cTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACa                                          >  4:117582/1‑149 (MQ=255)
                                                                                                         cTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACa                                          >  4:112407/1‑149 (MQ=255)
                                                                                                              tGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACg                                     >  1:36844/1‑149 (MQ=255)
                                                                                                                ttGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGaa                                            >  1:284419/1‑140 (MQ=255)
                                                                                                                  gAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAAc                                           <  1:137265/139‑1 (MQ=255)
                                                                                                                  gAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAAc                                           >  2:137265/1‑139 (MQ=255)
                                                                                                                   aaGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCtt                                 >  1:374023/1‑148 (MQ=255)
                                                                                                                     gTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCaa                              >  8:304189/1‑149 (MQ=255)
                                                                                                                        gAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAatat                           <  6:197887/149‑1 (MQ=255)
                                                                                                                        gAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAatat                           >  4:96415/1‑149 (MQ=255)
                                                                                                                         aaCGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATg                          <  1:407233/149‑1 (MQ=255)
                                                                                                                          aCGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATgg                         <  1:314442/149‑1 (MQ=255)
                                                                                                                           cGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATgg                         >  2:76954/1‑148 (MQ=255)
                                                                                                                            gggAAGTGCGCCCAGGAAGGCTCACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGaa                       <  4:272688/149‑1 (MQ=255)
                                                                                                                            gggAAGTGCGCCCAGGAAGGCTCACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGGCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGaa                       <  4:267620/149‑1 (MQ=255)
                                                                                                                             ggAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAAc                      <  8:298936/149‑1 (MQ=255)
                                                                                                                                aGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCt                   >  7:103894/1‑149 (MQ=255)
                                                                                                                                   gcgcCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGc                <  3:112407/149‑1 (MQ=255)
                                                                                                                                    cgcCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGc                <  3:117582/148‑1 (MQ=255)
                                                                                                                                    cgcCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCg               <  3:96415/149‑1 (MQ=255)
                                                                                                                                     gcCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCgg              >  1:153131/1‑149 (MQ=255)
                                                                                                                                     gcCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCgg              >  7:102332/1‑149 (MQ=255)
                                                                                                                                        cAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATgg                         >  4:177696/1‑135 (MQ=255)
                                                                                                                                        cAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATgg                         <  3:177696/135‑1 (MQ=255)
                                                                                                                                        cAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTc           <  2:326227/149‑1 (MQ=255)
                                                                                                                                           gAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCa                                                                                               >  8:45707/1‑62 (MQ=255)
                                                                                                                                           gAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCa                                                                                               <  7:45707/62‑1 (MQ=255)
                                                                                                                                               gCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGc                >  2:361983/1‑137 (MQ=255)
                                                                                                                                               gCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGc                <  1:361983/137‑1 (MQ=255)
                                                                                                                                                 tGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCcgcgcg  >  1:114378/1‑149 (MQ=255)
                                                                                                                                                  gACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGaa                                            <  8:276162/106‑1 (MQ=255)
                                                                                                                                                  gACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGaa                                            >  7:276162/1‑106 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCG  >  NC_000913/3161212‑3161505

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: