Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,336,974 A→T 8.3% S20T (TCG→ACG)  adiC ← arginine:agmatine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,336,9740AT 8.3% 245.2 / 8.9 101S20T (TCG→ACG) adiCarginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base A (44/46);  new base T (6/2);  total (53/48)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 103 total aligned reads displayed.

CTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGT  >  NC_000913/4336826‑4337121
                                                                                                                                                    |                                                                                                                                                   
cTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGac                                                                                                                                                    <  2:1009742/150‑1 (MQ=255)
  tGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGacac                                                                                                                                                  <  2:340521/150‑1 (MQ=255)
     gTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCAt                                                                                                                                               >  1:320507/1‑150 (MQ=255)
       cGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCa                                                                                                                                             >  2:1011748/1‑150 (MQ=255)
        gAGGAACGACGTTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCAACAACCATCCATAAATGGCAATCCCGCCAGAAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCAGACACCATCAg                                                                                                                                            >  1:857604/1‑150 (MQ=255)
        gAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAg                                                                                                                                            <  1:1192153/150‑1 (MQ=255)
          tgAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCTCCAACCATCCATAAATGTCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAATACACCTGACCCCATAATATTCCCCGACACCATCAggg                                                                                                                                          <  2:1224539/149‑1 (MQ=255)
             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)
                          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:464662/150‑1 (MQ=255)
                                     agagCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGag                                                                                                               <  1:552144/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)
                                                            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                             >  1:266176/1‑150 (MQ=255)
                                                                                                                       aaaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAtt                             >  2:848295/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  <  1:425655/150‑1 (MQ=255)
                                                                                                                                                  cGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGt  <  2:494044/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGT  >  NC_000913/4336826‑4337121

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: