Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,390,985 C→T E175E (GAG→GAA rsgA ← ribosome small subunit‑dependent GTPase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,390,9850CT98.7% 310.4 / NA 77E175E (GAG→GAArsgAribosome small subunit‑dependent GTPase A
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (32/44);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCAGGCTGGATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCGGAGACAAT  >  NC_000913/4390837‑4391128
                                                                                                                                                    |                                                                                                                                               
tCAGCAGGCTGGATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTc                                                                                                                                               >  2:1026625/1‑151 (MQ=255)
 cagcagGCTGGATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCg                                                                                                                                              <  1:1060564/151‑1 (MQ=255)
    cagGCTGGATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTg                                                                                                                                           <  2:299622/151‑1 (MQ=255)
    cagGCTGGATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTCAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTg                                                                                                                                           >  1:122583/1‑151 (MQ=255)
      ggCTGGATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCGATCTGCTCGTTTAc                                                                                                                                         >  2:891022/1‑151 (MQ=255)
      ggCTGGATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGAc                                                                                                                                         >  1:482215/1‑151 (MQ=255)
        cTGGATTTGCCGACGCCAGACTGCCCGGCAACAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACg                                                                                                                                        >  1:865689/1‑150 (MQ=255)
               tGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCa                                                                                                                                 <  2:128434/150‑1 (MQ=255)
               tGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCAt                                                                                                                                >  1:711994/1‑151 (MQ=255)
                gCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATg                                                                                                                               <  1:167324/151‑1 (MQ=255)
                 ccGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGc                                                                                                                              >  1:738523/1‑151 (MQ=255)
                 ccGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGc                                                                                                                              >  2:408981/1‑151 (MQ=255)
                    aCGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCtt                                                                                                                           >  2:1015830/1‑151 (MQ=255)
                     cGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCtt                                                                                                                           <  2:580427/150‑1 (MQ=255)
                     cGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTtc                                                                                                                          >  1:868671/1‑151 (MQ=255)
                          gACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTcgtcgt                                                                                                                     <  1:619076/151‑1 (MQ=255)
                           aCTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCgtcgtc                                                                                                                    >  1:476764/1‑151 (MQ=255)
                            cTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTcc                                                                                                                   >  1:701539/1‑151 (MQ=255)
                              gCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTcc                                                                                                                   <  2:508518/149‑1 (MQ=255)
                              gCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCca                                                                                                                  <  1:843148/150‑1 (MQ=255)
                              gCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCca                                                                                                                  <  2:407447/150‑1 (MQ=255)
                               cccGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCcagc                                                                                                                <  1:666107/151‑1 (MQ=255)
                               cccGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCcagc                                                                                                                >  2:401006/1‑151 (MQ=255)
                                ccGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCcagca                                                                                                               <  1:310304/151‑1 (MQ=255)
                                 cGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCcagcag                                                                                                              <  1:991631/151‑1 (MQ=255)
                                 cGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCcagcag                                                                                                              <  1:343018/151‑1 (MQ=255)
                                   gCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGt                                                                                                            <  1:447917/151‑1 (MQ=255)
                                   gCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGt                                                                                                            <  1:952143/151‑1 (MQ=255)
                                    cAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTc                                                                                                           <  1:36971/151‑1 (MQ=255)
                                      aaaaTGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCt                                                                                                          >  2:226149/1‑150 (MQ=255)
                                           gCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCttg                                                                                                    <  1:933704/151‑1 (MQ=255)
                                              gATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCttgtt                                                                                                  <  1:164669/150‑1 (MQ=255)
                                              gATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCttgt                                                                                                   <  1:484781/149‑1 (MQ=255)
                                                tGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCttgttg                                                                                                 >  1:361855/1‑149 (MQ=255)
                                                   gcCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGcac                                                                                            <  2:743231/151‑1 (MQ=255)
                                                         tCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACaataat                                                                                      >  1:968380/1‑151 (MQ=255)
                                                            aCGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCgg                                                                                   >  1:480345/1‑151 (MQ=255)
                                                            aCGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCgg                                                                                   >  2:691086/1‑151 (MQ=255)
                                                                ctctTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTc                                                                                >  1:490858/1‑150 (MQ=255)
                                                                    tCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCaa                                                                              <  2:845469/148‑1 (MQ=255)
                                                                    tCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAAt                                                                             <  2:63662/149‑1 (MQ=255)
                                                                     ccAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCt                                                                           <  2:953240/150‑1 (MQ=255)
                                                                      cAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCt                                                                           <  2:764605/149‑1 (MQ=255)
                                                                      cAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGc                                                                         >  1:252915/1‑151 (MQ=255)
                                                                      cAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGc                                                                         <  2:37810/151‑1 (MQ=255)
                                                                        gCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCaa                                                                       <  1:1030625/151‑1 (MQ=255)
                                                                         cGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAg                                                                      <  1:987146/151‑1 (MQ=255)
                                                                             tttAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGtt                                                                   <  2:1020662/150‑1 (MQ=255)
                                                                              ttAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGttt                                                                  <  1:899149/150‑1 (MQ=255)
                                                                                    cATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGc                                                           <  1:695547/151‑1 (MQ=255)
                                                                                     aTCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGcc                                                          <  2:445838/151‑1 (MQ=255)
                                                                                      tCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGcc                                                          >  1:289903/1‑150 (MQ=255)
                                                                                       ccTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGcca                                                         <  2:777591/150‑1 (MQ=255)
                                                                                       ccTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGcca                                                         <  2:680470/150‑1 (MQ=255)
                                                                                           aGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAgg                                                    >  1:984511/1‑151 (MQ=255)
                                                                                            gTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGt                                                   <  1:825337/151‑1 (MQ=255)
                                                                                              aTGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTaa                                                 <  1:297040/151‑1 (MQ=255)
                                                                                              aTGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTa                                                  <  2:151611/150‑1 (MQ=255)
                                                                                               tGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAAc                                                >  1:500507/1‑151 (MQ=255)
                                                                                                    ggAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCg                                           >  1:505770/1‑151 (MQ=255)
                                                                                                     gAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCg                                           >  2:1000067/1‑150 (MQ=255)
                                                                                                        aCCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGataa                                       <  2:187335/151‑1 (MQ=255)
                                                                                                        aCCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGataa                                       <  1:142246/151‑1 (MQ=255)
                                                                                                             cAATACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTg                                  >  2:1059888/1‑151 (MQ=255)
                                                                                                              aaTACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTg                                  >  1:1037323/1‑150 (MQ=255)
                                                                                                              aaTACGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGa                                 <  2:738498/151‑1 (MQ=255)
                                                                                                                 aCGCGATAACCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGc                               <  1:401016/150‑1 (MQ=255)
                                                                                                                         aCCGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTcc                      <  1:457562/151‑1 (MQ=255)
                                                                                                                           cGATATTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCgg                    >  1:524956/1‑151 (MQ=255)
                                                                                                                             atatTGCGGTAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGc                   >  1:818110/1‑150 (MQ=255)
                                                                                                                                      tAGATATCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCgg         >  2:837637/1‑151 (MQ=255)
                                                                                                                                         atatCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCGgag       <  2:653763/150‑1 (MQ=255)
                                                                                                                                           atCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCGGAgaca    <  2:794511/151‑1 (MQ=255)
                                                                                                                                           atCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCGGAgac     >  1:189494/1‑150 (MQ=255)
                                                                                                                                            tCCATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCGGAgaca    >  2:123354/1‑150 (MQ=255)
                                                                                                                                             ccATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCGGAgacaat  <  2:681845/151‑1 (MQ=255)
                                                                                                                                             ccATCTGTTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCGGAgacaat  >  2:348110/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TCAGCAGGCTGGATTTGCCGACGCCAGACTGCCCGGCAAAAATGCTGATGCGCCCGGTCAACGCCTCTTCCAGCGGTTTTAGCCCATCCTGAGTATGGCTGGAAACCATCAATACGCGATAACCGATATTGCGGTAGATATCCATCTGCTCGTTGACGAACGCCATGCCTTCGTCGTCCAGCAGGTCTATCTTGTTGAGCACAATAATCGGCTCAATCTGCAAGGTTTCGCAGGCCACCAGGTAACGGTCGATAATATTGAGCGACAGCTCCGGCAAAATGGCGGAGACAAT  >  NC_000913/4390837‑4391128

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