Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 650,532 T→C 5.9% A84A (GCA→GCG citD ← citrate lyase acyl carrier protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,5320TC 5.9% 203.9 / 2.6 69A84A (GCA→GCGcitDcitrate lyase acyl carrier protein
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base C (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACGCTGC  >  NC_000913/650384‑650680
                                                                                                                                                    |                                                                                                                                                    
gCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTg                                                                                                                                                     <  1:664215/150‑1 (MQ=255)
  cGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCAATGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGcc                                                                                                                                                   <  2:557694/150‑1 (MQ=255)
    ggTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcag                                                                                                                                                 >  2:971068/1‑150 (MQ=255)
     gTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagc                                                                                                                                                <  1:1298081/150‑1 (MQ=255)
      tAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagca                                                                                                                                               <  1:465051/150‑1 (MQ=255)
        gATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCACGCGCcagcggt                                                                                                                                             <  1:232618/150‑4 (MQ=255)
        gATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCACGCGCcagcggt                                                                                                                                             <  1:172017/150‑4 (MQ=255)
          tGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGc                                                                                                                                           >  2:777663/1‑150 (MQ=255)
            aaGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCtt                                                                                                                                         >  1:901721/1‑150 (MQ=255)
            aaGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCtt                                                                                                                                         <  1:415555/150‑1 (MQ=255)
            aaGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCtt                                                                                                                                         <  1:1393616/150‑1 (MQ=255)
              cgatctGCTGTCCATCGCGGCATTGGCACCAGGCACAAACAACATGCTTCGGCGGGGGCGTGTTTTACGTTGTTGCAGCGAAGGGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTcc                                                                                                                                       <  1:235539/145‑1 (MQ=255)
                   tGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCg                                                                                                                                  >  1:106238/1‑150 (MQ=255)
                    gCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGt                                                                                                                                 <  1:1304695/150‑1 (MQ=255)
                    gCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGCCAATCAACCCATGGAACAGCACGCATACCGCTGGCCCGTCCCAGCACGGCTTCAAgaacg                                                                                                                                 >  2:235539/1‑146 (MQ=255)
                     cTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTg                                                                                                                                >  1:158526/1‑150 (MQ=255)
                     cTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTg                                                                                                                                <  1:845295/150‑1 (MQ=255)
                       gACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCa                                                                                                                              >  2:496926/1‑150 (MQ=255)
                         ccATCGCGGCATTGGCACCAGGCACAAACAATATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCAc                                                                                                                             <  2:786969/149‑1 (MQ=255)
                         ccATCGCGGCATTGGCACCAGGCACAAACAATATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCAc                                                                                                                             >  1:786969/1‑149 (MQ=255)
                               cgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGATATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAAt                                                                                                                      >  2:1035104/1‑150 (MQ=255)
                               cgGCATGGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCCGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAAt                                                                                                                      <  2:901721/150‑1 (MQ=255)
                                ggCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATg                                                                                                                     <  1:1240394/150‑1 (MQ=255)
                                  cATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCa                                                                                                                   >  1:1118467/1‑150 (MQ=255)
                                    ttGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGt                                                                                                                 >  1:1023877/1‑150 (MQ=255)
                                      ggCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTcc                                                                                                               <  2:1115026/150‑1 (MQ=255)
                                      ggCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACCTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGCATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTCCACGTAAAATGCAGTcc                                                                                                               >  2:1029410/1‑150 (MQ=255)
                                               cacaAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAg                                                                                                                                     <  2:1149235/119‑1 (MQ=255)
                                               cacaAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAg                                                                                                                                     >  1:1149235/1‑119 (MQ=255)
                                               cacaAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCt                                                                                                      >  2:1065192/1‑150 (MQ=255)
                                                  aaacaaCATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTg                                                                                                   >  2:1030602/1‑150 (MQ=255)
                                                   aacaacATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAACTGCAGTCCAGTGCGCCTTTGt                                                                                                  >  1:735824/1‑150 (MQ=255)
                                                   aacaacATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGt                                                                                                  <  2:158526/150‑1 (MQ=255)
                                                       acaTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATAGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATc                                                                                              >  1:1409101/1‑150 (MQ=255)
                                                              tcggcggGTGAGAGTTTTACGTTGTTGCAGCGATGCTGAAATCATTGGCAATCCTCCAATGGCAGAGCCGGGATACCGCTTGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCTGTCCAGTGCGCCTTTGTCATCGACATTc                                                                                       <  2:496886/149‑1 (MQ=255)
                                                                 gcggGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGc                                                                                    <  2:45693/150‑1 (MQ=255)
                                                                  cggGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCt                                                                                   >  1:1161398/1‑150 (MQ=255)
                                                                   gggTGCGAGTTTTACGTTGTTGCCGCGAAGCGGAAATCATTGTCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGGGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCTGTGCGCCTTTGTCATCGACATTCAGCTg                                                                                  <  2:181442/150‑1 (MQ=255)
                                                                           gTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcgccgc                                                                          >  1:497743/1‑150 (MQ=255)
                                                                           gTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAAAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcgccgc                                                                          <  2:358867/150‑1 (MQ=255)
                                                                             tttACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCgcgc                                                                        <  2:453201/150‑1 (MQ=255)
                                                                              ttACGTTGTTGCAGCGAAGCGGAAATCATTTGCAATCCTCCCATGGCACAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCa                                                                       >  2:42878/1‑150 (MQ=255)
                                                                               tACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCAc                                                                      >  1:456430/1‑150 (MQ=255)
                                                                                      ttgCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAg                                                               <  1:1065192/150‑1 (MQ=255)
                                                                                        gCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCg                                                             >  2:324634/1‑150 (MQ=255)
                                                                                            cGAAGCGGAAATCATTGGCAATCCTCCCCTGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCg                                                         <  2:735824/150‑1 (MQ=255)
                                                                                                  ggAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACg                                                   >  2:808758/1‑150 (MQ=255)
                                                                                                  ggAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACg                                                   >  2:1359940/1‑150 (MQ=255)
                                                                                                  ggAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGACAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACg                                                   >  2:865872/1‑150 (MQ=255)
                                                                                                   gAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGt                                                  >  1:880691/1‑150 (MQ=255)
                                                                                                    aaaTCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTc                                                 <  2:984052/150‑1 (MQ=255)
                                                                                                            ggCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAAtgg                                         <  2:382214/150‑1 (MQ=255)
                                                                                                                  cctccCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGc                                                              <  2:470233/123‑1 (MQ=255)
                                                                                                                  cctccCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGc                                                              >  1:470233/1‑123 (MQ=255)
                                                                                                                      ccATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAAt                               >  1:753492/1‑150 (MQ=255)
                                                                                                                          ggCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCTCGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCa                           >  1:463452/1‑150 (MQ=255)
                                                                                                                            cAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATc                         <  2:1161398/150‑1 (MQ=255)
                                                                                                                              gagCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGc                       >  2:486381/1‑150 (MQ=255)
                                                                                                                               agCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGcc                      <  2:437232/150‑1 (MQ=255)
                                                                                                                                  cGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCaaa                   >  1:1217031/1‑150 (MQ=255)
                                                                                                                                   gggATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGTGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAAc                  <  2:1023877/150‑1 (MQ=255)
                                                                                                                                   gggATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAAc                  <  1:777663/150‑1 (MQ=255)
                                                                                                                                        cctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGttt             >  1:356623/4‑150 (MQ=255)
                                                                                                                                        cctgCTGGCACGCGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGttt             >  1:1354383/4‑150 (MQ=255)
                                                                                                                                         ccGCTGGCGCGTGCCAGCAGGGCTTCCGGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTtc            >  1:350050/1‑150 (MQ=255)
                                                                                                                                             tGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCaa        <  2:350050/150‑1 (MQ=255)
                                                                                                                                              ggCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAAc       <  1:496926/150‑1 (MQ=255)
                                                                                                                                               gcgcGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACg      <  2:106238/150‑1 (MQ=255)
                                                                                                                                                   gTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACgctgc  <  2:1118467/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGTTGTAGCGGGCGAGAACGTCCAGAATGGTGGTGCGAATTGCATCGCCAAACTGTTTCTCAACGCTGC  >  NC_000913/650384‑650680

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: