Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,601,824 G→C 8.1% P88R (CCC→CGC) ‡ ftsX ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,601,8240GC 8.1% 240.1 / 12.7 76P88R (CCC→CGC) ‡ftsXputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (36/33);  new base C (3/3);  total (40/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 = 2.63e-01

CACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAGGATT  >  NC_000913/3601676‑3601971
                                                                                                                                                    |                                                                                                                                                   
cACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTgg                                                                                                                                                     >  1:652413/1‑149 (MQ=255)
   gCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCa                                                                                                                                                  >  4:748631/1‑149 (MQ=255)
     cTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGc                                                                                                                                                <  4:227658/149‑1 (MQ=255)
         cTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCa                                                                                                                                            <  4:463442/149‑1 (MQ=255)
           cGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCaga                                                                                                                                          <  4:24908/149‑1 (MQ=255)
            ggCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTctgtt                                                                                                                                         <  2:445299/149‑5 (MQ=255)
                tGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAAt                                                                                                                                     <  1:625862/149‑1 (MQ=255)
                 gCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATg                                                                                                                                    <  3:235139/149‑1 (MQ=255)
                     cTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGa                                                                                                                                >  4:261273/1‑149 (MQ=255)
                        tgCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTct                                                                                                                                            <  2:639940/134‑2 (MQ=255)
                        tgCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTct                                                                                                                                            >  1:639940/1‑133 (MQ=255)
                           caccacGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAATAGTGATTTGCTGTGACGGATAATACTGCGCCGCAGCCTTGTTAACGTGTTTGTACCCCATACAACAGGCGTGGGGCAGCGTCAGAGAAATGGCGATAACCa                                                                                                                          <  2:319968/149‑1 (MQ=255)
                            accacGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCAt                                                                                                                         <  1:494261/149‑1 (MQ=255)
                             ccacGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATc                                                                                                                        >  4:433781/1‑149 (MQ=255)
                             ccacGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATc                                                                                                                        >  3:595099/1‑149 (MQ=255)
                                    cGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTa                                                                                                                 <  2:268448/149‑1 (MQ=255)
                                            cGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTGTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACAATCACCTTTAAAAAtgtg                                                                                                         >  2:273763/1‑149 (MQ=255)
                                             gTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTgg                                                                                                        <  4:705666/149‑1 (MQ=255)
                                              tCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGc                                                                                                       >  3:488856/1‑149 (MQ=255)
                                                    tCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTaaaaa                                                                                                             <  3:621845/137‑1 (MQ=255)
                                                    tCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTaaaaa                                                                                                             >  4:621845/1‑137 (MQ=255)
                                                      cAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGtt                                                                                               >  1:675880/1‑149 (MQ=255)
                                                       aGCGTTTTTTGCAGATAAACAGTGATCTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTAAACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGttt                                                                                              >  3:635577/1‑149 (MQ=255)
                                                           ttttttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCtt                                                                                          <  2:470977/149‑1 (MQ=255)
                                                              tttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTc                                                                                       <  1:42665/149‑1 (MQ=255)
                                                              tttGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTc                                                                                       <  1:630826/149‑1 (MQ=255)
                                                                   aGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATc                                                                                  >  2:318787/1‑149 (MQ=255)
                                                                   aGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATc                                                                                  >  4:661780/1‑149 (MQ=255)
                                                                    gATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATcc                                                                                 >  4:21681/1‑149 (MQ=255)
                                                                     aTAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCt                                                                                >  2:195631/1‑149 (MQ=255)
                                                                       aaaCAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGc                                                                              >  2:673914/1‑149 (MQ=255)
                                                                          cAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAAt                                                                           <  1:761287/149‑1 (MQ=255)
                                                                            gTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTaaaa                                                                                                              >  4:148516/1‑112 (MQ=255)
                                                                            gTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTaaaa                                                                                                              <  3:148516/112‑1 (MQ=255)
                                                                              gATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGGAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGCTCCTGCAATgcgg                                                                       <  4:491364/149‑2 (MQ=255)
                                                                                  tggagagACGGATAATACTTCGCCGCCGGCTGGTTAACGTTTTTCTACACCATATAGAAGACCAACGGCAGCGTACGAGCAATAGCGATAACCATCACCGTTAAAAGCGTGGAGAACGGTTTGCTTTTCAGCTCCAGCAATGCGCCGTg                                                                   <  2:56368/143‑1 (MQ=255)
                                                                                     ggTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGaa                                                                >  1:172189/1‑149 (MQ=255)
                                                                                     ggTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGaa                                                                <  4:537233/149‑1 (MQ=255)
                                                                                       tGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGACAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAgg                                                              <  3:58414/149‑1 (MQ=255)
                                                                                        gACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGc                                                             <  2:744490/149‑1 (MQ=255)
                                                                                            gATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAg                                                         >  3:439102/1‑149 (MQ=255)
                                                                                            gATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAg                                                         >  1:173924/1‑149 (MQ=255)
                                                                                             ataataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAgc                                                        >  1:669767/1‑149 (MQ=255)
                                                                                                ataCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCa                                                     <  2:675880/149‑1 (MQ=255)
                                                                                                  aCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCAcc                                                   >  3:592413/1‑149 (MQ=255)
                                                                                                     gCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGt                                                <  4:592413/149‑1 (MQ=255)
                                                                                                          gccgccTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGtt                                           <  1:180521/149‑1 (MQ=255)
                                                                                                           ccgccTGGTAAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTg                                          <  1:528272/149‑1 (MQ=255)
                                                                                                                tGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAAc                                     >  1:368325/1‑149 (MQ=255)
                                                                                                                   ttAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATATCCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTc                                              <  4:244061/137‑1 (MQ=255)
                                                                                                                   ttAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAg                                                               >  3:308975/1‑120 (MQ=255)
                                                                                                                   ttAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAg                                                               <  4:308975/120‑1 (MQ=255)
                                                                                                                   ttAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTc                                              >  3:244061/1‑137 (MQ=255)
                                                                                                                   ttAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGtt                                  <  3:410597/149‑1 (MQ=255)
                                                                                                                        gTTTTTGTACACCATATAACAGACGCTGAGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGttt                             >  4:276772/1‑149 (MQ=255)
                                                                                                                               tACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTg                                                                                                         <  4:428135/66‑1 (MQ=255)
                                                                                                                               tACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTg                                                                                                         >  3:428135/1‑66 (MQ=255)
                                                                                                                                 cacCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCa                                                                                      <  2:347700/83‑1 (MQ=255)
                                                                                                                                 cacCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCa                                                                                      >  1:347700/1‑83 (MQ=255)
                                                                                                                                 cacCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTa                    <  4:615800/149‑1 (MQ=255)
                                                                                                                                    cATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCg                 >  2:575467/1‑149 (MQ=255)
                                                                                                                                     atatAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCgg                >  4:703529/1‑149 (MQ=255)
                                                                                                                                     atatAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCgg                >  4:369050/1‑149 (MQ=255)
                                                                                                                                      tttctctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGt               >  2:470088/8‑149 (MQ=255)
                                                                                                                                      tataACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGt               >  4:649232/1‑149 (MQ=255)
                                                                                                                                         aaCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGtttt            <  3:264301/149‑1 (MQ=255)
                                                                                                                                          tctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGtt                                           <  2:497962/114‑1 (MQ=255)
                                                                                                                                          tctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTAGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGtt                                           >  1:497962/4‑117 (MQ=255)
                                                                                                                                          aCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTg           >  2:381433/1‑149 (MQ=255)
                                                                                                                                          aCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTg           >  2:618911/1‑149 (MQ=255)
                                                                                                                                           cAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGtt                                           >  1:522314/1‑116 (MQ=255)
                                                                                                                                           cAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGtt                                           <  2:522314/116‑1 (MQ=255)
                                                                                                                                            aGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGc         >  1:142513/1‑149 (MQ=255)
                                                                                                                                            aGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGc         >  3:541362/1‑149 (MQ=255)
                                                                                                                                            aGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGc         >  3:56077/1‑149 (MQ=255)
                                                                                                                                             gACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCg        >  3:86352/1‑149 (MQ=255)
                                                                                                                                               cGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCCTTTCTGACCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAg      <  2:552171/149‑1 (MQ=255)
                                                                                                                                                   gggCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAGGAtt  >  4:62269/1‑149 (MQ=255)
                                                                                                                                                   gggCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAGGAtt  >  4:101269/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAGGATT  >  NC_000913/3601676‑3601971

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