Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,838,231 A→T 8.3% intergenic (‑126/+23) 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,2310AT 8.3% 284.3 / 14.1 84intergenic (‑126/+23)hydN/ascGformate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/asc operon transcriptional repressor; prpBC operon repressor
Reads supporting (aligned to +/- strand):  ref base A (37/40);  new base T (5/2);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01

CAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATC  >  NC_000913/2838084‑2838374
                                                                                                                                                   |                                                                                                                                               
cAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGAc                                                                                                                                                <  1:126937/149‑1 (MQ=255)
 gaatttgatgaACCGTTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACg                                                                                                                                               >  2:709443/6‑149 (MQ=255)
    aaTGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGa                                                                                                                                              >  1:144885/1‑147 (MQ=255)
    aaTGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGa                                                                                                                                              <  2:144885/147‑1 (MQ=255)
    aaTGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATg                                                                                                                                            >  3:277203/1‑149 (MQ=255)
    aaTGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATg                                                                                                                                            >  1:660862/1‑149 (MQ=255)
    aaTGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATg                                                                                                                                            <  3:170295/149‑1 (MQ=255)
       gatgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACa                                                                                                                                         <  4:114501/149‑1 (MQ=255)
        atgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACCAGTGACGATGACAt                                                                                                                                        <  3:54233/149‑1 (MQ=255)
          gaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgt                                                                                                                                      <  3:675560/149‑1 (MQ=255)
           aaaCGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgtg                                                                                                                                     <  3:251552/149‑1 (MQ=255)
            aaCGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtga                                                                                                                                    <  3:118500/149‑1 (MQ=255)
                gTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGAc                                                                                                                                >  4:385957/1‑149 (MQ=255)
                  tcatcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGAt                                                                                                                                             >  3:71772/1‑134 (MQ=255)
                  tcatcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGAt                                                                                                                                             <  4:71772/134‑1 (MQ=255)
                  tcatcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACaa                                                                                                                              >  3:289141/1‑149 (MQ=255)
                   catcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAAt                                                                                                                             <  3:549197/149‑1 (MQ=255)
                       aaaTTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGAt                                                                                                                                             <  1:658095/129‑1 (MQ=255)
                       aaaTTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGAt                                                                                                                                             >  2:658095/1‑129 (MQ=255)
                          ttCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATc                                                                                                                      <  3:90106/149‑1 (MQ=255)
                             tcCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATcgcg                                                                                                                   >  4:590634/1‑149 (MQ=255)
                              ccAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGa                                                                                                                                           <  4:743327/124‑1 (MQ=255)
                              ccAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGa                                                                                                                                           >  3:743327/1‑124 (MQ=255)
                                 gCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAgg                                                                                                               >  3:694905/1‑149 (MQ=255)
                                  cAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGa                                                                                                              >  3:370253/1‑149 (MQ=255)
                                   aaTGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAg                                                                                                             <  3:382082/149‑1 (MQ=255)
                                       aCAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAt                                                                                                         <  4:272588/149‑1 (MQ=255)
                                       aCAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAt                                                                                                         >  1:87763/1‑149 (MQ=255)
                                       aCAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAt                                                                                                         <  3:640803/149‑1 (MQ=255)
                                       aCAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGACGACATGTGATGACAATGTTTATCGCGAAGGAGCAAt                                                                                                         <  4:551316/149‑1 (MQ=255)
                                         aGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtga                                                                                                       <  2:87763/149‑1 (MQ=255)
                                           ttAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagt                                                                                                     <  1:517282/149‑1 (MQ=255)
                                           ttAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagt                                                                                                     <  2:261765/149‑1 (MQ=255)
                                           ttAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagt                                                                                                     <  2:500161/149‑1 (MQ=255)
                                             aaTGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagtga                                                                                                   <  2:290595/149‑1 (MQ=255)
                                                       tACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGcggcgg                                                                                         <  2:308272/149‑1 (MQ=255)
                                                          aTTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATc                                                                                      >  1:213600/1‑149 (MQ=255)
                                                             cacaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGt                                                                                   >  3:737665/1‑149 (MQ=255)
                                                              acaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGtt                                                                                  <  4:694905/149‑1 (MQ=255)
                                                              acaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGtt                                                                                  >  4:583466/1‑149 (MQ=255)
                                                              acaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGtt                                                                                  >  3:416540/1‑149 (MQ=255)
                                                                  aTCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGtt                                                                                  <  1:189440/145‑1 (MQ=255)
                                                                  aTCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGtt                                                                                  >  2:189440/1‑145 (MQ=255)
                                                                        ccAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGaa                                                                        <  2:367363/149‑1 (MQ=255)
                                                                         cAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCAGCTGAAg                                                                       <  2:213600/149‑1 (MQ=255)
                                                                           gTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAg                                                                                    >  1:131457/1‑134 (MQ=255)
                                                                           gTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAg                                                                                    <  2:131457/134‑1 (MQ=255)
                                                                           gTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGt                                                                     >  2:86274/1‑149 (MQ=255)
                                                                            tttttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGtt                                                                    >  1:538126/1‑149 (MQ=255)
                                                                             ttttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGttt                                                                   <  4:289141/149‑1 (MQ=255)
                                                                               ttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTc                                                                 <  3:188772/149‑1 (MQ=255)
                                                                               ttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTc                                                                 <  3:385957/149‑1 (MQ=255)
                                                                                 aaTTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTcgg                                                               >  1:651574/1‑149 (MQ=255)
                                                                                      aCTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTg                                                          <  4:277203/149‑1 (MQ=255)
                                                                                        tGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTgag                                                        >  3:41803/1‑149 (MQ=255)
                                                                                          ttATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGaa                                                      >  3:314196/1‑149 (MQ=255)
                                                                                            aTTTAAGGAAATTAATTTATGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAAt                                                    >  2:601178/1‑149 (MQ=255)
                                                                                                aaGGACATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATcccc                                                >  4:55417/1‑149 (MQ=255)
                                                                                                aaGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATcccc                                                >  3:613291/1‑149 (MQ=255)
                                                                                                  ggAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATcc                                                  <  2:629731/145‑1 (MQ=255)
                                                                                                  ggAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATcc                                                  >  1:629731/1‑145 (MQ=255)
                                                                                                     aattaattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCAt                                           >  1:155583/1‑149 (MQ=255)
                                                                                                       ttaattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCg                                         >  2:144019/1‑149 (MQ=255)
                                                                                                        taattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGa                                        <  3:692794/149‑1 (MQ=255)
                                                                                                         aattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAg                                       <  3:583466/149‑1 (MQ=255)
                                                                                                               tGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATaaa                                 <  2:306927/149‑1 (MQ=255)
                                                                                                                gTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATaaaa                                >  2:580382/1‑149 (MQ=255)
                                                                                                                    tGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCa                            <  3:172368/149‑1 (MQ=255)
                                                                                                                    tGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCa                            <  3:341522/149‑1 (MQ=255)
                                                                                                                     gCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAg                           >  2:183617/1‑149 (MQ=255)
                                                                                                                      cAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGc                          >  4:570531/1‑149 (MQ=255)
                                                                                                                      cAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGATTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGc                          >  4:423535/1‑149 (MQ=255)
                                                                                                                       aGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGcc                         >  1:635604/1‑149 (MQ=255)
                                                                                                                         gAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACATTGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGt                       >  4:112366/1‑149 (MQ=255)
                                                                                                                          tcatcacaTGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTc                      >  3:365669/8‑149 (MQ=255)
                                                                                                                             aaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCaa                   <  2:651574/149‑1 (MQ=255)
                                                                                                                                caTGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCaataa                >  1:705574/2‑149 (MQ=255)
                                                                                                                                 aTGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGATTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCaataat               <  3:423535/149‑1 (MQ=255)
                                                                                                                                     cATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTcc           >  1:595286/1‑149 (MQ=255)
                                                                                                                                      aTCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCt          >  1:150626/1‑149 (MQ=255)
                                                                                                                                      aTCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCt          >  3:16630/1‑149 (MQ=255)
                                                                                                                                            acTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGGTTTCCGCTGAAGGTTTTCGGCGTTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAatca    >  2:471121/1‑149 (MQ=255)
                                                                                                                                              tAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAatcatc  <  2:595286/149‑1 (MQ=255)
                                                                                                                                              tAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAAGTTCCTGAatcatc  <  1:659052/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATC  >  NC_000913/2838084‑2838374

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: