Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,015,912 T→A 21.6% intergenic (+30/+40) rmf → / ← fabA ribosome modulation factor/beta‑hydroxydecanoyl thioester dehydrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,015,9120TA21.6% 166.8 / 42.9 74intergenic (+30/+40)rmf/fabAribosome modulation factor/beta‑hydroxydecanoyl thioester dehydrase
Reads supporting (aligned to +/- strand):  ref base T (27/31);  new base A (10/6);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CAGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATC  >  NC_000913/1015766‑1016054
                                                                                                                                                  |                                                                                                                                              
cAGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAtt                                                                                                                                               >  2:89893/1‑148 (MQ=255)
cAGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTg                                                                                                                                              >  4:248255/1‑149 (MQ=255)
 aGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATACTCTCTCTTTAAAAAGAAACCTCCGCATTGc                                                                                                                                             <  4:752475/149‑1 (MQ=255)
  ggCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGc                                                                                                                                             >  1:356467/1‑148 (MQ=255)
    ccGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCgg                                                                                                                                           <  2:696537/148‑1 (MQ=255)
    ccGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCgg                                                                                                                                           <  2:249210/148‑1 (MQ=255)
    ccGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGa                                                                                                                                          <  4:829726/149‑1 (MQ=255)
          tCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTc                                                                                                                                    <  3:159836/149‑1 (MQ=255)
            gCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCg                                                                                                                                   <  1:832210/148‑1 (MQ=255)
             ccGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCgg                                                                                                                                           <  4:593847/139‑1 (MQ=255)
             ccGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCgg                                                                                                                                           >  3:593847/1‑139 (MQ=255)
                     tCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGAt                                                                                                                          >  1:250065/1‑148 (MQ=255)
                        aaaGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATAct                                                                                                                       >  2:242100/1‑148 (MQ=255)
                        aaaGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAGGtt                                                                                                                                      <  1:129149/133‑1 (MQ=255)
                        aaaGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAGGtt                                                                                                                                      >  2:129149/1‑133 (MQ=255)
                            aaaTGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATCCtctgt                                                                                                                   <  2:491721/148‑1 (MQ=255)
                            aaaTGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgt                                                                                                                   >  1:491721/1‑148 (MQ=255)
                            aaaTGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgt                                                                                                                   <  1:89893/148‑1 (MQ=255)
                               tgtgTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGa                                                                                                               <  4:872194/149‑1 (MQ=255)
                                 tgtCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGAtt                                                                                                             <  4:223406/149‑1 (MQ=255)
                                   tCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAGGt                                                                                                                                       <  4:217511/121‑1 (MQ=255)
                                   tCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAGGt                                                                                                                                       >  3:217511/1‑121 (MQ=255)
                                     ccTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATaa                                                                                                         >  3:169763/1‑149 (MQ=255)
                                        aTCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATc                                                                                                       >  1:503942/1‑148 (MQ=255)
                                           aGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGaa                                                                                                   >  3:408867/1‑149 (MQ=255)
                                           aGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGa                                                                                                    <  1:431287/148‑1 (MQ=255)
                                           aGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGa                                                                                                    >  2:382609/1‑148 (MQ=255)
                                                 tGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGa                                                                                             >  4:411927/1‑149 (MQ=255)
                                                  gAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGa                                                                                             >  1:394449/1‑148 (MQ=255)
                                                  gAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGa                                                                                             <  2:895612/148‑1 (MQ=255)
                                                    aTCAAAGGGCACAATGGCTGGGAGGCTGGCGAGAAGCCATAGCGGACAGGTTAGTAATGGCATGATTCTGTCTCTTTAAAAAGAAACCTCAGCATTGCGGAGGTCTCGCCCTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGt                                                                                          <  3:562056/149‑1 (MQ=255)
                                                      cAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGa                                                                                                                                          <  2:645586/99‑1 (MQ=255)
                                                      cAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGa                                                                                                                                          >  1:645586/1‑99 (MQ=255)
                                                        aaGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATcc                                                                                      >  3:893014/1‑149 (MQ=255)
                                                        aaGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATcc                                                                                      >  3:290116/1‑149 (MQ=255)
                                                         aGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAACCCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGaa                                                                                                   >  2:800859/1‑135 (MQ=255)
                                                         aGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGaa                                                                                                   <  1:800859/135‑1 (MQ=255)
                                                           gTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCGGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAgg                                                                                                                                        <  1:89132/96‑1 (MQ=255)
                                                           gTCACAATGGCTGGGAGCCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAgg                                                                                                                                        >  2:89132/1‑96 (MQ=255)
                                                              acaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAAc                                                                                >  3:236453/1‑149 (MQ=255)
                                                                gaTGGCTGGAAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACac                                                                              >  4:780192/2‑148 (MQ=255)
                                                                 aTGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAg                                                                              >  2:598578/1‑148 (MQ=255)
                                                                     cTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCg                                                                          >  1:657389/1‑148 (MQ=255)
                                                                      tGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGa                                                                         <  1:242100/148‑1 (MQ=255)
                                                                          aGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACttt                                                                     <  1:598578/148‑1 (MQ=255)
                                                                           ggCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTc                                                                    <  1:853387/148‑1 (MQ=255)
                                                                            gCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCa                                                                   <  1:512520/148‑1 (MQ=255)
                                                                            gCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAg                                                                  >  3:899339/1‑149 (MQ=255)
                                                                               ggCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTc                                                               <  3:411927/149‑1 (MQ=255)
                                                                                      aGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGc                                                         <  2:657389/148‑1 (MQ=255)
                                                                                      aGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCg                                                        <  3:248255/149‑1 (MQ=255)
                                                                                               ggACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGAt                                                <  2:394449/148‑1 (MQ=255)
                                                                                               ggACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGAt                                                <  2:250065/148‑1 (MQ=255)
                                                                                                gACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTGCGCATTGCAGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCa                                              <  3:780192/149‑1 (MQ=255)
                                                                                                             tGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAAc                                 <  4:408867/149‑1 (MQ=255)
                                                                                                                 cTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAg                              >  1:752371/1‑148 (MQ=255)
                                                                                                                     ttCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACtt                         >  3:474940/1‑149 (MQ=255)
                                                                                                                            tctTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTCCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATc                  >  4:880566/1‑149 (MQ=255)
                                                                                                                             ctTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATc                  >  1:450075/1‑148 (MQ=255)
                                                                                                                               ttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGc               <  4:169763/149‑1 (MQ=255)
                                                                                                                                  aaaaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAg            <  4:290116/149‑1 (MQ=255)
                                                                                                                                  aaaaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAg            <  4:893014/149‑1 (MQ=255)
                                                                                                                                     aGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAgg           <  3:622379/147‑1 (MQ=255)
                                                                                                                                     aGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAgg           >  4:622379/1‑147 (MQ=255)
                                                                                                                                       aaaCCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCa       <  4:899339/149‑1 (MQ=255)
                                                                                                                                        aaCCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCa       >  1:169490/1‑148 (MQ=255)
                                                                                                                                        aaCCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCa       >  1:830098/1‑148 (MQ=255)
                                                                                                                                         aCCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATa     >  4:727070/1‑149 (MQ=255)
                                                                                                                                           cTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAAt   >  4:95602/1‑149 (MQ=255)
                                                                                                                                            tCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGAGCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGcc                     >  4:42954/1‑130 (MQ=255)
                                                                                                                                            tCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGAGCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGcc                     <  3:42954/130‑1 (MQ=255)
                                                                                                                                            tCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGTTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATc  >  3:336058/1‑149 (MQ=255)
                                                                                                                                            tCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAAt   >  2:66433/1‑148 (MQ=255)
                                                                                                                                             ccGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATc  >  2:223709/1‑148 (MQ=255)
                                                                                                                                             ccGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATc  <  2:169490/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CAGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATC  >  NC_000913/1015766‑1016054

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: