Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,673,861 A→G intergenic (+123/+59) yhjE → / ← yhjG putative transport protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,673,8610AG84.0% 815.9 / 138.7 331intergenic (+123/+59)yhjE/yhjGputative transport protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/21);  new base G (149/129);  total (181/150)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

Only 100 of 342 total aligned reads displayed.

TCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGAC  >  NC_000913/3673770‑3673951
                                                                                           |                                                                                          
 cAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCt                                                                                          >  2:3012578/1‑93 (MQ=255)
  aGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACtt                                                                                         <  1:3824452/93‑1 (MQ=255)
  aGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACtt                                                                                         >  2:3519319/1‑93 (MQ=255)
   ggCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCttt                                                                                        >  1:4613427/1‑93 (MQ=255)
   ggCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGCTCACGCCGCATCCGGCATGAACAAAGCACGCttt                                                                                        >  1:4484560/1‑93 (MQ=255)
    gCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTg                                                                                       >  1:5439854/1‑93 (MQ=255)
    gCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTg                                                                                       >  2:4502749/1‑93 (MQ=255)
     ccTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCttt                                                                                        <  1:6618906/91‑1 (MQ=255)
     ccTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGc                                                                                      >  1:5255635/1‑93 (MQ=255)
      cTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGTCATGAACAAAGCACACTTTGcc                                                                                     <  1:3270379/93‑1 (MQ=255)
      cTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGcc                                                                                     <  1:2462275/93‑1 (MQ=255)
        gAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCcaa                                                                                   >  2:1576672/1‑93 (MQ=255)
        gAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCcaa                                                                                   <  1:6192346/93‑1 (MQ=255)
        gAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCcaa                                                                                   >  2:5481536/1‑93 (MQ=255)
           tttCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCATGAACAAAGCACGCTTTGCcaacaa                                                                                >  1:5305389/1‑93 (MQ=255)
            ttCCATACAATATATTTAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAAt                                                                               <  1:2394881/93‑1 (MQ=255)
            ttCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAAt                                                                               <  1:2561718/93‑1 (MQ=255)
            ttCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAAt                                                                               <  2:5747189/93‑1 (MQ=255)
               cATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTa                                                                            >  1:764899/1‑93 (MQ=255)
                 tACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTaaa                                                                          >  2:5834956/1‑93 (MQ=255)
                   cAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATc                                                                        <  2:3907762/93‑1 (MQ=255)
                   cAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATc                                                                        <  1:974397/93‑1 (MQ=255)
                    aaTATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCt                                                                       <  2:1238991/93‑1 (MQ=255)
                      tatatTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGt                                                                     <  1:2364704/93‑1 (MQ=255)
                       atatTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATATGtt                                                                    <  2:2954160/93‑1 (MQ=255)
                        tatTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGttt                                                                   >  1:3661048/1‑93 (MQ=255)
                        tatTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGttt                                                                   >  2:599199/1‑93 (MQ=255)
                         atTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTa                                                                  >  2:3093901/1‑93 (MQ=255)
                         atTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTa                                                                  >  2:2483770/1‑93 (MQ=255)
                          ttGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTaa                                                                 >  2:5364712/1‑93 (MQ=255)
                           tGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAAt                                                                <  2:877526/93‑1 (MQ=255)
                             aaTTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTg                                                              >  2:4225786/1‑93 (MQ=255)
                              aTTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGc                                                             <  1:164464/93‑1 (MQ=255)
                              aTTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGc                                                             <  2:5541813/93‑1 (MQ=255)
                               ttCGAAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGcc                                                            <  2:4434838/93‑1 (MQ=255)
                                tCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCg                                                           <  2:3458856/93‑1 (MQ=255)
                                   cAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGAt                                                        <  1:4088198/93‑1 (MQ=255)
                                   cAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGAt                                                        <  1:3167736/93‑1 (MQ=255)
                                   cAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGAt                                                        <  1:6351699/93‑1 (MQ=255)
                                      ggATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTc                                                     >  2:2118545/1‑93 (MQ=255)
                                       gATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCa                                                    >  1:5642160/1‑93 (MQ=255)
                                         tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGa                                                  <  2:2226911/93‑1 (MQ=255)
                                         tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGa                                                  <  1:1068134/93‑1 (MQ=255)
                                          ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGAc                                                 >  2:3733581/1‑93 (MQ=255)
                                           tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACa                                                >  1:1602542/1‑93 (MQ=255)
                                            gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACAt                                               >  1:1028506/1‑93 (MQ=255)
                                            gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACAt                                               >  1:6464124/1‑93 (MQ=255)
                                            gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACAt                                               >  2:493961/1‑93 (MQ=255)
                                              aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATcc                                             <  2:6190877/93‑1 (MQ=255)
                                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTGAGACATCCg                                            >  1:2135215/1‑93 (MQ=255)
                                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCg                                            >  1:852903/1‑93 (MQ=255)
                                                 ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGc                                          <  1:5537986/93‑1 (MQ=255)
                                                  cGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCa                                         >  2:4597429/1‑93 (MQ=255)
                                                   ggATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCaa                                        >  1:5172650/1‑93 (MQ=255)
                                                   ggATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCaa                                        <  1:1145560/93‑1 (MQ=255)
                                                    gATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAAt                                       <  2:6249183/93‑1 (MQ=255)
                                                    gATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAAt                                       >  1:4193073/1‑93 (MQ=255)
                                                    gATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAAt                                       >  1:1290405/1‑93 (MQ=255)
                                                    gATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAAt                                       >  1:6558199/1‑93 (MQ=255)
                                                     aTAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAAtt                                      >  1:5030322/1‑93 (MQ=255)
                                                     aTAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAAtt                                      >  2:2771032/1‑93 (MQ=255)
                                                      tAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTa                                     <  1:4793907/93‑1 (MQ=255)
                                                         ggCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTatca                                  >  1:4654994/1‑93 (MQ=255)
                                                         ggCGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTatca                                  >  2:6605561/1‑93 (MQ=255)
                                                           cGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTatcatc                                >  1:2679491/1‑93 (MQ=255)
                                                           cGTTCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTatcatc                                <  1:2794707/93‑1 (MQ=255)
                                                            gTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTAtcatca                               >  1:4738419/1‑93 (MQ=255)
                                                            gTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTAtcatca                               <  2:3257267/93‑1 (MQ=255)
                                                             ttCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCAc                              <  2:6049681/93‑1 (MQ=255)
                                                             ttCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCCGCAATTATCATCAc                              >  2:4716920/1‑93 (MQ=255)
                                                              tCACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCAct                             <  1:3524694/93‑1 (MQ=255)
                                                               cACGCCGCATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATTActt                            >  1:3423069/1‑93 (MQ=255)
                                                                aCGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCActtt                           >  1:6066641/1‑92 (MQ=255)
                                                                 cgccgcATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCActtct                          <  1:5676107/93‑1 (MQ=255)
                                                                  gccgcATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCActtctt                         >  2:6411984/1‑93 (MQ=255)
                                                                  gccgcATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCActtctt                         >  2:1700964/1‑93 (MQ=255)
                                                                   ccgcATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACttcttc                        <  2:1470691/93‑1 (MQ=255)
                                                                   ccgcATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACttcttc                        <  1:1777690/93‑1 (MQ=255)
                                                                   ccgcATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACttcttc                        <  1:1956556/93‑1 (MQ=255)
                                                                    cgcATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCa                       >  1:831265/1‑93 (MQ=255)
                                                                     gcATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCAt                      >  2:6095661/1‑93 (MQ=255)
                                                                     gcATCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCAt                      >  1:4025594/1‑93 (MQ=255)
                                                                     gcATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCAt                      >  2:1893880/1‑93 (MQ=255)
                                                                       aTCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCt                    <  2:2346764/93‑1 (MQ=255)
                                                                        tCCGGCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTg                   >  1:5768226/1‑93 (MQ=255)
                                                                           ggCATGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGa                <  1:2305403/93‑1 (MQ=255)
                                                                           ggCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGa                >  2:6525959/1‑93 (MQ=255)
                                                                             cATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACa              >  1:436808/1‑93 (MQ=255)
                                                                              aTGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACaa             >  1:1408299/1‑93 (MQ=255)
                                                                               tGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACaaa            >  2:1047000/1‑93 (MQ=255)
                                                                               tGAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACaaa            >  2:3134115/1‑93 (MQ=255)
                                                                                gAACAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACaaaa           >  2:4891513/1‑93 (MQ=255)
                                                                                 aaCAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAAt          <  2:4041228/93‑1 (MQ=255)
                                                                                 aaCAAAGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAAt          <  1:5914122/93‑1 (MQ=255)
                                                                                 aaCAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAAt          >  1:4540468/1‑93 (MQ=255)
                                                                                 aaCAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAAt          <  1:2957637/93‑1 (MQ=255)
                                                                                  aCAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATc         <  1:2197603/93‑1 (MQ=255)
                                                                                     aaGCACGCTTTGCAAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTc      <  2:3216709/93‑1 (MQ=255)
                                                                                      aacacgCTTTGCCAACAATCTAAATCTGTTTAATTGCCCGATGACAGACATCCGACAATTATCATCACTTCTTCTTCTGCGACAAAATCGTcc     >  1:323439/3‑93 (MQ=255)
                                                                                      aGCACGCTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTcc     <  1:492778/93‑1 (MQ=255)
                                                                                           |                                                                                          
TCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACATCCGGCAATTATCATCACTTCTTCATCTGCGACAAAATCGTCCGAC  >  NC_000913/3673770‑3673951

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