Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,857,264 A→G G160G (GGT→GGC yidK ← putative cotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,857,2640AG99.2% 1179.5 / ‑2.8 383G160G (GGT→GGCyidKputative cotransporter
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (208/170);  total (209/174)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 405 total aligned reads displayed.

GTCCGCCACTGCCATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCACCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATCGGCAGAAAGCAGAC  >  NC_000913/3857174‑3857353
                                                                                          |                                                                                         
gTCCGCCACTGCCATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGcc                                                                                         >  1:4088761/1‑93 (MQ=255)
gTCCGCCACTGCCATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGcc                                                                                         <  2:105254/93‑1 (MQ=255)
gTCCGCCACTGCCATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGcc                                                                                         >  2:2186091/1‑93 (MQ=255)
 tCCGCCACTGCCATTGCGCGCAGTCTGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCg                                                                                        >  1:5947162/1‑93 (MQ=255)
  ccGCCACTGCCATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGt                                                                                       >  1:91002/1‑93 (MQ=255)
   cGCCACTGCCATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTg                                                                                      <  2:5747942/93‑1 (MQ=255)
   cGCCACTGCCATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTg                                                                                      <  1:4497620/93‑1 (MQ=255)
    cccACTGCCATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTga                                                                                     <  2:3354972/92‑1 (MQ=255)
    cccACTGCCATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTga                                                                                     <  1:1506645/92‑1 (MQ=255)
     ccACTGCCATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTgag                                                                                    >  1:5854544/1‑93 (MQ=255)
       aCTGCCATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGaa                                                                                  >  1:83177/1‑93 (MQ=255)
       aCTGCCATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGaa                                                                                  >  2:2515080/1‑93 (MQ=255)
        cTGCCATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGACaaa                                                                                 <  2:3051491/93‑1 (MQ=255)
         tGCCATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAAt                                                                                >  2:5203647/1‑93 (MQ=255)
          gCCATTTCGCGCAGTTCGCAGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATc                                                                               <  1:2213501/93‑1 (MQ=255)
           ccATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCt                                                                              <  1:7848167/93‑1 (MQ=255)
            cATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTg                                                                             >  2:2792672/1‑93 (MQ=255)
             aTTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGt                                                                            >  2:6894720/1‑93 (MQ=255)
             aTTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGt                                                                            <  1:6622995/93‑1 (MQ=255)
             aTTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGt                                                                            >  2:6404524/1‑93 (MQ=255)
              ttGGGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTa                                                                           >  1:3304452/1‑93 (MQ=255)
              ttGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGAGAGCCGCGCCGTGAGAAATCTGTa                                                                           >  2:5662415/1‑93 (MQ=255)
                  gcgcAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGa                                                                       >  2:5452504/1‑93 (MQ=255)
                    gcAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGAtt                                                                     >  2:171617/1‑93 (MQ=255)
                    gcAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGAtt                                                                     <  2:5947085/93‑1 (MQ=255)
                     gaGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTc                                                                    <  1:1622423/92‑1 (MQ=255)
                     cAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTc                                                                    >  2:1792742/1‑93 (MQ=255)
                       gTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTccc                                                                  >  1:2649967/1‑93 (MQ=255)
                        tCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTcccc                                                                 <  1:2477344/93‑1 (MQ=255)
                        tCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTcccc                                                                 >  1:4942795/1‑93 (MQ=255)
                          cgccgATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGa                                                               <  1:5387733/93‑1 (MQ=255)
                            ccgATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACg                                                             >  2:7393817/1‑93 (MQ=255)
                             cgATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGt                                                            <  1:6886095/93‑1 (MQ=255)
                             cgATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGt                                                            >  1:6379032/1‑93 (MQ=255)
                              gATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTg                                                           >  2:6309244/1‑93 (MQ=255)
                              gATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTg                                                           >  2:6800479/1‑93 (MQ=255)
                              gATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTg                                                           <  2:1279264/93‑1 (MQ=255)
                               aTCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGa                                                          >  1:7066189/1‑93 (MQ=255)
                                tCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGaa                                                         >  2:4221203/1‑93 (MQ=255)
                                 cACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGaaa                                                        >  1:57157/1‑93 (MQ=255)
                                 cACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGaaa                                                        <  2:40749/93‑1 (MQ=255)
                                  aCCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAAc                                                       >  2:2571084/1‑93 (MQ=255)
                                   ccGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACa                                                      <  2:3730205/93‑1 (MQ=255)
                                   ccGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACa                                                      <  1:4873587/93‑1 (MQ=255)
                                     gCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAgg                                                    >  2:193726/1‑93 (MQ=255)
                                      cATACAAAATTCCCGCCAGACCAAGCAAAATTCCCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGc                                                   <  2:7810289/93‑1 (MQ=255)
                                      cATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATGGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGc                                                   <  2:4868261/93‑1 (MQ=255)
                                      cATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGc                                                   >  2:7840333/1‑93 (MQ=255)
                                       aTACAAAATTCCCGCCAGATCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCt                                                  >  1:7339874/1‑93 (MQ=255)
                                        tACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGCAATCTGTAGCGATTCCCCGACGTGAAACAGGCTg                                                 >  1:2038933/1‑93 (MQ=255)
                                        tACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTg                                                 >  1:1375776/1‑93 (MQ=255)
                                        tACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTg                                                 >  1:3622213/1‑93 (MQ=255)
                                        tACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTg                                                 >  1:4694716/1‑93 (MQ=255)
                                          cAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGtt                                               <  2:657484/93‑1 (MQ=255)
                                          cAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGtt                                               <  1:599227/93‑1 (MQ=255)
                                           aaaaTTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGctgctg                                              >  2:7043035/1‑93 (MQ=255)
                                             actccccGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAg                                            <  2:6972367/89‑1 (MQ=255)
                                               tccccGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAgcg                                          <  2:4469772/91‑1 (MQ=255)
                                                  ccGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCa                                       <  2:4669594/93‑1 (MQ=255)
                                                   ccccAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCaa                                      >  2:3497453/3‑93 (MQ=255)
                                                   cGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCaa                                      <  2:767614/93‑1 (MQ=255)
                                                    gCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAAc                                     >  1:2843754/1‑93 (MQ=255)
                                                    cccAGCCCAAGCAAAATTACCAGTAGCCAGATAGCCGCTCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAAc                                     <  2:5489400/92‑1 (MQ=255)
                                                     ccAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAAcg                                    >  2:1568302/1‑93 (MQ=255)
                                                      cAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAAcgc                                   >  2:6159856/1‑93 (MQ=255)
                                                       aGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAAcgcg                                  >  1:3723255/1‑93 (MQ=255)
                                                        gACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAAcgcac                                 <  1:7423962/93‑3 (MQ=255)
                                                        gACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAAcgcac                                 <  2:5364361/93‑3 (MQ=255)
                                                        gACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACgcgc                                 >  2:1195989/1‑93 (MQ=255)
                                                        gACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACgcgc                                 >  1:4252757/1‑93 (MQ=255)
                                                        gACCAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACgcgc                                 >  2:3213128/1‑93 (MQ=255)
                                                          ccAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGCGCGCCGACGCGCCg                               >  2:876059/1‑93 (MQ=255)
                                                           cAAGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCgg                              >  1:6749293/1‑93 (MQ=255)
                                                             aGCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAg                            >  1:1267893/1‑93 (MQ=255)
                                                              gCAAAATTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGt                           <  1:2771220/93‑1 (MQ=255)
                                                                aaaaTTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTag                         <  1:1057621/93‑1 (MQ=255)
                                                                aaaaTTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTag                         <  1:6435779/93‑1 (MQ=255)
                                                                aaaaTTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTag                         >  2:4989924/1‑93 (MQ=255)
                                                                aaaaTTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTag                         >  1:3464178/1‑93 (MQ=255)
                                                                 aaaTTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTaga                        <  1:2964057/93‑1 (MQ=255)
                                                                  aaTTACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTagag                       >  2:5222626/1‑93 (MQ=255)
                                                                    ttACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGaa                     >  1:1812231/1‑93 (MQ=255)
                                                                     tACCAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAAc                    <  1:531259/93‑1 (MQ=255)
                                                                       ccAGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACaa                  >  1:7625051/1‑93 (MQ=255)
                                                                         aGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATc                <  1:3070530/93‑1 (MQ=255)
                                                                         aGTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATc                <  2:7095867/93‑1 (MQ=255)
                                                                          gTAGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATCg               >  1:3995677/1‑93 (MQ=255)
                                                                            aGCCAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATCGGc             <  2:5906365/93‑1 (MQ=255)
                                                                               cAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATCGGCAGa          <  1:3186557/93‑1 (MQ=255)
                                                                               cAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATCGGCAGa          <  1:6250535/93‑1 (MQ=255)
                                                                               cAGATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATCGGCAGa          >  2:1452654/1‑93 (MQ=255)
                                                                                 gATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATCGGCAGaaa        >  1:5222415/1‑93 (MQ=255)
                                                                                 gATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATCGGCAGaaa        >  1:696020/1‑93 (MQ=255)
                                                                                 gATAGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGATGAGCGCCAACGCGCCGGAGTAGAGAACAATCGGCAGaaa        >  1:7763918/1‑93 (MQ=255)
                                                                                   tAGCGGCGCCGTGATAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATCGGCAGAAAGc      <  1:6109568/93‑1 (MQ=255)
                                                                                    aGCCGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATCGGCAGAAAGCa     >  2:2363899/1‑93 (MQ=255)
                                                                                      ccGCGCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATCGGCAGAAAGCAGa   <  2:456851/93‑1 (MQ=255)
                                                                                       cgcgCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATCGGCAGAAAGCAGAc  <  1:3954438/93‑1 (MQ=255)
                                                                                       cgcgCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATCGGCAGAAAGCAGAc  >  1:4321680/1‑93 (MQ=255)
                                                                                       cgcgCCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATCGGCAGAAAGCAGAc  <  1:1675742/93‑1 (MQ=255)
                                                                                          |                                                                                         
GTCCGCCACTGCCATTGCGCGCAGTCCGCCGATCACCGCATACAAAATTCCCGCCAGACCAAGCAAAATTACCAGTAGCCAGATAGCCGCACCGTGAGAAATCTGTAGCGATTCCCCGACGTGAAACAGGCTGTTGAGCGCCAACGCGCCGGAGTAGAGAACAATCGGCAGAAAGCAGAC  >  NC_000913/3857174‑3857353

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