Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,594,232 T→C 5.9% N417D (AAC→GAC)  livM ← branched‑chain amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,594,2320TC 5.9% 208.6 / 3.6 68N417D (AAC→GAC) livMbranched‑chain amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (34/30);  new base C (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

AGGGCCGATTAACGAGACGATCTCCTGCGGGTACAGTTCAAGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATCACGCATCAACTCGCGCGACACCAC  >  NC_000913/3594090‑3594364
                                                                                                                                              |                                                                                                                                    
aGGGCCGATTAACGAGACGATCTCCTGCGGGTACAGTTCAAGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGttttt                                                                                                                                   >  3:148939/1‑146 (MQ=255)
     cGATTAACGAGACGATCTCCTGCGGGTACAGTTCAAGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTc                                                                                                                            <  5:322700/148‑1 (MQ=255)
      gATTAACGAGACGATCTCCTGCGGGTACAGTTCAAGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCtt                                                                                                                             >  5:95930/1‑146 (MQ=255)
      gATTAACGAGACGATCTCCTGCGGGTACAGTTCAAGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCtt                                                                                                                             >  2:240691/1‑146 (MQ=255)
       aTTAACGAGACGATCTCCTGCGGGTACAGTTCAAGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAg                                                                                                                          >  2:324883/1‑148 (MQ=255)
           aCGAGACGATCTCCTGCGGGTACAGTTCAAGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGc                                                                                                                      <  3:142747/148‑1 (MQ=255)
                aCGATCTCCTGCGGGTACAGTTCAAGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGgcg                                                                                                                 <  4:148939/148‑1 (MQ=255)
                     ctcCTGCGGGTACAGTTCAAGATTGACGTTGTTCACCGCGAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTAAGCTTCAGTTGCGGGCGCGTc                                                                                                             <  8:36384/147‑1 (MQ=255)
                        cTGCGGGTACAGTTCAAGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTc                                                                                                             <  1:240691/144‑1 (MQ=255)
                        cTGCGGGTACAGTTCAAGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTc                                                                                                             <  3:188114/144‑1 (MQ=255)
                               tACAGTTCAAGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCaa                                                                                                  <  6:95930/148‑1 (MQ=255)
                                aCAGTTCAAGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAg                                                                                                 >  2:207146/1‑148 (MQ=255)
                                    ttCAAGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATg                                                                                                          >  8:384165/1‑135 (MQ=255)
                                    ttCAAGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATg                                                                                                          <  7:384165/135‑1 (MQ=255)
                                     tCAAGATTGACGTTGTTCCCCGCCAGCAGGCCGCCGAACCGCATCAACAGGCCGATAACAGATAAAAATGGCTGACTCCTGCCTGCTCTCCTTTCGCTGCGCGGTTTTTCAGCTTCAGATGCGGGCGCGTCATGGGCAGCAAGCCCTg                                                                                            >  1:139047/1‑148 (MQ=255)
                                        aGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCgg                                                                                         <  1:207146/148‑1 (MQ=255)
                                                   gttCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGatc                                                                               >  6:80681/1‑147 (MQ=255)
                                                      cACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGatcatc                                                                            <  1:324883/147‑1 (MQ=255)
                                                            cagcagGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCa                                                                     >  5:353121/1‑148 (MQ=255)
                                                               cagGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCa                                                                  <  6:353121/148‑1 (MQ=255)
                                                               cagGCCGCCGAAGCGCATCAGCAGGCCGTTAACAGATAATAATGGCTGACGCATGCGTGGACTCCTTTCTCTGCGCCGTTTTTAAGCTTCAGTTGCGGGCGCGTCATGGGCAGCACGCCATGCGGAGGGCAGatcataatctgaccca                                                                  >  6:377860/1‑141 (MQ=255)
                                                                agGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCCGCAAGCCCTGCGGACGCCAGATCCTCCTCAGCACCAt                                                                 >  1:351198/1‑148 (MQ=255)
                                                                                 catcaGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGc                                                                                                >  5:323028/1‑100 (MQ=255)
                                                                                 catcaGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGc                                                                                                <  6:323028/100‑1 (MQ=255)
                                                                                 catcaGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGc                                                                                     >  1:352358/1‑111 (MQ=255)
                                                                                 catcaGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGc                                                                                     <  2:352358/111‑1 (MQ=255)
                                                                                 catcaGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCAtt                                                 >  1:93094/1‑147 (MQ=255)
                                                                                 catcaGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCAtt                                                 >  8:337359/1‑147 (MQ=255)
                                                                                 catcaGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCAtt                                                 >  3:51289/1‑147 (MQ=255)
                                                                                  atcaGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTaa                                               <  5:80681/148‑1 (MQ=255)
                                                                                    caGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTAAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGCCGCCAGATCATCATCAGCACAATCAAACCACCGAGAAGTAACa                                             >  6:198595/1‑148 (MQ=255)
                                                                                          gTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAaccac                                                         <  2:241445/130‑1 (MQ=255)
                                                                                          gTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAaccac                                                         >  1:241445/1‑130 (MQ=255)
                                                                                           ttAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTa                                      >  8:114423/1‑148 (MQ=255)
                                                                                            tAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACa                                             >  5:318254/1‑140 (MQ=255)
                                                                                            tAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTAt                                     >  2:347937/1‑148 (MQ=255)
                                                                                             aaCAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTc                                                                                                             >  4:237886/1‑75 (MQ=255)
                                                                                             aaCAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTc                                                                                                             <  3:237886/75‑1 (MQ=255)
                                                                                                   taataaTGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAAc                                              <  6:199618/132‑1 (MQ=255)
                                                                                                   taataaTGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAAc                                              >  5:199618/1‑132 (MQ=255)
                                                                                                   taataaTGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTg                               >  1:154940/1‑147 (MQ=255)
                                                                                                   taataaTGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTCTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTg                               >  4:276219/1‑147 (MQ=255)
                                                                                                   taataaTGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTCTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTg                               >  4:271402/1‑147 (MQ=255)
                                                                                                     ataatGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGaa                             <  2:93094/147‑1 (MQ=255)
                                                                                                     ataatGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTCTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGaa                             <  3:271402/147‑1 (MQ=255)
                                                                                                     ataatGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTCTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGaa                             <  3:276219/147‑1 (MQ=255)
                                                                                                               acgcATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATCACGCAt                    <  2:139047/143‑1 (MQ=255)
                                                                                                                cTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGAtcatcatca                                                                        <  5:16809/93‑1 (MQ=255)
                                                                                                                cTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGAtcatcatca                                                                        >  6:16809/1‑93 (MQ=255)
                                                                                                                  cATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGc                                                    >  3:16483/1‑111 (MQ=255)
                                                                                                                  cATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGc                                                    <  4:16483/111‑1 (MQ=255)
                                                                                                                   aTGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATCACGCATCAAc                <  7:337359/146‑1 (MQ=255)
                                                                                                                   aTGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATCACGCATCAAc                <  7:114423/146‑1 (MQ=255)
                                                                                                                   aTGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATCACGCATCAAc                <  1:347937/146‑1 (MQ=255)
                                                                                                                   aTGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATCACGCATCAACTc              >  6:276562/1‑148 (MQ=255)
                                                                                                                     gCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATCACGCATCAACTcgc            >  4:257847/1‑148 (MQ=255)
                                                                                                                     gCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATCACGCATCAACTcgc            >  7:222977/1‑148 (MQ=255)
                                                                                                                      ccTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATCACGCATCAACTcgc            <  3:257847/147‑1 (MQ=255)
                                                                                                                      acTGATCTACTTTAGATGCTCCTTTTTTCATATTCAGTTGAGGGCTCGTCATGGGAAGCAAGCCCTGCGGACTCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACCTGCTGTATTAGTTGCAATACCGCATCAACTcgcg           >  4:288349/2‑148 (MQ=255)
                                                                                                                         gCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTAt                                     <  8:152961/119‑1 (MQ=255)
                                                                                                                         gCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTAt                                     >  7:152961/1‑119 (MQ=255)
                                                                                                                         gCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGaa                             >  2:108899/1‑127 (MQ=255)
                                                                                                                         gCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGaa                             <  1:108899/127‑1 (MQ=255)
                                                                                                                          ctctCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATCACGCATCAACTCGCGCGac       <  8:222977/148‑1 (MQ=255)
                                                                                                                           tctcCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATCACGCATCAACTCGCGCGaca      >  7:272902/1‑148 (MQ=255)
                                                                                                                               cTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATCACGCATCAACTCGCGCGAcacca   <  4:51289/147‑1 (MQ=255)
                                                                                                                                tttCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATCACGCATCAACTCGCGCGAcaccac  >  1:335309/1‑147 (MQ=255)
                                                                                                                                   cGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATc                          <  3:312980/120‑1 (MQ=255)
                                                                                                                                   cGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATc                          >  4:312980/1‑120 (MQ=255)
                                                                                                                                              ttttCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAAt                           <  3:42444/108‑1 (MQ=255)
                                                                                                                                              ttttCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAAt                           >  4:42444/1‑108 (MQ=255)
                                                                                                                                              |                                                                                                                                    
AGGGCCGATTAACGAGACGATCTCCTGCGGGTACAGTTCAAGATTGACGTTGTTCACCGCCAGCAGGCCGCCGAAGCGCATCATCAGGCCGTTAACAGATAATAATGGCTGACTCATGCCTGCTCTCCTTTCGCTGCGCCGTTTTTCAGCTTCAGTTGCGGGCGCGTCATGGGCAGCAAGCCCTGCGGACGCCAGATCATCATCAGCACCATCAAACCACCGAGCATTAACATGCTGTATTCGTTGAAATCACGCATCAACTCGCGCGACACCAC  >  NC_000913/3594090‑3594364

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