Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,850,750 A→G intergenic (‑135/‑793) ivbL ← / → emrD ilvB operon leader peptide/2‑module integral membrane pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,850,7500AG96.6% 622.6 / 1.2 207intergenic (‑135/‑793)ivbL/emrDilvB operon leader peptide/2‑module integral membrane pump
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base G (99/100);  total (103/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 220 total aligned reads displayed.

CAGCTTTACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACGTGTTCC  >  NC_000913/3850660‑3850839
                                                                                          |                                                                                         
caGCTTTACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGcc                                                                                         <  2:1676993/93‑1 (MQ=255)
caGCTTTACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGcc                                                                                         <  1:2727195/93‑1 (MQ=255)
 aGCTTTACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGccc                                                                                        >  1:1936547/1‑93 (MQ=255)
  gCTTTACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGcccc                                                                                       <  1:4345016/93‑1 (MQ=255)
   cTTTACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCccg                                                                                      >  1:4070741/1‑93 (MQ=255)
   cTTTACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCccg                                                                                      >  2:4081878/1‑93 (MQ=255)
    tttACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCccgc                                                                                     >  1:4219850/1‑93 (MQ=255)
    tttACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCccgc                                                                                     >  1:569858/1‑93 (MQ=255)
      tACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCccgccg                                                                                   >  1:929645/1‑93 (MQ=255)
        cAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGa                                                                                 >  1:1621052/1‑93 (MQ=255)
        cAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGa                                                                                 <  2:134871/93‑1 (MQ=255)
         aGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAg                                                                                >  2:79832/1‑93 (MQ=255)
         aGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAg                                                                                >  1:4419744/1‑93 (MQ=255)
          ggggAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATTACAAAAAGCCCCGCCGGAgg                                                                               >  2:47892/1‑92 (MQ=255)
          ggggAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGc                                                                               <  1:3534431/93‑1 (MQ=255)
              agaCAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAgg                                                                           >  2:3277093/1‑93 (MQ=255)
               gaCAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGt                                                                          >  1:4819/1‑93 (MQ=255)
               gaCAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGt                                                                          >  2:2659598/1‑93 (MQ=255)
               gaCAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGt                                                                          >  2:4115189/1‑93 (MQ=255)
               gaCAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGt                                                                          <  2:3597357/93‑1 (MQ=255)
                aCAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGtt                                                                         >  2:2090903/1‑93 (MQ=255)
                  aaTGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTc                                                                       <  1:3654757/93‑1 (MQ=255)
                  aaTGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTc                                                                       >  2:1238988/1‑93 (MQ=255)
                   aTGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCg                                                                      >  2:3120810/1‑93 (MQ=255)
                   aTGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCg                                                                      >  1:1009575/1‑93 (MQ=255)
                    tGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGt                                                                     >  2:4652348/1‑93 (MQ=255)
                     ggAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTc                                                                    <  2:1312422/93‑1 (MQ=255)
                      gAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCa                                                                   <  2:2694213/93‑1 (MQ=255)
                        aaaTTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGt                                                                 <  1:83191/93‑1 (MQ=255)
                        aaaTTTTTCAGCAAGGGAAAATTAAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGt                                                                 <  2:3227335/93‑1 (MQ=255)
                         aaTTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTc                                                                <  2:4593559/93‑1 (MQ=255)
                           tttttCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGc                                                              <  1:4028601/93‑1 (MQ=255)
                              ttCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTg                                                           >  1:685088/1‑93 (MQ=255)
                               tCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGc                                                          <  1:3264187/93‑1 (MQ=255)
                               tCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGc                                                          <  2:3420049/93‑1 (MQ=255)
                                cAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCg                                                         >  1:594564/1‑93 (MQ=255)
                                  gCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGc                                                       >  2:2229624/1‑93 (MQ=255)
                                    aaGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTg                                                     <  1:1058000/93‑1 (MQ=255)
                                     aGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTgg                                                    >  1:3809327/1‑93 (MQ=255)
                                     aGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTgg                                                    <  1:472455/93‑1 (MQ=255)
                                     aGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTgg                                                    >  2:3327731/1‑93 (MQ=255)
                                     aGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTgg                                                    <  1:2832869/93‑1 (MQ=255)
                                     aGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGGTTCGTCAGTCGCCTGCGGCTgg                                                    >  2:2278636/1‑93 (MQ=255)
                                      gggAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGt                                                   >  1:1708234/1‑93 (MQ=255)
                                      gggAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGt                                                   >  1:241900/1‑93 (MQ=255)
                                      gggAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGt                                                   <  2:1583171/93‑1 (MQ=255)
                                        gAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTaa                                                 <  2:2381319/93‑1 (MQ=255)
                                         aaaaTTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAAc                                                >  1:4167882/1‑93 (MQ=255)
                                         aaaaTTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAAc                                                <  2:3912797/93‑1 (MQ=255)
                                         aaaaTTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAAc                                                <  1:4455989/93‑1 (MQ=255)
                                         aaaaTTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAAc                                                <  2:3768559/93‑1 (MQ=255)
                                          aaaTTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAAcc                                               >  2:2266070/1‑93 (MQ=255)
                                           aaTTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCGGTCGCCTGCGGCTGGTAACCg                                              <  1:4021865/93‑1 (MQ=255)
                                           aaTTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCg                                              >  2:2016936/1‑93 (MQ=255)
                                           aaTTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCg                                              >  1:3739146/1‑93 (MQ=255)
                                             ttGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAATCGCa                                            >  1:3596122/1‑93 (MQ=255)
                                               gAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCaaa                                          <  1:2866430/93‑1 (MQ=255)
                                                aGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAg                                         <  1:1504285/93‑1 (MQ=255)
                                                aGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAg                                         >  2:4251803/1‑93 (MQ=255)
                                                 ggggTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGc                                        >  2:2556657/1‑93 (MQ=255)
                                                   ggTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGcac                                      >  2:101713/1‑93 (MQ=255)
                                                    gTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGcaca                                     >  2:4323277/1‑93 (MQ=255)
                                                    gTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGcaca                                     >  1:3241333/1‑93 (MQ=255)
                                                     ttGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCacac                                    <  1:1252457/93‑1 (MQ=255)
                                                      tGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACt                                   >  2:300137/1‑93 (MQ=255)
                                                      tGATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACt                                   <  2:2772206/93‑1 (MQ=255)
                                                       gATCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTg                                  >  1:265191/1‑93 (MQ=255)
                                                        aTCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGt                                 >  2:4425158/1‑93 (MQ=255)
                                                         tCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGta                                <  2:3518303/93‑1 (MQ=255)
                                                         tCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGta                                <  1:2226073/93‑1 (MQ=255)
                                                         tCACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGta                                <  1:299382/93‑1 (MQ=255)
                                                          cACGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGtat                               >  1:309037/1‑93 (MQ=255)
                                                           aCGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGtatt                              >  1:3406273/1‑93 (MQ=255)
                                                            cGTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGtatta                             <  1:2975765/93‑1 (MQ=255)
                                                            cGTTTTGTACTGAATTGCAGATAACAAAAAGCCACGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGtatta                             >  1:862759/1‑93 (MQ=255)
                                                             gTTTTGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGtattat                            <  2:622246/93‑1 (MQ=255)
                                                              ttttGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATg                           >  2:825999/1‑93 (MQ=255)
                                                                ttGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTc                         >  1:2672209/1‑93 (MQ=255)
                                                                 tGTACTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCa                        >  2:749141/1‑93 (MQ=255)
                                                                    aCTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCAAACTGTATTATGTCAaca                     >  2:1223834/1‑93 (MQ=255)
                                                                     cTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAacac                    <  2:673873/93‑1 (MQ=255)
                                                                     cTGAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAacac                    >  1:204614/1‑93 (MQ=255)
                                                                      tGAATTGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACt                   >  2:4526716/1‑93 (MQ=255)
                                                                       gAATTGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACta                  >  1:2287259/1‑92 (MQ=255)
                                                                          ttGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGaaa               >  2:2431542/1‑93 (MQ=255)
                                                                          ttGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGaaa               <  1:3322249/93‑1 (MQ=255)
                                                                          ttGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGaaa               <  2:4737497/93‑1 (MQ=255)
                                                                           tGCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAg              <  2:1638557/93‑1 (MQ=255)
                                                                            gCAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGt             <  1:3722368/93‑1 (MQ=255)
                                                                             cAGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGta            <  1:1666104/93‑1 (MQ=255)
                                                                             cAGATAACAAAAAACCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGta            <  1:853214/93‑1 (MQ=255)
                                                                              aGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGtat           >  2:1877562/1‑93 (MQ=255)
                                                                              aGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGtat           >  1:4669027/1‑93 (MQ=255)
                                                                              aGATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGtat           <  2:1797844/93‑1 (MQ=255)
                                                                               gATAACAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGtata          <  1:1125999/93‑1 (MQ=255)
                                                                                  aaCAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACgt       <  1:217389/93‑1 (MQ=255)
                                                                                  aaCAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTAGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACAATGAAAGGATACgt       >  1:1498543/1‑93 (MQ=39)
                                                                                   aCAAAAAGCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACgtg      <  1:4484889/93‑1 (MQ=255)
                                                                                      aaaagCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACGTGTTc   <  1:2223310/93‑1 (MQ=255)
                                                                                       aaagCCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACGTGTTcc  >  1:1309643/1‑93 (MQ=255)
                                                                                          |                                                                                         
CAGCTTTACAGGGGAGACAATGGAAAATTTTTCAGCAAGGGAAAATTGAGGGGTTGATCACGTTTTGTACTGAATTGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTCGTCAGTCGCCTGCGGCTGGTAACCGCAAAGCACACTGTATTATGTCAACACTGAAAGTATACGTGTTCC  >  NC_000913/3850660‑3850839

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