Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,322,394 C→G 16.9% D185H (GAC→CAC)  trpE ← component I of anthranilate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,322,3940CG16.9% 467.5 / 89.0 191D185H (GAC→CAC) trpEcomponent I of anthranilate synthase
Reads supporting (aligned to +/- strand):  ref base C (82/71);  new base G (24/7);  total (112/79)
Fisher's exact test for biased strand distribution p-value = 1.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.94e-02

Only 100 of 197 total aligned reads displayed.

GCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATA  >  NC_000913/1322295‑1322493
                                                                                                   |                                                                                                   
gCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTc                                                                                                     <  2:63201/100‑1 (MQ=255)
 cTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCa                                                                                                    >  1:754134/1‑100 (MQ=255)
  tGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCaa                                                                                                   <  2:955097/100‑1 (MQ=255)
    acttAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAt                                                                                                  <  1:1949719/96‑1 (MQ=255)
    aCGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATc                                                                                                 <  1:596669/100‑1 (MQ=255)
     cGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATc                                                                                                 >  2:2402325/1‑99 (MQ=255)
     cGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCa                                                                                                <  2:3959395/100‑1 (MQ=255)
        aGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCa                                                                                             >  2:1240352/1‑100 (MQ=255)
          ttcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATc                                                                                           <  2:2941920/100‑1 (MQ=255)
          ttcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATc                                                                                           <  2:1820180/100‑1 (MQ=255)
            cgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAg                                                                                         <  2:3727539/100‑1 (MQ=255)
            cgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAg                                                                                         <  2:2297107/100‑1 (MQ=255)
             gttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGc                                                                                        <  1:3363235/100‑1 (MQ=255)
              ttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAgaa                                                                                       <  1:1115269/100‑3 (MQ=255)
               tcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAgaaa                                                                                      <  1:3650129/100‑4 (MQ=255)
                  ggCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTc                                                                                   <  1:1364227/100‑1 (MQ=255)
                  ggCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTc                                                                                   >  2:2839665/1‑100 (MQ=255)
                    cGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAgaaaaaaag                                                                                 <  1:4091385/100‑9 (MQ=38)
                      aGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCg                                                                               <  1:2704891/100‑1 (MQ=255)
                        cAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCgag                                                                             <  2:2862382/100‑1 (MQ=255)
                        cAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCg                                                                               >  1:946190/1‑98 (MQ=255)
                         aGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCgaga                                                                            <  1:3177949/100‑1 (MQ=255)
                          gTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTTAATCACCATCAGCGTTTTAGCGAGAt                                                                           >  1:1451642/1‑100 (MQ=255)
                          gTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGAt                                                                           <  2:607060/100‑1 (MQ=255)
                          gTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGAt                                                                           >  2:3145422/1‑100 (MQ=255)
                          gTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAAAAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACGATCAGCGTTTCAGCGAGAt                                                                           >  1:817935/1‑100 (MQ=255)
                           tGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATa                                                                          <  1:4062719/100‑1 (MQ=255)
                                 gttgttTTTCTTCTTCATTCGGATCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAAc                                                                     <  2:2568268/99‑1 (MQ=255)
                                 gttgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAAcn                                                                    <  2:3141383/100‑2 (MQ=255)
                                  ttgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAg                                                                   <  1:3043314/100‑1 (MQ=255)
                                   tgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGa                                                                  >  1:2029871/1‑100 (MQ=255)
                                    gttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGaa                                                                 >  2:1688582/1‑100 (MQ=255)
                                     tttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAACACAGaa                                                                 >  2:325913/1‑99 (MQ=255)
                                     tttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGaaa                                                                <  1:3457735/100‑1 (MQ=255)
                                     tttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGaaa                                                                <  2:1604693/100‑1 (MQ=255)
                                       tttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATc                                                              >  2:1123996/1‑100 (MQ=255)
                                          cttcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAgg                                                           >  2:547837/1‑100 (MQ=255)
                                          cttcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAgg                                                           <  2:701124/100‑1 (MQ=255)
                                          cttcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAgg                                                           >  1:1575078/1‑100 (MQ=255)
                                           ttcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAggg                                                          <  2:3600068/100‑1 (MQ=255)
                                             cgtcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCa                                                        >  2:1913799/3‑100 (MQ=255)
                                              ttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAg                                                       >  1:1868589/1‑100 (MQ=255)
                                              ttcATTCGGACCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTCTCAGCGAGATAAAAACAGAAATCAGGGCAg                                                       <  1:2888394/100‑1 (MQ=255)
                                               tcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGt                                                      >  1:2988493/1‑100 (MQ=255)
                                               tcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGt                                                      >  2:3582293/1‑100 (MQ=255)
                                                 aTTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTa                                                    >  1:3139707/1‑100 (MQ=255)
                                                   tCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTTTTTTTTTCAATCACCATAAGCGTTTCAGCGACCTAAAAACAGAAATCAGGGCAGTTAtt                                                  >  1:1683646/1‑100 (MQ=255)
                                                      gAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTTTTTTTTTTTATTACCAACAGCGGTTTAACGGGATTACGAAAGAAATGAGGGCAGGTAttt                                                 >  1:2556566/1‑98 (MQ=255)
                                                       aGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCCGGGCCGATATTTTc                                               >  2:4118286/1‑99 (MQ=255)
                                                       aGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTcc                                              >  2:4001113/1‑100 (MQ=255)
                                                        gCAAGCAGGCTGGCTCGCGTACGTCTGCATGTGTTCAGATGGCAATTCGCCATCAGCGTGACAGCGAGATAAAAACAGAAATCAGGGCAGTTAtactcc                                              <  1:1978755/99‑6 (MQ=255)
                                                         cAAACAGGCTGGCCTTAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGc                                            <  2:4329044/100‑1 (MQ=255)
                                                         cAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGc                                            <  2:4189531/100‑1 (MQ=255)
                                                         cAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGc                                            >  2:945125/1‑100 (MQ=255)
                                                           aaCAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTg                                          >  2:2716915/1‑100 (MQ=255)
                                                             cAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGAc                                        <  1:4341122/100‑1 (MQ=255)
                                                               ggCTGGCCTGAATACGGGTGCTTTTTTTTTTTTTTTTTATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACa                                       >  1:4122161/1‑99 (MQ=255)
                                                               ggCTGGCCTGAATACGGGTGCTTTTTTTTTTTTTGTTTATCACCATCAGCGTTTCAGCGAGATAAATACAGAAATCAGGGCAGTTATTTTCCGCTGACAg                                      >  1:402401/1‑100 (MQ=255)
                                                               ggCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCCTCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGAc                                        <  2:1586039/98‑1 (MQ=255)
                                                                gCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGAc                                        <  1:3566257/97‑1 (MQ=255)
                                                                   ggCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGc                                  >  1:232679/1‑100 (MQ=255)
                                                                     ccTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCgg                                >  2:716805/1‑100 (MQ=255)
                                                                      cTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGt                               >  1:660759/1‑100 (MQ=255)
                                                                            aCGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCt                         <  2:3352400/100‑1 (MQ=255)
                                                                             cGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCtt                        >  1:2875388/1‑100 (MQ=255)
                                                                               ggTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCa                      <  2:900902/100‑1 (MQ=255)
                                                                                gTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCaa                     <  1:971298/100‑1 (MQ=255)
                                                                                 tGCTTTTTTTCTGATAGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCaaa                    <  2:1176010/100‑1 (MQ=255)
                                                                                  gctgaaacgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAAt                   >  2:2202813/10‑100 (MQ=37)
                                                                                   cTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGCGATAAAAACAGAAATCAGCGCAGTTATATTCCGCTCACAGTTGCGGTCAATCTTCAAATc                  >  1:3078681/1‑100 (MQ=255)
                                                                                   cTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATc                  >  2:456430/1‑100 (MQ=255)
                                                                                    tttttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAAttc                 >  1:2429331/1‑98 (MQ=255)
                                                                                    tttttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATcc                 >  2:2087032/1‑100 (MQ=255)
                                                                                    tttttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATcc                 >  1:3171442/1‑100 (MQ=255)
                                                                                     ttttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATccc                <  2:1677569/100‑1 (MQ=255)
                                                                                      tttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCg               <  1:4197643/100‑1 (MQ=255)
                                                                                        acgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGCTAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCAGcc             >  1:2252503/4‑100 (MQ=255)
                                                                                         ttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGCAATCAGGGCCGTTATTATCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCca            >  1:1929975/1‑100 (MQ=255)
                                                                                         ttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCca            >  2:2616257/1‑100 (MQ=255)
                                                                                         cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGAGAGTTGCGGTACCCATTCAAATCCCg               >  1:1471534/3‑97 (MQ=255)
                                                                                         cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGcc             >  1:132572/3‑99 (MQ=255)
                                                                                         cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCca            >  2:3407708/3‑100 (MQ=255)
                                                                                         cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCca            >  2:3523595/3‑100 (MQ=255)
                                                                                          gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac           >  2:3030250/2‑100 (MQ=255)
                                                                                          gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac           >  1:3411988/2‑100 (MQ=255)
                                                                                          gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac           >  2:3709192/2‑100 (MQ=255)
                                                                                          gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac           >  1:1112344/2‑100 (MQ=255)
                                                                                          gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac           >  2:2585797/2‑100 (MQ=255)
                                                                                          gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATATTCAAATCCCGCcac           >  2:4341747/2‑100 (MQ=255)
                                                                                           cTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcaca          <  1:666346/100‑1 (MQ=255)
                                                                                           cTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcaca          <  2:1520969/100‑1 (MQ=255)
                                                                                           cTGATGGACAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCGGACAGTTGCGGTAGATCTTCAAATCCAGCcac           <  1:239940/99‑1 (MQ=255)
                                                                                            tGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACaa         >  1:851321/1‑100 (MQ=255)
                                                                                              aTGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAgg       <  2:3422415/100‑1 (MQ=255)
                                                                                              aTGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAgg       <  2:1879976/100‑1 (MQ=255)
                                                                                               tGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGt      <  2:3437848/100‑1 (MQ=255)
                                                                                               tGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGt      >  2:174950/1‑100 (MQ=255)
                                                                                                  tCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCAt   <  1:716803/100‑1 (MQ=255)
                                                                                                   cAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAACTCTTCAAATCCCGCCACAAGGTCATa  <  2:2675533/100‑1 (MQ=255)
                                                                                                   cAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATa  <  2:2029879/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATA  >  NC_000913/1322295‑1322493

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: