Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,336,988 A→T 7.8% V15D (GTC→GAC)  adiC ← arginine:agmatine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,336,9880AT 7.8% 261.6 / 8.0 98V15D (GTC→GAC) adiCarginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base A (42/47);  new base T (6/2);  total (49/49)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTT  >  NC_000913/4336839‑4337137
                                                                                                                                                     |                                                                                                                                                     
aCGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCAGCAGGGTGa                                                                                                                                                       >  1:494044/1‑150 (MQ=255)
  gACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCCTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACCCCTGACCCCATAATATTCCCCGGCACCATCCGGGTGCCg                                                                                                                                                     >  2:452385/1‑150 (MQ=255)
    cATTTTTGCGTATACCATCTAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTTATGGTGTCggg                                                                                                                                                   <  1:1232079/150‑1 (MQ=255)
    cATTTTTGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCAAAATATTCCCCGACACCATCAGGGTGACggg                                                                                                                                                   <  1:1157204/150‑1 (MQ=255)
     aTTTTGTCGTATACCATCGAGCGCGCCAGTGCACCTATAATCGTCACCACCAATGCATAAATTGCACTCCCGCGAGTAGAGGCCAGGTATGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGTCACCCTGATGGTGTCgggg                                                                                                                                                  <  2:206785/150‑1 (MQ=255)
       ttacgCGTATACCAACCACTGCCCCACTTCACCGATAATCGTCAACTCCCATCCATAAATGGCAATCCCGCCAGTAGAGGCGAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATACTATTCCCCCACCCCATCAGGGTGACGGGGAt                                                                                                                                                <  2:78013/146‑1 (MQ=255)
         tGGCGTATACCATCTAGATCCCCTGTGCCCCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTAGCAGGTAACAGAAATACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTa                                                                                                                                              <  2:1212275/150‑1 (MQ=255)
         tGGCGTATACCATCGACAGCCCCAGTGCACCGATAATCGTCACCCACCATCCATAAATGGCAAGCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTa                                                                                                                                              <  1:761272/150‑1 (MQ=255)
          ggCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACCCCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTaa                                                                                                                                             >  2:26035/1‑150 (MQ=255)
          ggCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTaa                                                                                                                                             >  2:21571/1‑150 (MQ=255)
           gCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCTCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAg                                                                                                                                            >  2:1093474/1‑150 (MQ=255)
            cGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCACCCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGc                                                                                                                                           >  1:630933/1‑150 (MQ=255)
             gTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGcc                                                                                                                                          <  1:441990/150‑1 (MQ=255)
              tataCCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGccc                                                                                                                                         <  2:441567/150‑1 (MQ=255)
               ataCCATCGAGAGCCCCAGTGCACCGATAATCGTCACAAACCATCCATAAATTGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCa                                                                                                                                        <  2:1187167/150‑1 (MQ=255)
                taCCATCGAGCGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCATTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACCACTGACCCCATAAGATTCCCCGACACCATCAGGGTGACGGGGGTTACGccccc                                                                                                                                       >  1:176808/1‑148 (MQ=255)
                taCCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCAc                                                                                                                                       <  1:429539/150‑1 (MQ=255)
                   cATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACttt                                                                                                                                    >  2:425655/1‑150 (MQ=255)
                      cGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgtg                                                                                                                                 <  2:989301/150‑1 (MQ=255)
                        agagCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGag                                                                                                                               <  1:552144/150‑1 (MQ=255)
                        agagCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGag                                                                                                                               <  1:464662/150‑1 (MQ=255)
                        agagCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGACAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGag                                                                                                                               <  1:603815/150‑1 (MQ=255)
                          agCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagca                                                                                                                             >  1:198041/1‑150 (MQ=255)
                            ccccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatc                                                                                                                           >  2:770668/1‑150 (MQ=255)
                            ccccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatc                                                                                                                           <  2:1120734/150‑1 (MQ=255)
                                        gATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCACCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGa                                                                                                               >  2:184843/1‑150 (MQ=255)
                                            aTCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATa                                                                                                           <  2:383550/150‑1 (MQ=255)
                                            aTCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGGATGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGCAGACATa                                                                                                           <  2:1277909/150‑1 (MQ=255)
                                              cGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCCCTTTGTGAGCATCAGCATCCGAAGACATAAt                                                                                                         >  2:1265869/1‑150 (MQ=255)
                                              cGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAt                                                                                                         <  2:1266006/150‑1 (MQ=255)
                                               gTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAACACACCTGACCCCATAATATTCCCCGCCACCATCAGGGGGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAtt                                                                                                        >  2:164890/1‑150 (MQ=255)
                                               gTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGCCCCCATAATATTCCCCGACCCCATCAGGGAGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACAAAAtt                                                                                                        >  1:107558/1‑150 (MQ=255)
                                               gTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAtt                                                                                                        >  1:368882/1‑150 (MQ=255)
                                                tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTa                                                                                                       <  1:296941/150‑1 (MQ=255)
                                                 cagcaaccaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTaa                                                                                                      <  1:652186/147‑1 (MQ=255)
                                                 cACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTaa                                                                                                      >  1:359798/1‑150 (MQ=255)
                                                  accaaccaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTaaa                                                                                                     >  2:260785/1‑150 (MQ=255)
                                                    caaccaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACt                                                                                                   >  1:1051418/1‑150 (MQ=255)
                                                      accaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACGTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTcc                                                                                                 <  2:198041/150‑1 (MQ=255)
                                                        catccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTg                                                                                               <  1:1011752/150‑1 (MQ=255)
                                                         atccatAAATGGCCATCCCGCCAGTAGAGTCCAGGGTTGCAGGTAACAGAAACACACCTGACCCCATAATATTCCCCGACACCCTCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGc                                                                                              <  2:1110034/150‑1 (MQ=255)
                                                          tccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCg                                                                                             <  1:606408/150‑1 (MQ=255)
                                                            catAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGaa                                                                                           >  1:928339/1‑150 (MQ=255)
                                                            catAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGAAGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGaa                                                                                           >  1:686353/1‑150 (MQ=255)
                                                             atAAATGGCAATCCCGCCAGTAGCGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAg                                                                                          >  2:610841/1‑150 (MQ=255)
                                                                aaTGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCg                                                                                       <  2:1203151/150‑1 (MQ=255)
                                                                   ggCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGc                                                                                    <  1:974303/150‑1 (MQ=255)
                                                                    gCAATCCCGCCATTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCt                                                                                   <  1:150846/150‑1 (MQ=255)
                                                                     cAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCtt                                                                                  >  1:967381/1‑150 (MQ=255)
                                                                       aTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCg                                                                                <  2:718694/150‑1 (MQ=255)
                                                                        tCCCGCCAGTAGATGCCAGGTTTGCAGGTAACAGTAAAACACCTGACCCCATAATCTTCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                                              <  2:594389/150‑1 (MQ=255)
                                                                        tCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGt                                                                               <  2:651799/150‑1 (MQ=255)
                                                                        tCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGt                                                                               <  1:1102052/150‑1 (MQ=255)
                                                                           cGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                                            <  2:598378/150‑1 (MQ=255)
                                                                           cGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                                            <  1:876347/150‑1 (MQ=255)
                                                                              cAGTAGAGGCCAGGTTTGCAGGTAACAGATCAACACCTGACCCCATAATATTCCCCGACACCATCAGGGGGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAAt                                                                         <  1:1070465/150‑1 (MQ=255)
                                                                               aGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATa                                                                        >  2:946530/1‑150 (MQ=255)
                                                                                  agagGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGCCACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAAt                                                                     >  2:865010/1‑150 (MQ=255)
                                                                                     ggCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGAGAGCTTCGTGACAATAAATTAc                                                                  >  1:8443/1‑150 (MQ=255)
                                                                                      gCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACg                                                                 <  2:630933/150‑1 (MQ=255)
                                                                                       ccAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGt                                                                >  1:137544/1‑150 (MQ=255)
                                                                                       ccAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGt                                                                >  1:1054820/1‑150 (MQ=255)
                                                                                        cAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                                               <  2:1039932/150‑1 (MQ=255)
                                                                                          ggTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcat                                                             >  1:1258186/1‑150 (MQ=255)
                                                                                                cAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcata                                                       >  1:1109544/1‑150 (MQ=255)
                                                                                                  ggTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGCCATAAGTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACg                                                     >  1:336499/1‑150 (MQ=255)
                                                                                                  ggTAACAGAAAAACACCAGACCCCATAATATTCCCCGACACCAGCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATAGGCATCCGAACACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCAAATCATACg                                                     >  2:68640/1‑150 (MQ=255)
                                                                                                   gTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGc                                                    <  2:857604/150‑1 (MQ=255)
                                                                                                    tAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                                   <  2:975866/150‑1 (MQ=255)
                                                                                                     aaCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCt                                                  >  2:1070518/1‑150 (MQ=255)
                                                                                                      aCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTg                                                 <  2:1152764/150‑1 (MQ=255)
                                                                                                      aCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCCTCGTGACAATAAATTACGTCATATCATACGCCTg                                                 >  2:154304/1‑150 (MQ=255)
                                                                                                        aGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGTTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCa                                               <  2:318718/150‑1 (MQ=255)
                                                                                                        aGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCa                                               >  2:773966/1‑150 (MQ=255)
                                                                                                         gAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACACCTGCAt                                              >  1:482279/1‑150 (MQ=255)
                                                                                                          aaaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAtt                                             >  2:991100/1‑150 (MQ=255)
                                                                                                          aaaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAtt                                             >  2:848295/1‑150 (MQ=255)
                                                                                                          aaaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAtt                                             >  1:266176/1‑150 (MQ=255)
                                                                                                              acacCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGAt                                         <  2:634374/150‑1 (MQ=255)
                                                                                                                acCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAg                                       >  1:171272/1‑150 (MQ=255)
                                                                                                                      ccccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGTGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAat                                 <  1:841118/150‑1 (MQ=255)
                                                                                                                       cccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt                                >  1:452854/1‑150 (MQ=255)
                                                                                                                              tatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCCCTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCg                         >  1:50053/1‑150 (MQ=255)
                                                                                                                              tatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACACCTGCATTCGATAGTTGAATTATTATCg                         <  2:482279/150‑1 (MQ=255)
                                                                                                                               atTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGa                        >  2:1034724/1‑150 (MQ=255)
                                                                                                                                ttCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGAc                       >  2:1074229/1‑150 (MQ=255)
                                                                                                                                  ccccGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATtattatatacag                     >  1:75001/1‑144 (MQ=255)
                                                                                                                                   cccGTCACCCTGATGGTGTCGGGGATTTAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCAAACGCCTGCATACGATAGTTGAATTATTATCGACCgg                    >  2:1014594/1‑150 (MQ=255)
                                                                                                                                   cccGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAACCTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATCCGCCTGCATTCGATAGTTGAATTATTATCGACCgg                    >  1:95260/1‑150 (MQ=255)
                                                                                                                                   cccGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCgg                    >  2:409529/1‑150 (MQ=255)
                                                                                                                                   cccGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCgg                    >  2:1010274/1‑150 (MQ=255)
                                                                                                                                   cccGTCACCCTGATGGGGTCGAGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGAAGTCGATAGATGAATTATTATCGACCgg                    >  2:579272/1‑150 (MQ=255)
                                                                                                                                   cccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCgg                    >  1:802282/1‑150 (MQ=255)
                                                                                                                                    ccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCggg                   <  1:304093/150‑1 (MQ=255)
                                                                                                                                     cGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGt                  <  2:494044/150‑1 (MQ=255)
                                                                                                                                     cGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGt                  <  1:425655/150‑1 (MQ=255)
                                                                                                                                           cATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGc            <  1:855842/150‑1 (MQ=255)
                                                                                                                                           cATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGc            <  2:320507/150‑1 (MQ=255)
                                                                                                                                                   tGACGGGGATTAAGCCCACTTTGGGAGCATCAGAATACGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAACTTACGGCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCg    <  1:787908/150‑1 (MQ=255)
                                                                                                                                                     aCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCATGCTAAGGAGAGATTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTTAATTATTATCGACCGGGTTATCGCTGATGCCGtt  <  1:818443/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTT  >  NC_000913/4336839‑4337137

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