Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,118,531 A→G 5.7% E6G (GAA→GGA)  zapB → FtsZ stabilizer, septal ring assembly factor, cell division stimulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,118,5310AG 5.7% 230.6 / 2.4 72E6G (GAA→GGA) zapBFtsZ stabilizer, septal ring assembly factor, cell division stimulator
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base G (2/2);  total (36/36)
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.73e-01

GGACAGCGATGGCGCGGCTTCATACAATCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGA  >  NC_000913/4118392‑4118675
                                                                                                                                           |                                                                                                                                                
ggACAGCGATGGCGCGGCTTCATACAATCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTga                                                                                                                                         <  4:382238/149‑1 (MQ=255)
 gACAGCGATGGCGCGGCTTCATACAATCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTgag                                                                                                                                        <  3:194038/149‑1 (MQ=255)
        aTGGCGCGGCTTCATACAATCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTgg                                                                                                                                 <  1:328068/149‑1 (MQ=255)
          ggCGCGGCTTCATACAATCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGaa                                                                                                                               <  3:529688/149‑1 (MQ=255)
              cgGCTTCATACAATCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCaa                                                                                                                           <  4:167969/149‑1 (MQ=255)
               ggCTTCATACAATCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCaaa                                                                                                                          >  2:6112/1‑149 (MQ=255)
                gCTTCATACAATCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCaaaa                                                                                                                         <  2:480023/149‑1 (MQ=255)
                gCTTCATACAATCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCaaaa                                                                                                                         <  1:406995/149‑1 (MQ=255)
                  ttCATACAATCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGt                                                                                                                       >  1:299344/1‑149 (MQ=255)
                  ttCATACAATCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGt                                                                                                                       >  3:561478/1‑149 (MQ=255)
                     aTACAATCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTAca                                                                                                                    <  3:97765/149‑1 (MQ=255)
                      tACAATCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGC‑AAAGTAcagc                                                                                                                  >  2:412827/1‑149 (MQ=255)
                          aTCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAgg                                                                                                               <  1:467604/149‑1 (MQ=255)
                           tCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGc                                                                                                              <  4:78819/149‑1 (MQ=255)
                             ggAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGTAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGa                                                                                                            <  2:462392/149‑1 (MQ=255)
                               aGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCACTTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGAtt                                                                                                          >  4:98501/1‑149 (MQ=255)
                                   aactaaAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATa                                                                                                      >  4:304749/1‑149 (MQ=255)
                                    actaaAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATAc                                                                                                     >  1:284755/1‑149 (MQ=255)
                                    actaaAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATAc                                                                                                     >  2:545230/1‑149 (MQ=255)
                                         aGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCa                                                                                                <  1:545230/149‑1 (MQ=255)
                                           tGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCAct                                                                                              >  3:435612/1‑149 (MQ=255)
                                              gcTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTg                                                                                           >  3:266654/1‑149 (MQ=255)
                                                  gTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGc                                                                                       >  1:294219/1‑149 (MQ=255)
                                                   tatttattAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCa                                                                                      >  2:460652/1‑149 (MQ=255)
                                                     tttattAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGttt                                                                                                                                           >  4:221640/1‑94 (MQ=255)
                                                     tttattAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGttt                                                                                                                                           <  3:221640/94‑1 (MQ=255)
                                                           aaGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGaaa                                                                              <  3:101995/149‑1 (MQ=255)
                                                             ggCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATc                                                                            <  2:488675/149‑1 (MQ=255)
                                                                gTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGaa                                                                         <  4:23108/149‑1 (MQ=255)
                                                                gTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGaa                                                                         <  4:184125/149‑1 (MQ=255)
                                                                   aCCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAagag                                                                      >  4:474386/1‑149 (MQ=255)
                                                                       gTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTga                                                                  <  2:108448/149‑1 (MQ=255)
                                                                         aaTCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTgaaa                                                                <  3:45700/149‑1 (MQ=255)
                                                                                aCGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGaaaaaa                                                         <  4:398127/149‑1 (MQ=255)
                                                                                 cGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGaaaaaaa                                                        <  4:435612/149‑1 (MQ=255)
                                                                                  gAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAaac                                                       >  2:66111/1‑149 (MQ=255)
                                                                                  gAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAaac                                                       >  4:518746/1‑149 (MQ=255)
                                                                                   aGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAaaca                                                      >  2:145544/1‑149 (MQ=255)
                                                                                    ggATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAaacaa                                                     <  1:482368/149‑1 (MQ=255)
                                                                                    ggATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAaacaa                                                     <  1:594209/149‑1 (MQ=255)
                                                                                     gATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAaacaac                                                    >  4:232181/1‑149 (MQ=255)
                                                                                       tttttATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATg                                                                                  >  3:103030/1‑117 (MQ=255)
                                                                                       tttttATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATg                                                                                  <  4:103030/117‑1 (MQ=255)
                                                                                         tttATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCAc                                                <  2:297244/149‑1 (MQ=255)
                                                                                              ccATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGGTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCg                                           >  2:131426/1‑149 (MQ=255)
                                                                                               cATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGc                                          >  4:388663/1‑149 (MQ=255)
                                                                                                aTCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTgaaa                                                                <  4:297435/126‑1 (MQ=255)
                                                                                                aTCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCa                                         >  2:491454/1‑149 (MQ=255)
                                                                                                aTCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATCAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTgaaa                                                                >  3:297435/1‑126 (MQ=255)
                                                                                                 tCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAg                                        >  4:351311/1‑149 (MQ=255)
                                                                                                 tCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAg                                        >  4:569510/1‑149 (MQ=255)
                                                                                                 tCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAg                                        >  4:238519/1‑149 (MQ=255)
                                                                                                   aaCGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGa                                      >  2:206822/1‑149 (MQ=255)
                                                                                                    aCGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGaa                                     <  2:284755/149‑1 (MQ=255)
                                                                                                    aCGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGaa                                     <  1:145544/149‑1 (MQ=255)
                                                                                                    aCGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGaa                                     >  4:287578/1‑149 (MQ=255)
                                                                                                        cTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTc                                 >  2:459778/1‑149 (MQ=255)
                                                                                                               ttCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGc                          <  3:287578/149‑1 (MQ=255)
                                                                                                                 cAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAG‑AAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCa                       <  3:474386/149‑1 (MQ=255)
                                                                                                                  aGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCa                       >  2:332693/1‑149 (MQ=255)
                                                                                                                       agGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATc                  <  1:454236/149‑1 (MQ=255)
                                                                                                                         gTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAg                >  2:303329/1‑149 (MQ=255)
                                                                                                                          tATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAgc               >  3:299447/1‑149 (MQ=255)
                                                                                                                           aTGACAATGTCATTAGGAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAaacaa                                                     >  3:94735/1‑110 (MQ=255)
                                                                                                                           aTGACAATGTCATTAGGAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAaacaa                                                     >  1:534882/1‑110 (MQ=255)
                                                                                                                           aTGACAATGTCATTAGGAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAaacaa                                                     <  2:534882/110‑1 (MQ=255)
                                                                                                                           aTGACAATGTCATTAGGAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAaacaa                                                     <  4:94735/110‑1 (MQ=255)
                                                                                                                                 aTGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAaga        >  2:402575/1‑149 (MQ=255)
                                                                                                                                  tGTCATTAGAAGTGTTTGAGAAACTGGAAGC‑AAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGc      <  1:412827/149‑1 (MQ=255)
                                                                                                                                     cATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTg    <  4:561478/149‑1 (MQ=255)
                                                                                                                                      aTTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTgg   <  1:459778/149‑1 (MQ=255)
                                                                                                                                       ttAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGa  >  2:374225/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                
GGACAGCGATGGCGCGGCTTCATACAATCGGAGCTAACTAAAGTGCGCTCGTATTTATTAAGGCGTCACCGGTAATCGGGACGAGGATTTTTATCCATCAACGCCTTGCAATTCAGGAGAGGTATGACAATGTCATTAGAAGTGTTTGAGAAACTGGAAGCAAAAGTACAGCAGGCGATTGATACCATCACTCTGTTGCAGATGGAAATCGAAGAGCTGAAAGAAAAAAACAACTCACTGTCGCAGGAAGTTCAAAATGCCCAGCATCAGCGCGAAGAGCTGGA  >  NC_000913/4118392‑4118675

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