Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,838,222 A→T 6.0% intergenic (‑117/+32) hydN ← / ← ascG formate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/asc operon transcriptional repressor; prpBC operon repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,838,2220AT 6.0% 353.8 / 9.4 100intergenic (‑117/+32)hydN/ascGformate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/asc operon transcriptional repressor; prpBC operon repressor
Reads supporting (aligned to +/- strand):  ref base A (50/44);  new base T (2/4);  total (52/48)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCT  >  NC_000913/2838081‑2838366
                                                                                                                                             |                                                                                                                                                
cgTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAg                                                                                                                                            <  1:379676/148‑1 (MQ=255)
  tCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAAAGACAGTTAATGCCCCCATACATTCACAAATTCAGCCAGATTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTAATCGACACTAGTGa                                                                                                                                         <  3:64230/149‑1 (MQ=255)
   cAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGa                                                                                                                                         >  1:676599/1‑148 (MQ=255)
     gCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACg                                                                                                                                       <  2:858224/148‑1 (MQ=255)
     gCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACg                                                                                                                                       <  2:378984/148‑1 (MQ=255)
      cAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGAt                                                                                                                                     <  4:68359/149‑1 (MQ=255)
      cAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGAt                                                                                                                                     <  4:25571/149‑1 (MQ=255)
       aaTGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATg                                                                                                                                    <  3:132744/149‑1 (MQ=255)
         tgatgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGa                                                                                                                                   <  1:502569/148‑1 (MQ=255)
          gatgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGAc                                                                                                                                  >  1:342328/1‑148 (MQ=255)
            tgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGACAtc                                                                                                                               <  4:897815/149‑2 (MQ=255)
             gaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtg                                                                                                                               >  2:561975/1‑148 (MQ=255)
             gaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtg                                                                                                                               <  1:732568/148‑1 (MQ=255)
             gaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATACCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgt                                                                                                                              <  3:897730/149‑1 (MQ=255)
              aaaCGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgtg                                                                                                                             >  3:198167/1‑149 (MQ=255)
                  ggTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtgatg                                                                                                                          <  2:130938/148‑1 (MQ=255)
                   gTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGAc                                                                                                                        >  3:678365/1‑149 (MQ=255)
                    ttCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGAc                                                                                                                        >  1:456215/1‑148 (MQ=255)
                        tcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATg                                                                                                                    >  1:244534/1‑148 (MQ=255)
                         caAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGt                                                                                                                   >  1:498329/1‑148 (MQ=255)
                          aaaTTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGttt                                                                                                                 >  4:646879/1‑149 (MQ=255)
                          aaaTTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGttt                                                                                                                 <  4:112147/149‑1 (MQ=255)
                           aaTTCTCCAGCAATGACAGTTAATGCGTCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTa                                                                                                                >  3:103338/1‑149 (MQ=255)
                            aTTCTCCAGCAATGGCAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTa                                                                                                                <  2:539598/148‑1 (MQ=255)
                            aTTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTAt                                                                                                               >  4:676134/1‑149 (MQ=255)
                             ttCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATc                                                                                                              <  3:676134/149‑1 (MQ=255)
                                  cAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGa                                                                                                          <  1:545972/148‑1 (MQ=255)
                                    gCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGAc                                                                                                                        >  4:73357/1‑132 (MQ=255)
                                    gCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGAc                                                                                                                        <  3:73357/132‑1 (MQ=255)
                                    gCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAg                                                                                                        <  2:864094/148‑1 (MQ=255)
                                        tGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCa                                                                                                   >  4:659702/1‑149 (MQ=255)
                                         gACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCaa                                                                                                  >  4:741331/1‑149 (MQ=255)
                                          aCAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGtt                                                                                                                  <  2:322279/132‑1 (MQ=255)
                                          aCAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGtt                                                                                                                  >  1:322279/1‑132 (MQ=255)
                                          aCAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCaa                                                                                                  <  2:342328/148‑1 (MQ=255)
                                          aCAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCaa                                                                                                  >  2:304279/1‑148 (MQ=255)
                                              ttAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgag                                                                                              >  1:529279/1‑148 (MQ=255)
                                              ttAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgag                                                                                              >  2:133595/1‑148 (MQ=255)
                                              ttAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATGTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgag                                                                                              >  2:803650/1‑148 (MQ=255)
                                                 aTGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagtgag                                                                                          >  3:124696/1‑149 (MQ=255)
                                                 aTGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagtgag                                                                                          >  4:468213/1‑149 (MQ=255)
                                                  tGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagtgag                                                                                          <  2:147443/148‑1 (MQ=255)
                                                       cGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCgcgg                                                                                    >  3:729943/1‑149 (MQ=255)
                                                       cGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCgcgg                                                                                    >  3:276235/1‑149 (MQ=255)
                                                       cGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCgcgg                                                                                    <  3:263406/149‑1 (MQ=255)
                                                            cATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGAt                                                                               <  4:276235/149‑1 (MQ=255)
                                                             aTTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGAt                                                                               <  2:498329/148‑1 (MQ=255)
                                                              ttCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCa                                                                             >  3:447996/1‑149 (MQ=255)
                                                              ttCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACCAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATc                                                                              >  1:637753/1‑148 (MQ=255)
                                                               tCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCa                                                                             >  1:291634/1‑148 (MQ=255)
                                                               tCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAg                                                                            >  4:97211/1‑149 (MQ=255)
                                                                cacaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAg                                                                            <  2:244534/148‑1 (MQ=255)
                                                                 acaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGtt                                                                          >  3:901180/1‑149 (MQ=255)
                                                                  caAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGtt                                                                          >  2:191311/1‑148 (MQ=255)
                                                                     aTCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTc                                                                       >  1:80844/1‑148 (MQ=255)
                                                                           ccAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGaa                                                                <  3:97211/149‑1 (MQ=255)
                                                                            cAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAg                                                               >  4:258767/1‑149 (MQ=255)
                                                                            cAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAg                                                               >  3:478477/1‑149 (MQ=255)
                                                                              gTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGt                                                             >  3:128576/1‑149 (MQ=255)
                                                                                  ttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTc                                                         <  3:484105/149‑1 (MQ=255)
                                                                                        tACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATcgc                                                                                                            <  3:196743/92‑1 (MQ=255)
                                                                                        tACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATcgc                                                                                                            >  4:196743/1‑92 (MQ=255)
                                                                                          cTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGTGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTga                                                 <  3:509683/149‑1 (MQ=255)
                                                                                          cTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTga                                                 <  4:124696/149‑1 (MQ=255)
                                                                                          cTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTg                                                  <  1:191311/148‑1 (MQ=255)
                                                                                           tGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTgag                                                >  3:874758/1‑149 (MQ=255)
                                                                                                  tAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATccc                                         <  4:901180/149‑1 (MQ=255)
                                                                                                    aGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGt                                                                                                                   <  4:269599/73‑1 (MQ=255)
                                                                                                    aGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGt                                                                                                                   >  3:269599/1‑73 (MQ=255)
                                                                                                    aGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCg                                       >  3:836603/1‑149 (MQ=255)
                                                                                                    aGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCg                                       >  3:263302/1‑149 (MQ=255)
                                                                                                      gAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGcc                                     >  4:95007/1‑149 (MQ=255)
                                                                                                      gAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGcc                                     <  4:103338/149‑1 (MQ=255)
                                                                                                       aaaTTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTa                                                                                                                >  1:174862/1‑73 (MQ=255)
                                                                                                       aaaTTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTa                                                                                                                <  2:174862/73‑1 (MQ=255)
                                                                                                         attaattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATc                                  >  3:185837/1‑149 (MQ=255)
                                                                                                          ttaattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGttt                                                           >  1:435069/1‑123 (MQ=255)
                                                                                                          ttaattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGttt                                                           <  2:435069/123‑1 (MQ=255)
                                                                                                           taattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATc                                                                              >  3:367880/1‑103 (MQ=255)
                                                                                                           taattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATc                                                                              <  4:367880/103‑1 (MQ=255)
                                                                                                           taattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGa                                <  4:128576/149‑1 (MQ=255)
                                                                                                           taattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGa                                <  4:644252/149‑1 (MQ=255)
                                                                                                            aattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAg                               <  3:468213/149‑1 (MQ=255)
                                                                                                             attTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAg                               >  2:106706/1‑148 (MQ=255)
                                                                                                                   gTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATaaa                         >  1:828032/1‑148 (MQ=255)
                                                                                                                       tGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCACGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATc                     <  2:529279/148‑1 (MQ=255)
                                                                                                                         cAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAg                   >  1:293975/1‑148 (MQ=255)
                                                                                                                          aGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGcc                 >  4:7094/1‑149 (MQ=255)
                                                                                                                          aGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGcc                 <  4:804453/149‑1 (MQ=255)
                                                                                                                          aGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGc                  <  1:106706/148‑1 (MQ=255)
                                                                                                                          aGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGc                  <  1:803650/148‑1 (MQ=255)
                                                                                                                             aaaaaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGt               >  1:724121/1‑148 (MQ=255)
                                                                                                                                aaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCaa           >  3:87784/1‑149 (MQ=255)
                                                                                                                                aaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCa            <  2:522396/148‑1 (MQ=255)
                                                                                                                                      gTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAAtt      >  1:126845/1‑148 (MQ=255)
                                                                                                                                       tCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTc    <  4:478477/149‑1 (MQ=255)
                                                                                                                                        cATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTcc   <  4:447996/149‑1 (MQ=255)
                                                                                                                                          tCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGtttt                                                          >  2:456820/1‑92 (MQ=255)
                                                                                                                                          tCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGtttt                                                          <  1:456820/92‑1 (MQ=255)
                                                                                                                                          tCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCt  >  2:55318/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                
CGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCT  >  NC_000913/2838081‑2838366

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: