Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,666,468 A→T 6.4% Q385L (CAG→CTG)  clcB → H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,666,4680AT 6.4% 244.7 / 4.8 88Q385L (CAG→CTG) clcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base A (38/44);  new base T (2/4);  total (40/48)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

GTATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATAA  >  NC_000913/1666319‑1666617
                                                                                                                                                     |                                                                                                                                                     
gtATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGTCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGAGGCTTTGTGATATGACCGGGGAGTATCa                                                                                                                                                       >  2:677247/1‑150 (MQ=255)
gtATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCa                                                                                                                                                       <  1:465862/150‑1 (MQ=255)
gtATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCa                                                                                                                                                       <  1:599108/150‑1 (MQ=255)
gtATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCa                                                                                                                                                       <  1:781253/150‑1 (MQ=255)
 tATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAg                                                                                                                                                      >  1:474721/1‑150 (MQ=255)
 tATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAg                                                                                                                                                      <  1:796940/150‑1 (MQ=255)
  aTGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGc                                                                                                                                                     >  1:416557/1‑150 (MQ=255)
    ggTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTa                                                                                                                                                   >  2:1227476/1‑150 (MQ=255)
     ggcGTAGCCTGGGATTATGGTTCCCTGAGGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCAGCTAc                                                                                                                                                  <  2:499126/148‑1 (MQ=255)
        gTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcc                                                                                                                                               >  1:878497/1‑150 (MQ=255)
          aGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcccc                                                                                                                                             <  2:1532825/150‑1 (MQ=255)
          aGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcccc                                                                                                                                             <  2:1384513/150‑1 (MQ=255)
          aGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTcccc                                                                                                                                             >  2:1704572/1‑150 (MQ=255)
          aGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTcccc                                                                                                                                             <  2:228701/150‑1 (MQ=255)
            ccTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCgg                                                                                                                                           >  2:376760/1‑150 (MQ=255)
            ccTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                                                                           <  1:64792/150‑1 (MQ=255)
            ccTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCgg                                                                                                                                           <  1:727002/150‑1 (MQ=255)
                 gATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTat                                                                                                                                      >  2:1395576/1‑150 (MQ=255)
                 gATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATGGACCGGGATGGCGACACTGTTGGCGGCGCCCCCGCCCGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAGGCCCCGGTTTat                                                                                                                                      >  2:1448496/1‑150 (MQ=255)
                   ttATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattg                                                                                                                                    >  1:1666180/1‑150 (MQ=255)
                   ttATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGtcatttc                                                                                                                                    <  1:520252/150‑7 (MQ=255)
                     aTGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgat                                                                                                                                  <  2:1048958/150‑1 (MQ=255)
                       ggTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                                                                                <  1:1245305/150‑1 (MQ=255)
                       ggTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgattg                                                                                                                                <  1:508371/150‑1 (MQ=255)
                           cccTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTg                                                                                                                            >  1:1428862/1‑150 (MQ=255)
                             cTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCg                                                                                                                          >  1:580240/1‑150 (MQ=255)
                                   gCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTg                                                                                                                    <  1:1733528/150‑1 (MQ=255)
                                   gCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTg                                                                                                                    >  1:545344/1‑150 (MQ=255)
                                   ccGAAGAAATTACACTTTTACTCGGATTGACCGGGATTGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTg                                                                                                                    <  1:1241028/149‑1 (MQ=255)
                                     gaagaaATTACACTTTTACTCGGATTGCCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCTCCGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAATCCACGGTTTATTTATTTCCTTAGTACTTGCg                                                                                                                  >  2:490063/1‑150 (MQ=255)
                                       agaaATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCgcc                                                                                                                <  2:397268/150‑3 (MQ=255)
                                       agaaATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTc                                                                                                                <  1:207936/150‑1 (MQ=255)
                                        gaaATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCg                                                                                                               <  2:474721/150‑1 (MQ=255)
                                         aaaTTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCgg                                                                                                              >  1:785809/1‑150 (MQ=255)
                                         aaaTTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCgg                                                                                                              >  1:722776/1‑150 (MQ=255)
                                           aTTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTa                                                                                                            <  1:1599018/150‑1 (MQ=255)
                                           aTTACACTTTAACTCGGATTGAGCGGGAGGGCGACAATGTTGTCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGATATGACCGGGGAGTTTCAGCTAGTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTa                                                                                                            <  1:827635/150‑1 (MQ=255)
                                               cacTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAAttt                                                                                                        <  2:1247666/150‑1 (MQ=255)
                                                   tttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                                                                    <  1:1047193/150‑1 (MQ=255)
                                                    ttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                                                                                   <  2:416557/150‑1 (MQ=255)
                                                    ttACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCgg                                                                                                   <  1:1337319/150‑1 (MQ=255)
                                                       cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACg                                                                                                <  1:1605946/150‑1 (MQ=255)
                                                       cTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGGTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACg                                                                                                >  2:587184/1‑150 (MQ=255)
                                                           gcttGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTac                                                                                            <  2:1297895/148‑1 (MQ=255)
                                                           gATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGGTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGtttc                                                                                            >  1:1435799/1‑148 (MQ=255)
                                                               gACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGGTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCg                                                                                        >  1:1468044/1‑150 (MQ=255)
                                                                 ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTGCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTTCACCGTg                                                                                      >  1:125146/1‑150 (MQ=255)
                                                                 ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTg                                                                                      >  2:700969/1‑150 (MQ=255)
                                                                  cGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTgc                                                                                     >  2:920055/1‑149 (MQ=255)
                                                                   gggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                                                                    >  2:1725823/1‑150 (MQ=255)
                                                                   gggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                                                                    <  1:194814/150‑1 (MQ=255)
                                                                   gggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                                                                    <  1:1123756/150‑1 (MQ=255)
                                                                   gggATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAc                                                                                    >  2:1025655/1‑150 (MQ=255)
                                                                     gATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTCTCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGActc                                                                                  <  2:1615749/150‑1 (MQ=255)
                                                                      aTGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGActct                                                                                 <  2:1666180/150‑1 (MQ=255)
                                                                       tGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATTACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtcta                                                                                >  2:583760/1‑150 (MQ=255)
                                                                        ggCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctat                                                                               >  2:1582486/1‑150 (MQ=255)
                                                                          cGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatct                                                                             >  1:1666032/1‑150 (MQ=255)
                                                                           gACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatcta                                                                            <  1:822717/150‑1 (MQ=255)
                                                                            acacTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAc                                                                           <  2:106824/150‑1 (MQ=255)
                                                                                tGTTGGCGTCGACCACGCTGGCGCCGATTATGTCGACTTTGATGATATGTGAAATTAACGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCATGCGTATTTGCGTCGGTAATTTCGCGGACGTTACACCGTTACTCTATCTAcctca                                                                       <  2:752068/150‑4 (MQ=255)
                                                                                   tGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGc                                                                    <  1:1561822/150‑1 (MQ=255)
                                                                                         cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg                                                              >  1:1557018/1‑150 (MQ=255)
                                                                                         cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTg                                                              <  1:1345698/150‑1 (MQ=255)
                                                                                         cgACCACGCACGCGCCGATTATGTCGACGTTCATGATATGTGAAATGACCCGGGAGTATCAGCTACTCCCCGGTTTATAGATTGCCTGCCTAATTGCGTAGGTAATTTCGCGGACCATACACCTTCACTCTAACAACCGCACACACACtc                                                              >  1:320530/1‑149 (MQ=255)
                                                                                          gACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgg                                                             <  1:700969/150‑2 (MQ=255)
                                                                                              acgcacgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag                                                         <  1:376760/150‑1 (MQ=255)
                                                                                                     cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCt                                                  >  1:1271606/1‑150 (MQ=255)
                                                                                                          aTTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCg                                             >  2:563136/1‑150 (MQ=255)
                                                                                                           ttATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGGCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCgg                                            <  2:64251/150‑1 (MQ=255)
                                                                                                              tGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATg                                         >  1:608720/1‑150 (MQ=255)
                                                                                                               gTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCTGATGt                                        >  2:1645580/1‑150 (MQ=255)
                                                                                                                    cGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGc                                   <  2:1422087/150‑1 (MQ=255)
                                                                                                                     gTTGATGATATGTGACAAGACCGGGGAGAATCAGCAACTCCCCGGTTTATTGATTGCCTGCGTACTTGCGTCGGTAATATCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGcc                                  <  2:199052/150‑1 (MQ=255)
                                                                                                                     gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGcc                                  <  1:1239215/150‑1 (MQ=255)
                                                                                                                      ttGATGATATGTGAAATGACCGGGGAGGAGCAGCTACTCCCCGGTTTATTGATTGCCGGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCCCACTGCGCAGCATAGCTAAGCGGATGTATTGccc                                 <  1:490063/150‑1 (MQ=255)
                                                                                                                       tgatgaTATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGcccc                                <  1:833571/150‑1 (MQ=255)
                                                                                                                       tgatgaTATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGcccc                                >  2:548482/1‑150 (MQ=255)
                                                                                                                          tgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGt                             >  1:1157529/1‑150 (MQ=255)
                                                                                                                              atGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCCGCACCCTGCGCAGCAGAGCTAAGCGGATGTATTGCCCCAGTTCAc                         >  2:844813/1‑150 (MQ=255)
                                                                                                                                gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGc                       <  1:1025655/150‑1 (MQ=255)
                                                                                                                                gtgAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGAATAGCTAAGCGGATGTATTGCCCCAGTTCACGc                       <  2:591605/150‑1 (MQ=255)
                                                                                                                                   aaaTGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGt                    <  1:717920/150‑1 (MQ=255)
                                                                                                                                      tGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCg                 <  1:157379/150‑1 (MQ=255)
                                                                                                                                       gACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCgg                <  1:1536610/150‑1 (MQ=255)
                                                                                                                                         ccGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgcg              >  1:1307510/1‑150 (MQ=255)
                                                                                                                                           ggggAGTAGCTGATACTCCCCGGTTTATTCATTCCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTCCACCGTGACTCTATCTACCGCAAGCACACTGCGCAGCTTAGATAAGCCGATCTATTGCCCAAGTTCACGCTGTTCGgcgcg            >  2:276653/1‑150 (MQ=255)
                                                                                                                                               aGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGc        >  1:965087/1‑150 (MQ=255)
                                                                                                                                                gTAGCAGCTACTCCGCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCa       <  1:1663866/150‑1 (MQ=255)
                                                                                                                                                  aTCAGCTACTCCCCGGTTTATTGAATGCCTGCGTAATTGCGTCGGTAAATTCCCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCACAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGa     >  1:197410/1‑150 (MQ=255)
                                                                                                                                                   tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGAt    >  1:1212155/1‑150 (MQ=255)
                                                                                                                                                   tCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGAt    <  1:1476183/150‑1 (MQ=255)
                                                                                                                                                     aGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATaa  >  2:458121/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTATGGTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATAA  >  NC_000913/1666319‑1666617

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: