Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,322,394 C→G 15.5% 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,3940CG15.5% 393.6 / 62.6 157D185H (GAC→CAC) trpEcomponent I of anthranilate synthase
Reads supporting (aligned to +/- strand):  ref base C (69/62);  new base G (17/7);  total (88/69)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

Only 100 of 162 total aligned reads displayed.

GCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATA  >  NC_000913/1322295‑1322493
                                                                                                   |                                                                                                   
gCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTc                                                                                                     >  2:3488827/1‑100 (MQ=255)
gCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTc                                                                                                     <  2:2869850/100‑1 (MQ=255)
 cTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCa                                                                                                    <  2:1101446/100‑1 (MQ=255)
 cTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCa                                                                                                    >  2:1852866/1‑100 (MQ=255)
  tGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCaa                                                                                                   <  1:125706/100‑1 (MQ=255)
   gACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAt                                                                                                  >  1:469581/1‑100 (MQ=255)
      gTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCAc                                                                                               <  2:2122591/100‑1 (MQ=255)
        aGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCa                                                                                             >  2:3314040/1‑100 (MQ=255)
         gttcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCAt                                                                                            >  1:102715/1‑100 (MQ=255)
         gttcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCAt                                                                                            >  1:3015544/1‑100 (MQ=255)
         gttcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCAt                                                                                            >  2:958206/1‑100 (MQ=255)
         gttcgttcAGGCGAGCAGTGAGACGTTGGTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCa                                                                                             <  2:2441834/99‑1 (MQ=255)
           tcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAATAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCa                                                                                          <  1:1650262/100‑1 (MQ=255)
            tgttcAGTCGAGCAGTGAGACTTTGTTTTTCTTCTTCATTCTGGGCAAACAGTCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAg                                                                                         <  2:1603646/99‑1 (MQ=255)
             gttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGc                                                                                        >  2:2958342/1‑100 (MQ=255)
             gttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGc                                                                                        <  1:901407/100‑1 (MQ=255)
              ttcATGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAga                                                                                        <  2:3522388/99‑2 (MQ=255)
              ttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAgaa                                                                                       <  1:2378835/100‑3 (MQ=255)
              ttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCg                                                                                       <  2:270163/100‑1 (MQ=255)
                cAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGtt                                                                                     >  1:2820437/1‑100 (MQ=255)
                  ggCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTc                                                                                   >  1:2446716/1‑100 (MQ=255)
                    cGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAg                                                                                 >  2:1866971/1‑100 (MQ=255)
                     gAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGc                                                                                >  1:522306/1‑100 (MQ=255)
                        cAGTGAGACGTTGTTTTTCTTCTTCTTTCGGAGCAAACAGGCTGGCCTGAATATGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCgag                                                                             <  2:2025011/100‑1 (MQ=255)
                         attGAGACGTTGTGTTTCTTCTTCATGCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCga                                                                              <  1:3432233/96‑1 (MQ=255)
                         aGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCgaga                                                                            <  1:586569/100‑1 (MQ=255)
                          gTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGAt                                                                           <  2:2484714/100‑1 (MQ=255)
                            gagaCGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATaa                                                                         <  1:3493068/100‑1 (MQ=255)
                            gagaCGTTGTTTTTCTTCTTCATTCAGAGCAAACAGGCTGGCCTGAATACGTGTACTTTTTTTCTGATGGTTAATCACCATCAGCGTTTTAGCCAGCTa                                                                          >  2:1709463/1‑99 (MQ=255)
                             agaCGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGTTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATaaa                                                                        >  2:3234038/1‑100 (MQ=255)
                              gaCGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATaaaa                                                                       >  2:531448/1‑100 (MQ=255)
                                cGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCCGCGAGATaaaatc                                                                     <  1:2470157/100‑3 (MQ=255)
                                cGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAAc                                                                     <  2:3015402/100‑1 (MQ=255)
                                 gttgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACa                                                                    >  1:1127749/1‑100 (MQ=255)
                                 gttgttTTTCTCCTTCATTCGGAGCAAACAGGCTGGTCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCCCCATCAGCGTTTCAGCGAGATAAAAACa                                                                    <  2:2141107/100‑1 (MQ=255)
                                  ttgttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAg                                                                   <  2:246694/100‑1 (MQ=255)
                                     tttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGaaa                                                                >  1:2161301/1‑100 (MQ=255)
                                         tcttcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAg                                                            >  1:1061736/1‑100 (MQ=255)
                                          cttcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAgg                                                           <  1:808853/100‑1 (MQ=255)
                                             cttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATAGTCAATCACCATCAGCTTTTCAGCGAGATAAAAACAGAAATCAggg                                                          >  2:3258458/1‑98 (MQ=255)
                                                 aTTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGt                                                      <  2:375803/98‑1 (MQ=255)
                                                 aTTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTa                                                    >  2:3635150/1‑100 (MQ=255)
                                                 aTTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTa                                                    <  1:1852875/100‑1 (MQ=255)
                                                 aTTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTa                                                    <  1:1852736/100‑1 (MQ=255)
                                                     ggAGCAAACAGGCTGGCCTGAATAGGGGTGCTTTTTTTAGGATGGTCAATCACCATCATTGTTTCAGCGAGATATATACAGAAGTCAGGGCAGTTAttt                                                 <  2:1154826/99‑1 (MQ=255)
                                                     ggAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTCATGGTCAATCACCATCAGCGTTTCAGCGATATAAAAACAGAAATCAGGGCAGTTAtttt                                                >  2:240489/1‑100 (MQ=255)
                                                      gAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCCGTTAtttt                                                >  2:1804148/1‑99 (MQ=255)
                                                       aGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTcc                                              >  1:3564097/1‑100 (MQ=255)
                                                       aGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTcc                                              <  1:958856/100‑1 (MQ=255)
                                                       aGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTcc                                              >  1:1432916/1‑100 (MQ=255)
                                                        gCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCg                                             >  1:2036131/1‑100 (MQ=255)
                                                          aaaCAGGCTGGCCTGAATACGGGTGCTTTTTTTTTTTTTTTCAATCAACATCAGCGTTTCAGCGATATAAAAACAGAAATCAGGGCAGTTAtttt                                                >  2:178545/1‑95 (MQ=255)
                                                          aaaCAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCt                                           >  1:180622/1‑100 (MQ=255)
                                                           aaCAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCt                                           >  2:1128767/1‑99 (MQ=255)
                                                           aaCAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTg                                          >  1:645005/1‑100 (MQ=255)
                                                           aaCAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTg                                          <  1:117477/100‑1 (MQ=255)
                                                               ggCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATACAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAg                                      >  1:1076887/1‑100 (MQ=255)
                                                                gCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGt                                     >  1:3164803/1‑100 (MQ=255)
                                                                 cTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAgct                                    >  1:610019/1‑98 (MQ=255)
                                                                  ctgCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTg                                   <  1:1473329/98‑1 (MQ=255)
                                                                    gCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTgag                                 >  2:1473321/1‑98 (MQ=255)
                                                                    gCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGc                                  <  2:3305135/99‑1 (MQ=255)
                                                                     ccTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCgg                                <  1:648294/100‑1 (MQ=255)
                                                                     ccTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCgg                                <  1:967981/100‑1 (MQ=255)
                                                                        gAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTaa                             >  1:3301723/1‑100 (MQ=255)
                                                                         aaTACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTaaa                            <  2:610017/100‑1 (MQ=255)
                                                                             cGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCtt                        >  1:1867571/1‑100 (MQ=255)
                                                                               ggggCTCTTTGTCTGATGGTCGTTCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATc                          <  2:1676503/93‑1 (MQ=255)
                                                                               ggTTCTTTTTTTCTGATCTTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCtt                        <  2:1202023/98‑1 (MQ=255)
                                                                               ggTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCa                      <  1:47959/100‑1 (MQ=255)
                                                                                gTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCaa                     <  2:1520221/100‑1 (MQ=255)
                                                                                gTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATATTCaa                     >  1:317032/1‑100 (MQ=255)
                                                                                  gctgaaacgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAAt                   >  1:2766770/10‑100 (MQ=37)
                                                                                  gCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAAt                   >  1:2006107/1‑100 (MQ=255)
                                                                                   ctgaaacgcTGATGGTGATTGACCATCAGCGGTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATc                  >  1:3389099/9‑100 (MQ=37)
                                                                                   cTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATc                  >  2:3382530/1‑100 (MQ=255)
                                                                                   cTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATc                  >  1:1450307/1‑100 (MQ=255)
                                                                                    tttttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATcc                 <  2:55113/100‑1 (MQ=255)
                                                                                    tttttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATcc                 >  1:2355660/1‑100 (MQ=255)
                                                                                      tttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCg               <  2:2494468/100‑1 (MQ=255)
                                                                                       aacgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGc              >  1:1271901/5‑100 (MQ=255)
                                                                                        acgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTCACAGTTGCGGTAAATCTTCAAATCCCGcc             >  2:2190983/4‑100 (MQ=255)
                                                                                         cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCg               >  2:3130889/3‑97 (MQ=255)
                                                                                         cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCca            >  2:3546369/3‑100 (MQ=255)
                                                                                         cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCca            >  1:2312037/3‑100 (MQ=255)
                                                                                          tCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac           <  1:1551446/100‑1 (MQ=255)
                                                                                          gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac           >  1:1803078/2‑100 (MQ=255)
                                                                                          gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac           >  1:2954953/2‑100 (MQ=255)
                                                                                          gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac           >  2:2764679/2‑100 (MQ=255)
                                                                                          gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac           >  1:3325370/2‑100 (MQ=255)
                                                                                           cTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcaca          >  2:2941279/1‑100 (MQ=255)
                                                                                             gATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAg        <  1:2266852/100‑1 (MQ=255)
                                                                                              aTGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAgg       <  2:1640305/100‑1 (MQ=255)
                                                                                               tGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGt      >  1:1426998/1‑100 (MQ=255)
                                                                                                ggTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTc     >  1:773484/1‑100 (MQ=255)
                                                                                                 gTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACTAg        >  2:52333/1‑96 (MQ=255)
                                                                                                 gTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCa    >  2:3054917/1‑100 (MQ=255)
                                                                                                 gTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCa    <  2:1309009/100‑1 (MQ=255)
                                                                                                  tCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCAt   >  1:832370/1‑100 (MQ=255)
                                                                                                   cAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATa  <  2:832364/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: ‑