Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,736,097 T→A 7.2% intergenic (+138/+24) sodB → / ← ydhP superoxide dismutase, Fe/putative MFS transporter, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,736,0970TA 7.2% 219.7 / 7.3 83intergenic (+138/+24)sodB/ydhPsuperoxide dismutase, Fe/putative MFS transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base A (2/4);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

CTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATCCCAGCCCAGCG  >  NC_000913/1735952‑1736241
                                                                                                                                                 |                                                                                                                                                
cTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATc                                                                                                                                              >  1:1017898/1‑150 (MQ=255)
 tGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATcc                                                                                                                                             >  2:1021880/1‑150 (MQ=255)
 tGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAAcc                                                                                                                                             <  1:467008/150‑1 (MQ=255)
 tGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAAcc                                                                                                                                             <  1:479626/150‑1 (MQ=255)
 tGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAAcc                                                                                                                                             <  1:1195584/150‑1 (MQ=255)
    ataataACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTc                                                                                                                                          <  1:474169/150‑1 (MQ=255)
    ataataACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGGCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTc                                                                                                                                          <  1:841937/150‑1 (MQ=255)
       ataaCTGCTGGCAACTGCAGCATTGCCTGTAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTCGGTCGGAGAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGc                                                                                                                                       <  1:204689/150‑1 (MQ=255)
       ataaCTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGc                                                                                                                                       <  2:338024/150‑1 (MQ=255)
       ataaCTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGAGCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGc                                                                                                                                       <  1:1261194/150‑1 (MQ=255)
        taaCTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCg                                                                                                                                      <  1:741497/150‑1 (MQ=255)
        taaCTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCg                                                                                                                                      <  1:674810/150‑1 (MQ=255)
           cTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGttt                                                                                                                                   >  2:344588/1‑150 (MQ=255)
            tGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATATATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGtttt                                                                                                                                  >  1:484954/1‑150 (MQ=255)
            tGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGtttt                                                                                                                                  <  2:377658/150‑1 (MQ=255)
              aTGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTg                                                                                                                                <  1:350345/150‑1 (MQ=255)
                ggCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGGTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGtt                                                                                                                              >  2:926702/1‑150 (MQ=255)
                        cagcaTTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGc                                                                                                                      <  2:1117492/150‑1 (MQ=255)
                            aTTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGtt                                                                                                                  <  2:17909/150‑1 (MQ=255)
                             ttGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTa                                                                                                                 <  2:509446/150‑1 (MQ=255)
                             ttGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTa                                                                                                                 <  1:971184/150‑1 (MQ=255)
                              tGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAATTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAg                                                                                                                >  1:925085/1‑150 (MQ=255)
                               gCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGc                                                                                                               >  1:472729/1‑150 (MQ=255)
                               gCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGc                                                                                                               <  2:913848/150‑1 (MQ=255)
                                 cTGAAGCGCTACGCTTATCATGCCTACGCGGATCAGAGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCaa                                                                                                             <  1:917260/150‑1 (MQ=255)
                                  tGAAGCTCTACGCTTATCAGGACTACGCGGATCATCGATGGAGGTCTGATAAGTCACTCGCCGCATCCGTCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAAc                                                                                                            <  1:350784/150‑1 (MQ=255)
                                     aGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTTTCCTTCTGCGTTTTTGTTTAATTATCTGTTAGCAACGCa                                                                                                         >  2:89131/1‑150 (MQ=255)
                                        gcTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAAc                                                                                                      >  2:917919/1‑150 (MQ=255)
                                          tACGCTTATCAGTCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTg                                                                                                    <  1:247731/150‑1 (MQ=255)
                                           aCGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGt                                                                                                   >  2:1173080/1‑150 (MQ=255)
                                            cGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGtt                                                                                                  <  1:432714/150‑1 (MQ=255)
                                              cTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTc                                                                                                >  1:178927/1‑150 (MQ=255)
                                              cTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTc                                                                                                >  1:813223/1‑150 (MQ=255)
                                                 aTCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTTATTAGCTGTTAGCAACGCAAACTGTTTCAgg                                                                                             >  2:1259182/1‑150 (MQ=255)
                                                 aTCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAgg                                                                                             >  2:1261372/1‑150 (MQ=255)
                                                 aTCACGCCTACGCGGATCAACGATGTAGGTCGGATAAGGCACTCGACGCATACGGCAAGATAAATCGCTCGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTATTTGTTAAATTAGCTGTTAGCAACGCAAACTGTTTCAgg                                                                                             <  1:56614/150‑1 (MQ=255)
                                                  tCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgt                                                                                            <  2:747532/150‑1 (MQ=255)
                                                  tCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgt                                                                                            <  1:759535/150‑1 (MQ=255)
                                                   cAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgtt                                                                                           >  2:45777/1‑150 (MQ=255)
                                                    aGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGAACCGCAAACTGTTTCAGgttg                                                                                          >  2:719282/1‑150 (MQ=255)
                                                      gCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttgtt                                                                                        >  2:596677/1‑150 (MQ=255)
                                                         tACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTc                                                                                     <  2:1252165/150‑1 (MQ=255)
                                                         tACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCt                                                                                    >  1:480323/1‑150 (MQ=255)
                                                          aCGCGGATCATCGATGTATGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCt                                                                                    <  1:93439/150‑1 (MQ=255)
                                                            gcgGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTgg                                                                                  >  1:1059190/1‑150 (MQ=255)
                                                               gATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTg                                                                               >  2:918442/1‑150 (MQ=255)
                                                                   atcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACAt                                                                           >  1:1194278/1‑150 (MQ=255)
                                                                      gATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATaaa                                                                        >  1:602827/1‑150 (MQ=255)
                                                                       aTGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTGTATTTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAac                                                                       >  2:1160064/1‑150 (MQ=255)
                                                                         gTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGAAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTAAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAacac                                                                     <  2:847001/150‑1 (MQ=255)
                                                                            ggTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAg                                                                  <  2:956799/150‑1 (MQ=255)
                                                                                  aTAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCaataa                                                            <  2:508098/150‑1 (MQ=255)
                                                                                  aTAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAAATAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCaataa                                                            >  2:873073/1‑150 (MQ=255)
                                                                                   tAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCaataat                                                           <  1:1205706/150‑1 (MQ=255)
                                                                                    aaGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGCAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCTATAATg                                                          >  1:189663/1‑150 (MQ=255)
                                                                                     agtcACTCGCCGCATCCGGCAAGTTAATTCGGACGTTGTCTGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAATCTGTTTCAGGTTGTTTTCTGGGTGACATAAACAGCTGCTATAATGc                                                         <  2:1200989/147‑1 (MQ=255)
                                                                                      ttcACTCTCCGCATCCGGCAATATAAATCTCACCTTGTCAGCAACTGTAACGAAGAAGGGTATCCTTCTGCGTTTTAGTTTAATTAGCAGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGcc                                                        <  2:638556/148‑1 (MQ=255)
                                                                                             gccgcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCACTAATGCCAGTCCcg                                                 >  1:863716/1‑150 (MQ=255)
                                                                                              ccgcACCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTTTCCTTCTGCGTTTTTGGTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCcgc                                                >  2:668579/1‑150 (MQ=255)
                                                                                                 cATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGCAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAATCACCAGCAATAATGCCAGTCCCGCGAc                                             >  2:1103944/1‑150 (MQ=255)
                                                                                                     cGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATc                                         <  1:628873/150‑1 (MQ=255)
                                                                                                         aaGATTAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTc                                     <  2:169274/150‑1 (MQ=255)
                                                                                                         aaGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTc                                     <  2:526767/150‑1 (MQ=255)
                                                                                                          aGATAAAAAGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGGTAGCAACGCAAACTTTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTcc                                    <  2:887466/150‑1 (MQ=255)
                                                                                                           gATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTccc                                   >  2:967155/1‑150 (MQ=255)
                                                                                                               aaTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCa                               >  2:644693/1‑150 (MQ=255)
                                                                                                                aTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCAc                              >  2:231600/1‑150 (MQ=255)
                                                                                                                  cGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCg                            >  1:765261/1‑150 (MQ=255)
                                                                                                                  cGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCg                            >  1:890441/1‑150 (MQ=255)
                                                                                                                           tCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGc                   <  1:344588/150‑1 (MQ=255)
                                                                                                                            cAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCt                  <  2:1068909/150‑1 (MQ=255)
                                                                                                                            cAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCt                  >  1:872319/1‑150 (MQ=255)
                                                                                                                             aGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTg                 >  2:619206/1‑150 (MQ=255)
                                                                                                                                aaCTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTAt              >  1:994776/1‑150 (MQ=255)
                                                                                                                                  caaaaaCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATcc            >  1:1111759/5‑150 (MQ=255)
                                                                                                                                  cTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCACCGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCCCCGGCACAAAGCTGTATcc            >  1:295928/1‑150 (MQ=255)
                                                                                                                                   tGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATccc           <  2:80337/150‑1 (MQ=255)
                                                                                                                                       aCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATcccagcc       >  1:371253/1‑150 (MQ=255)
                                                                                                                                        cGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATcccagccc      <  1:1313159/150‑1 (MQ=255)
                                                                                                                                        cGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATcccagccc      <  1:918442/150‑1 (MQ=255)
                                                                                                                                            gAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATCCCAGCCCAGCg  <  2:588772/150‑1 (MQ=255)
                                                                                                                                            gAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATCCCAGCCCAGCg  <  1:612356/150‑1 (MQ=255)
                                                                                                                                            gAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGATCCCATCCCCGGCACTAAGCTGTATCCCAGCCCAGCg  >  2:1288521/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATCCCAGCCCAGCG  >  NC_000913/1735952‑1736241

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: