Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,339,805 C→T 7.9% Q239Q (CAG→CAA)  adiA ← arginine decarboxylase, degradative

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,339,8050CT 7.9% 188.9 / 5.0 63Q239Q (CAG→CAA) adiAarginine decarboxylase, degradative
Reads supporting (aligned to +/- strand):  ref base C (34/24);  new base T (2/3);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

CGGATAGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCA  >  NC_000913/4339661‑4339953
                                                                                                                                                |                                                                                                                                                    
cGGATAGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCat                                                                                                                                                 >  2:199803/1‑150 (MQ=255)
    tAGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAGGCTTGCatgacc                                                                                                                                             <  2:256408/150‑3 (MQ=255)
    tAGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAGGCTTGCatgacc                                                                                                                                             <  1:1545154/150‑3 (MQ=255)
     aGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAGGCTTGCatgaccg                                                                                                                                            <  1:1213920/150‑4 (MQ=255)
     aGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACATTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATgg                                                                                                                                            <  1:82902/150‑1 (MQ=255)
      gATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGt                                                                                                                                           >  2:221645/1‑150 (MQ=255)
          ggCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCgg                                                                                                                                       >  2:277466/1‑150 (MQ=255)
          ggCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCgg                                                                                                                                       >  1:265952/1‑150 (MQ=255)
          ggCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCgg                                                                                                                                       >  1:1512087/1‑150 (MQ=255)
            cccGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGtt                                                                                                                                     <  2:1015584/150‑1 (MQ=255)
               gATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTaga                                                                                                                                  <  2:972720/150‑1 (MQ=255)
                  aaTGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGcc                                                                                                                               <  2:1512087/150‑1 (MQ=255)
                    tGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCgg                                                                                                                             <  1:952925/150‑1 (MQ=255)
                    tGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGGCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCgg                                                                                                                             <  2:1174350/150‑1 (MQ=255)
                     gCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGa                                                                                                                            >  1:299822/1‑150 (MQ=255)
                            cGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcg                                                                                                                     >  2:83107/1‑150 (MQ=255)
                              gTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAACCCTTGTTCGATGGATTTATGGCAGTTACGGGCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcgac                                                                                                                   <  1:701657/150‑1 (MQ=255)
                               ttGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcgacg                                                                                                                  <  2:1246539/150‑1 (MQ=255)
                                 gCGGCTTGGCACCATGTAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCgacgac                                                                                                                >  1:1269199/1‑150 (MQ=255)
                                  cGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCACACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGCCGACt                                                                                                               >  1:741140/1‑150 (MQ=255)
                                  cGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACt                                                                                                               >  1:569247/1‑150 (MQ=255)
                                   ggCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGCCGACTa                                                                                                              >  1:829402/1‑150 (MQ=255)
                                    gCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTAc                                                                                                             >  1:551989/1‑150 (MQ=255)
                                     cTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTAcc                                                                                                            <  2:1629548/150‑1 (MQ=255)
                                      ttGGCACCATATATACCGGTTTCGCGCCTGTAAGCATCAAACCTTGTTCGATGGATTTATGGCATTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGTATGTACCGACGACTACCg                                                                                                           <  2:1370885/150‑1 (MQ=255)
                                         gCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGAcc                                                                                                        <  1:1392889/150‑1 (MQ=255)
                                           aCCATATAGACCGGGTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAg                                                                                                      <  2:829402/150‑1 (MQ=255)
                                              atatAGACCGGTTTCGCGCCTGTCAGCATCAAACGTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCTCGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAg                                                                                                   <  2:475888/150‑1 (MQ=255)
                                              atatAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAg                                                                                                   <  2:465051/150‑1 (MQ=255)
                                               tataGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGc                                                                                                  >  2:954828/1‑150 (MQ=255)
                                                ataGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCg                                                                                                 >  1:242534/1‑150 (MQ=255)
                                                  aGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGAt                                                                                               >  2:1444182/1‑150 (MQ=255)
                                                   gACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATc                                                                                              >  2:320747/1‑150 (MQ=255)
                                                   gACCGGTTTCGCACCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATc                                                                                              <  1:1487609/150‑1 (MQ=255)
                                                        gTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCAc                                                                                         >  2:880068/1‑150 (MQ=255)
                                                                 cTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATAcgc                                                                                >  1:688371/1‑150 (MQ=255)
                                                                  tGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATAcgcg                                                                               <  1:880068/150‑1 (MQ=255)
                                                                  tGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGAATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATAcgcg                                                                               >  2:1294893/1‑150 (MQ=255)
                                                                   gTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGt                                                                              >  1:72333/1‑150 (MQ=255)
                                                                         aTCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCa                                                                        >  2:500260/1‑150 (MQ=255)
                                                                             aaCCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGCGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATAtt                                                                    >  2:1591608/1‑150 (MQ=255)
                                                                              aCCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATAttt                                                                   >  2:1507239/1‑150 (MQ=255)
                                                                                 ttGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttc                                                                <  2:1325909/150‑1 (MQ=255)
                                                                                   gTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgc                                                              <  1:320747/150‑1 (MQ=255)
                                                                                      cGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgcttt                                                           >  2:416010/1‑150 (MQ=255)
                                                                                        aTGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgctttcg                                                         >  2:643413/1‑150 (MQ=255)
                                                                                                   gCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGcc                                              <  1:199803/150‑1 (MQ=255)
                                                                                                          cGGTCAACGACCACGAAATCGATATCGGTCATGCAAGCCTGCATGATGGTGCGGTAAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATgg                                       >  2:1482023/1‑150 (MQ=255)
                                                                                                               aaCGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATCCGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAg                                  >  2:403642/1‑150 (MQ=255)
                                                                                                                aCGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGc                                 >  2:1281229/1‑150 (MQ=255)
                                                                                                                 cGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCa                                <  2:1269199/150‑1 (MQ=255)
                                                                                                                    ccACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAg                             >  2:796441/1‑150 (MQ=255)
                                                                                                                      aCGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGaa                           >  2:133127/1‑150 (MQ=255)
                                                                                                                       cGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAAc                          >  1:213062/1‑150 (MQ=255)
                                                                                                                            tCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAgg                     >  1:1340741/1‑150 (MQ=255)
                                                                                                                             cGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAggg                    <  1:1444182/150‑1 (MQ=255)
                                                                                                                                  accatCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAgtt               >  1:896180/5‑150 (MQ=255)
                                                                                                                                   ccatCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAgttc              >  1:360957/4‑150 (MQ=255)
                                                                                                                                    ggTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAgttcg             >  1:100517/1‑150 (MQ=255)
                                                                                                                                    ggTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGATTTCGCCAATTGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAgttcg             >  1:1339373/1‑150 (MQ=255)
                                                                                                                                       cATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGGCAAGCAAAGAACCGAGGGAAgttcgttc          <  2:1069341/150‑1 (MQ=255)
                                                                                                                                        aTGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGttcgttcg         <  1:796441/150‑1 (MQ=255)
                                                                                                                                              gCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGccc   <  1:1507239/150‑1 (MQ=255)
                                                                                                                                               ccTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCa  >  2:1238621/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CGGATAGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCCCA  >  NC_000913/4339661‑4339953

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: ‑