Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,322,394 C→G 16.6% 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.6% 228.8 / 37.5 91D185H (GAC→CAC) trpEcomponent I of anthranilate synthase
Reads supporting (aligned to +/- strand):  ref base C (37/37);  new base G (10/5);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

GCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATA  >  NC_000913/1322295‑1322493
                                                                                                   |                                                                                                   
gCTGACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTc                                                                                                     <  1:2081107/100‑1 (MQ=255)
   gACGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAAt                                                                                                  <  2:716216/100‑1 (MQ=255)
     cGTAGTTCTTTCAGGCGAGCAGTGAGACGTTGTTTTGCTTCTGCTTTCGGAGCACACAGGCTGGCCTGAATACTGGTGCTTTTTTTCTGCTGGTCAATCa                                                                                                <  2:198606/100‑1 (MQ=255)
     cGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCa                                                                                                >  2:1422952/1‑100 (MQ=255)
     cGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCa                                                                                                <  1:930619/100‑1 (MQ=255)
     cGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCa                                                                                                >  1:727187/1‑100 (MQ=255)
     cGTAGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCa                                                                                                <  2:2029414/100‑1 (MQ=255)
        aGTTCGTTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCa                                                                                             <  1:951337/100‑1 (MQ=255)
         gttcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGTGTGCTTTTTTTCTGATGGTCAATCACCAt                                                                                            >  2:391101/1‑100 (MQ=255)
          ttcgttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATc                                                                                           <  1:2176133/100‑1 (MQ=255)
          ttcggtcAGGCGAGCAGTGAGGCGTTGTTTTTCTTCTTGATTCGGAGCGAACAGGCTGGCCTGAATACGGTTGCTTTTTTTCTGATGGTCAATCACCATc                                                                                           <  2:268045/100‑1 (MQ=255)
             gttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAga                                                                                        <  1:74320/100‑2 (MQ=255)
             gttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGc                                                                                        >  1:1280020/1‑100 (MQ=255)
             gttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGc                                                                                        >  2:337354/1‑100 (MQ=255)
              ttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTTGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAgaa                                                                                       <  2:388522/100‑3 (MQ=255)
              ttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCg                                                                                       <  2:522728/100‑1 (MQ=255)
               tcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAgaaa                                                                                      <  1:1037976/100‑4 (MQ=255)
               tcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAgaaa                                                                                      <  1:1640586/100‑4 (MQ=255)
                cAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTGATTGACCATCAgaaaa                                                                                     <  2:1360796/100‑5 (MQ=255)
                  ggCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCGAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTc                                                                                   <  2:2479182/100‑1 (MQ=255)
                  ggCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTc                                                                                   <  2:765434/100‑1 (MQ=255)
                  ggCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTc                                                                                   >  2:228836/1‑100 (MQ=255)
                     gAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGc                                                                                >  1:1340291/1‑100 (MQ=255)
                      aGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCg                                                                               >  2:1465839/1‑100 (MQ=255)
                      aGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCg                                                                               >  1:2333125/1‑100 (MQ=255)
                        cAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCgag                                                                             >  2:2021377/1‑100 (MQ=255)
                         aGTGGGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCgaga                                                                            <  2:266099/100‑1 (MQ=255)
                          gTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATCCGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCgaga                                                                            >  2:839684/1‑99 (MQ=255)
                          gTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGAt                                                                           >  1:2341613/1‑100 (MQ=255)
                          gTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGAt                                                                           <  1:1760859/100‑1 (MQ=255)
                          gTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGAt                                                                           <  1:1760658/100‑1 (MQ=255)
                          gTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGCCAATCACCATCAGCGTTTCAGCGAGAt                                                                           >  1:700760/1‑100 (MQ=255)
                             ataCGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATaaa                                                                        <  1:865059/98‑1 (MQ=255)
                                cGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAAc                                                                     <  2:486237/100‑1 (MQ=255)
                                cGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAAc                                                                     >  1:2372382/1‑100 (MQ=255)
                                  atctttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAg                                                                   <  1:1325318/97‑1 (MQ=255)
                                  atctttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAg                                                                   <  1:1325401/97‑1 (MQ=255)
                                    gttTTTCTTCTTTATTCGGAGCAAACAGGCTGGCCTTAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATTAGCGTTTGAGCGAGATAAAAACata                                                                  >  2:2322842/1‑97 (MQ=255)
                                    gttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGaa                                                                 >  2:13772/1‑100 (MQ=255)
                                    gttTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGaa                                                                 <  2:2354556/100‑1 (MQ=255)
                                      ttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGagat                                                               >  2:1325392/1‑97 (MQ=255)
                                      ttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGagat                                                               >  2:1325309/1‑97 (MQ=255)
                                      ttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAAt                                                               >  1:2281349/1‑100 (MQ=255)
                                        ttcttcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCa                                                             >  1:1286824/1‑100 (MQ=255)
                                                      gAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTc                                               >  1:320815/1‑100 (MQ=255)
                                                        gCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCg                                             >  2:2110261/1‑100 (MQ=255)
                                                          aaaCAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCt                                           <  1:2252359/100‑1 (MQ=255)
                                                           aaCAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTg                                          <  1:544731/100‑1 (MQ=255)
                                                           aaCAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTg                                          <  2:156538/100‑1 (MQ=255)
                                                                 cTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGtt                                    >  1:1316684/1‑100 (MQ=255)
                                                                    gCCTGAATACGGGTGCTTTTTTTCTTATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCg                                 <  2:2407093/100‑1 (MQ=255)
                                                                     ccTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCgg                                <  1:791304/100‑1 (MQ=255)
                                                                     ccTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCgg                                <  2:2274690/100‑1 (MQ=255)
                                                                       tGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGTCAGCTGCGGTa                              <  2:1347157/100‑1 (MQ=255)
                                                                       tGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTa                              >  2:377701/1‑100 (MQ=255)
                                                                       tGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTa                              >  1:1174849/1‑100 (MQ=255)
                                                                       tGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTa                              >  1:1627354/1‑100 (MQ=255)
                                                                          aTACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAAt                           <  1:1615692/100‑1 (MQ=255)
                                                                              gggTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTc                       >  1:1758773/1‑100 (MQ=255)
                                                                               ggTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCa                      <  1:341354/100‑1 (MQ=255)
                                                                               ggTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTCGCGGTAAATCTTCa                      <  1:2110275/100‑1 (MQ=255)
                                                                                gTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCaa                     <  1:2477370/100‑1 (MQ=255)
                                                                                gTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCaa                     >  2:800174/1‑100 (MQ=255)
                                                                                 tGCTTTTTTTCTGATGGTAAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCaaa                    <  2:875430/100‑1 (MQ=255)
                                                                                  gCTTTTTTTCTGATGGTAAAAAACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAAt                   >  2:922946/1‑100 (MQ=255)
                                                                                   cTTTTTTTCTGATGGTGAATAACCATGAGGTTTTCAGCGAGATAAAAACAGAAATCATGGCAGTTATTTTACGCTGACAGTTGCGGTAAAtttaaa                      >  2:1627134/1‑91 (MQ=255)
                                                                                    tttttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATcc                 <  1:2508750/100‑1 (MQ=255)
                                                                                    tttttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATcc                 >  1:688558/1‑100 (MQ=255)
                                                                                     ttttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATccc                >  2:2132694/1‑100 (MQ=255)
                                                                                     gaaacgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAAt                   >  2:1620074/7‑97 (MQ=38)
                                                                                      tttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCg               <  2:1547995/100‑1 (MQ=255)
                                                                                       aacgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCGGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGc              >  2:1798727/5‑100 (MQ=255)
                                                                                        tttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGcc             <  2:1772798/100‑1 (MQ=255)
                                                                                        tttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGcc             <  1:1111745/100‑1 (MQ=255)
                                                                                        acgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGcc             >  1:2101461/4‑100 (MQ=255)
                                                                                         cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCGGCTGACAGTTGCGGCAAATCTTCAAATCCCGcc             >  1:622940/3‑99 (MQ=255)
                                                                                         cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGcc             >  2:1379218/3‑99 (MQ=255)
                                                                                         cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCca            >  1:1285866/3‑100 (MQ=255)
                                                                                         cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCca            >  1:619964/3‑100 (MQ=255)
                                                                                         cgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCca            >  1:236704/3‑100 (MQ=255)
                                                                                          tCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac           <  2:784761/100‑1 (MQ=255)
                                                                                          tCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac           <  1:2517755/100‑1 (MQ=255)
                                                                                          gcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac           >  2:649062/2‑100 (MQ=255)
                                                                                           cTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcaca          >  2:2367907/1‑100 (MQ=255)
                                                                                            tGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACaa         >  2:1159938/1‑100 (MQ=255)
                                                                                              aTGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAgg       <  1:1101109/100‑1 (MQ=255)
                                                                                               tGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGt      >  2:484251/1‑100 (MQ=255)
                                                                                               tGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGt      <  1:228838/100‑1 (MQ=255)
                                                                                               tGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAAACAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATATTCAAATCCCGCCACAAgg       >  1:976317/1‑99 (MQ=255)
                                                                                                ggTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTc     >  1:1908517/1‑100 (MQ=255)
                                                                                                 gTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCa    >  1:350258/1‑100 (MQ=255)
                                                                                                  tCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCAt   >  1:1497963/1‑100 (MQ=255)
                                                                                                   cAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATa  <  1:2322856/100‑1 (MQ=255)
                                                                                                   cAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATa  >  1:2166886/1‑100 (MQ=255)
                                                                                                   cAATCACCAGCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTGTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATa  <  2:2160091/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: