Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,275,682 (A)5→4 100% coding (401/1572 nt) garD → D‑galactarate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,275,6780A.74.4% 324.7 / 183.5 160coding (397/1572 nt)garDD‑galactarate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (18/22);  new base . (53/66);  total (72/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

CACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCG  >  NC_000913/3275579‑3275776
                                                                                                   |                                                                                                  
cacacGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCa                                                                                                    >  1:1283402/1‑100 (MQ=255)
  cacGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaaa                                                                                                  >  1:1624049/1‑100 (MQ=255)
  cacGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaaa                                                                                                  >  2:279371/1‑100 (MQ=255)
   acGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaaaa                                                                                                 <  1:3974977/100‑1 (MQ=255)
    cGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaaaac                                                                                                <  1:3152309/100‑2 (MQ=255)
     gCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaaaacc                                                                                               <  1:2699102/100‑3 (MQ=255)
     gCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaaaacc                                                                                               <  2:4052852/100‑3 (MQ=255)
     gCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaaaacc                                                                                               <  1:3406848/100‑3 (MQ=255)
      cTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAAcc                                                                                              <  1:370149/100‑1 (MQ=255)
       tGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaaaacctg                                                                                             <  1:1062457/100‑5 (MQ=255)
       tGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaaaacctg                                                                                             <  1:3862275/100‑5 (MQ=255)
        gCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGc                                                                                           >  2:180925/1‑100 (MQ=255)
        gCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGc                                                                                           >  2:1164098/1‑100 (MQ=255)
         ccACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCt                                                                                          >  1:1264518/1‑100 (MQ=255)
         ccACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCt                                                                                          <  1:2521322/100‑1 (MQ=255)
         ccACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGCATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCt                                                                                          <  2:2666237/100‑1 (MQ=255)
           aCTGGCAACCAAAGTCCCGGAACGCTTAGCGGCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCATCGTGGGCACCAAAACCTGCTc                                                                                         <  2:641881/99‑1 (MQ=255)
           aCTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTc                                                                                         <  1:2496965/100‑1 (MQ=255)
            cTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCACTGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCg                                                                                        >  2:2962574/1‑100 (MQ=255)
             tGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCg                                                                                        <  2:632308/98‑1 (MQ=255)
              ggCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTa                                                                                     <  2:992502/100‑1 (MQ=255)
               gCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTAt                                                                                    >  1:1278888/1‑100 (MQ=255)
                cAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATc                                                                                   >  2:1312059/1‑100 (MQ=255)
                 aaCCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATc                                                                                   >  1:3187494/1‑100 (MQ=255)
                 aaCCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATca                                                                                  <  2:1943915/100‑1 (MQ=255)
                 aaCCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATca                                                                                  <  2:3344325/100‑1 (MQ=255)
                 aaCCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATca                                                                                  >  1:1575941/1‑100 (MQ=255)
                  actaaaGTCCCTGAACCCTTACCGCCTCTGGAAGGATACACCTTTGACGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATc                                                                                   <  2:1979237/96‑1 (MQ=255)
                   ctaaaGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATcac                                                                                 <  1:2148220/98‑1 (MQ=255)
                   ccAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATGCACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATcac                                                                                 <  1:2931014/100‑1 (MQ=255)
                   ccAAAGGCCCGGAACCCTTACCGCCACTGGAAGGATACACCTTTGAGCGCGATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATcacc                                                                                <  2:1303478/100‑1 (MQ=255)
                    cAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGTAGCGTGGGCACCAAAACCTGCTCGGTATcacca                                                                               >  2:4121622/1‑100 (MQ=255)
                      aaGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCaccacc                                                                             <  2:1676317/100‑1 (MQ=255)
                      aaGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCaccacc                                                                             <  2:1908330/100‑1 (MQ=255)
                      aaGTCCCGGAACCCTGACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCaccacc                                                                             <  1:1864499/100‑1 (MQ=255)
                        gTCCGGGAGCCCTTATCGCGGCTGGAAGGGGACAGGGGTGAGGGCTAGCGCAATGCCGATGGCAGCGTGGGGACCTAAACCTGCTCGGTATgactagc                                                                             <  2:3651795/98‑5 (MQ=255)
                        gTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAg                                                                           <  2:3884305/100‑1 (MQ=255)
                         tCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAg                                                                           <  2:668607/100‑1 (MQ=255)
                          cccGGCACGCTTACCGCCGCTATAAGGATACACCTTTGATGGCTATCGCAATGCCGATGACAGCGTCAGCACCAAAAAACTGCGCGGTATCAccacca                                                                            >  2:1576393/1‑98 (MQ=255)
                            cGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTc                                                                       <  1:1212539/100‑1 (MQ=255)
                             ggAACCCTTACGGCCGCTGGAAGGATACACCTTTGAGGGCTCTCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTcc                                                                      <  2:1187721/100‑1 (MQ=255)
                             ggAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTcc                                                                      <  1:1820586/100‑1 (MQ=255)
                              gAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTcc                                                                      <  1:595587/100‑1 (MQ=255)
                                aCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCAc                                                                    <  1:3275482/100‑1 (MQ=255)
                                 cccTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACtg                                                                  >  1:1232161/1‑100 (MQ=255)
                                  ccTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACtg                                                                  >  1:442699/1‑100 (MQ=255)
                                  ccTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACGACCAGCGTCCACtg                                                                  >  1:1554544/1‑99 (MQ=255)
                                  ccTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACtgt                                                                 >  1:346092/1‑100 (MQ=255)
                                   cTTACCCCCGCTGGAAGGATACACCTTTGAGGCCTATAGCAATGCCGGTGGCAGCGTGGGCCCGCAAAACTGCTCCGTATCACCACCAGCGTCCAc                                                                    >  2:3667128/1‑96 (MQ=255)
                                     tACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCTCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACtgtgtg                                                              <  2:1606313/100‑1 (MQ=255)
                                     tACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAACACCTGCTCGGTATCACCACCAGCGTCCACtgtgcg                                                              >  2:3838110/1‑98 (MQ=255)
                                     tACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACtgtgtg                                                              >  2:374218/1‑100 (MQ=255)
                                     tACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACtgtgtg                                                              <  1:1212250/100‑1 (MQ=255)
                                      aCCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTgg                                                             >  2:3545707/1‑100 (MQ=255)
                                       ccgccgCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGc                                                            <  1:2636876/100‑1 (MQ=255)
                                           cgcTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAgg                                                         >  1:2017755/1‑100 (MQ=255)
                                            ggaggAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCg                                                       <  1:379804/97‑1 (MQ=255)
                                             cTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCg                                                       >  2:3516823/1‑100 (MQ=255)
                                               ggAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCgtgg                                                    >  2:4021514/1‑100 (MQ=255)
                                                  aGGATGCCCCGTTGAGTGCTAGCGGCATGGCGATCGCACCTTGGGCACCAAAACGTGCTCGGTATCACCGCCAGGGTCCACTGTGTGGCAgg                                                         <  2:320057/92‑1 (MQ=255)
                                                     aTACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGTAGGCGTGGTGGAc                                               <  2:2582964/100‑1 (MQ=255)
                                                     aTACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGAc                                               <  1:1314067/100‑1 (MQ=255)
                                                     aTACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACt                                              >  2:3717557/1‑100 (MQ=255)
                                                         acCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATgt                                          <  2:518582/100‑1 (MQ=255)
                                                         acCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATgt                                          <  1:2903529/100‑1 (MQ=255)
                                                         acCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATgt                                          <  1:1358856/100‑1 (MQ=255)
                                                           cTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACTTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATgtag                                        <  1:2461826/100‑1 (MQ=255)
                                                               gAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTaa                                     <  1:3039392/100‑1 (MQ=255)
                                                               gAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTaaa                                    >  1:1182660/1‑100 (MQ=255)
                                                                   gCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAAt                                 >  1:1038267/1‑100 (MQ=255)
                                                                    cTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCa                               >  2:3675607/1‑100 (MQ=255)
                                                                    cTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCa                               >  1:1862128/1‑100 (MQ=255)
                                                                    cTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCa                               >  1:167034/1‑100 (MQ=255)
                                                                     tATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCAt                              <  1:1753269/100‑1 (MQ=255)
                                                                      aTCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCAtt                             >  2:1303018/1‑100 (MQ=255)
                                                                      aTCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCAtt                             <  1:1658498/100‑1 (MQ=255)
                                                                      aTCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACGGTGTGGCAGGCGTGGTGGGCTATGTAGTAAAAATCAtt                             >  2:3095814/1‑100 (MQ=255)
                                                                       tCGCAATGGCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCAACTGTGTGGGAGGCGTGGTGGACTATGTAGTAAAAATCATTg                            <  2:1038034/100‑1 (MQ=255)
                                                                       tCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTg                            >  1:88312/1‑100 (MQ=255)
                                                                       tCGCAATGCCGATGGCAGCATGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAGAATCAttc                            <  2:4061447/100‑2 (MQ=255)
                                                                        cGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGa                           >  2:2787326/1‑100 (MQ=255)
                                                                         gCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGaa                          >  1:152273/1‑100 (MQ=255)
                                                                         gCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGa                           >  1:374584/1‑99 (MQ=255)
                                                                          cAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTg                            <  2:2527571/98‑1 (MQ=255)
                                                                          cAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGaa                          <  2:2858673/100‑1 (MQ=255)
                                                                          cAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAAc                         >  1:860739/1‑100 (MQ=255)
                                                                           aaTGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAAcg                        <  1:2263474/100‑1 (MQ=255)
                                                                              gCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGa                     <  1:3616944/100‑1 (MQ=255)
                                                                                 gATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCt                  <  2:623364/100‑1 (MQ=255)
                                                                                    ggCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTgn                <  2:307904/100‑2 (MQ=255)
                                                                                           tGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGaaa        <  1:51799/100‑1 (MQ=255)
                                                                                            gggCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGaaa        >  2:2709851/1‑100 (MQ=255)
                                                                                             ggCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAAt       >  1:1043113/1‑100 (MQ=255)
                                                                                             ggCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATa      <  2:3903725/100‑1 (MQ=255)
                                                                                             ggCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATa      >  2:1933932/1‑100 (MQ=255)
                                                                                             ggCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATa      >  1:2547530/1‑100 (MQ=255)
                                                                                              gCACCAAAACCTGCTCGGTATCACTACGAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATTTAGTAAAAATCATTGAACGCGATCTGGTCCCGACAt       >  2:3417741/1‑98 (MQ=255)
                                                                                              gCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATAc     <  2:3132009/100‑1 (MQ=255)
                                                                                               cACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATAc     <  1:927477/100‑1 (MQ=255)
                                                                                                 ccAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATAccc   >  2:1376828/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
CACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCG  >  NC_000913/3275579‑3275776

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