Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,071,901 G→A A197T (GCC→ACC)  yihV → putative kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,071,9010GA99.6% 901.1 / ‑5.6 249A197T (GCC→ACC) yihVputative kinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (113/135);  total (113/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 253 total aligned reads displayed.

CCGCAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGGCCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAATGGT  >  NC_000913/4071811‑4071990
                                                                                          |                                                                                         
ccgcAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCACGCTTAACGGGCGTGAAAGAGATGAcc                                                                                         >  2:6163744/1‑93 (MQ=255)
 cgcAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCa                                                                                        >  2:1025300/1‑93 (MQ=255)
   cAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGt                                                                                      <  2:3467110/93‑1 (MQ=255)
    aGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTg                                                                                     >  2:4356301/1‑93 (MQ=255)
    aGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTg                                                                                     >  1:2555278/1‑93 (MQ=255)
      gATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTgcg                                                                                   >  1:251579/1‑93 (MQ=255)
       atatCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTgcgc                                                                                  >  2:27700/1‑93 (MQ=255)
        tatCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCt                                                                                 >  2:5110958/1‑93 (MQ=255)
        tatCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCt                                                                                 >  2:3899075/1‑93 (MQ=255)
        tatCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCt                                                                                 >  2:1526909/1‑93 (MQ=255)
         atCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTa                                                                                <  2:5552944/93‑1 (MQ=255)
         atCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTa                                                                                >  1:1650243/1‑93 (MQ=255)
            aGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTTACGGGCGTGAAAGAGATGACCAGTGCGCTaaaa                                                                             >  1:1101463/1‑93 (MQ=255)
            aGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTaaaa                                                                             >  2:6569852/1‑93 (MQ=255)
              tGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACa                                                                           <  2:2462110/93‑1 (MQ=255)
              tGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACa                                                                           <  2:4256995/93‑1 (MQ=255)
               gAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAg                                                                          <  1:5247664/93‑1 (MQ=255)
               gAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAg                                                                          <  2:4770146/93‑1 (MQ=255)
                aGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAgg                                                                         <  2:6146808/93‑1 (MQ=255)
                aGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAgg                                                                         <  2:3844561/93‑1 (MQ=255)
                  cTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGca                                                                       >  2:1236574/1‑93 (MQ=255)
                  cTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGca                                                                       <  1:81787/93‑1 (MQ=255)
                    ggtggCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGcaca                                                                     >  2:6498447/1‑93 (MQ=255)
                      tggCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACaaa                                                                   >  2:3757314/1‑93 (MQ=255)
                      tggCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACaaa                                                                   <  2:6012083/93‑1 (MQ=255)
                      tggCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACaaa                                                                   <  1:3575640/93‑1 (MQ=255)
                          aTTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCt                                                               >  1:5256375/1‑93 (MQ=255)
                             aaGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTcac                                                            >  1:1936895/1‑93 (MQ=255)
                              aGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTcaca                                                           >  2:490896/1‑93 (MQ=255)
                                  aTCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATg                                                       >  2:3407036/1‑93 (MQ=255)
                                   tCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATgg                                                      >  2:2082006/1‑93 (MQ=255)
                                    cACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGt                                                     >  1:4205182/1‑93 (MQ=255)
                                    cACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGt                                                     >  1:6095425/1‑93 (MQ=255)
                                    cACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGt                                                     >  2:2604172/1‑93 (MQ=255)
                                         ggCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGt                                                >  2:3304743/1‑93 (MQ=255)
                                         ggCCGTTTCAGAACCGGCTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGt                                                <  1:455561/93‑1 (MQ=255)
                                          gCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTc                                               <  2:3549980/93‑1 (MQ=255)
                                          gCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTc                                               >  1:4726883/1‑93 (MQ=255)
                                          gCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTc                                               <  2:2537488/93‑1 (MQ=255)
                                          gCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTc                                               >  2:1294524/1‑93 (MQ=255)
                                          gCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTc                                               >  1:657884/1‑93 (MQ=255)
                                            cTTTTCAGCACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTa                                             >  1:6521127/1‑93 (MQ=255)
                                            cTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCGTGTCTa                                             >  2:3590608/1‑93 (MQ=255)
                                            cTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTa                                             >  1:1244632/1‑93 (MQ=255)
                                            cTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTa                                             >  2:183821/1‑93 (MQ=255)
                                            cTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTa                                             <  1:4352625/93‑1 (MQ=255)
                                             ttttCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTAt                                            >  1:4426253/1‑93 (MQ=255)
                                             ttttCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTAt                                            <  2:5941175/93‑1 (MQ=255)
                                             ttttCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTAt                                            >  1:1711357/1‑93 (MQ=255)
                                                tCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTAtgtg                                         <  1:6321513/93‑1 (MQ=255)
                                                 cAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGa                                        >  2:723107/1‑93 (MQ=255)
                                                   gAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGAcc                                      >  1:4467326/1‑93 (MQ=255)
                                                   gAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAAAGGTCATGTCTATGTGAcc                                      >  1:1471265/1‑93 (MQ=255)
                                                   gAACCGGGGCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGAcc                                      >  1:6531310/1‑93 (MQ=255)
                                                    aaCCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGAccc                                     <  2:81675/93‑1 (MQ=255)
                                                    aaCCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGAccc                                     <  1:3647211/93‑1 (MQ=255)
                                                     aCCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCa                                    >  2:4416493/1‑93 (MQ=255)
                                                      ccGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAg                                   <  2:3014113/93‑1 (MQ=255)
                                                       cGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAgg                                  <  2:5680742/93‑1 (MQ=255)
                                                       cGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAgg                                  >  1:468189/1‑93 (MQ=255)
                                                        gggTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAggg                                 <  2:3713065/93‑1 (MQ=255)
                                                        gggTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAggg                                 <  1:5393721/93‑1 (MQ=255)
                                                        gggTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAggg                                 <  1:5142571/93‑1 (MQ=255)
                                                        gggTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAggg                                 >  1:4593322/1‑93 (MQ=255)
                                                          gTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTa                               >  1:209071/1‑93 (MQ=255)
                                                          gTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTa                               <  2:190804/93‑1 (MQ=255)
                                                          gTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTa                               <  1:6261316/93‑1 (MQ=255)
                                                             tGGCGCGCTTAACGGGCGTGAAAGAGATGGCCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAgcg                            <  2:2317380/93‑1 (MQ=255)
                                                               gcacgcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCa                          <  1:5446905/90‑1 (MQ=255)
                                                                cgcgcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAg                         <  1:38082/93‑1 (MQ=255)
                                                                 gcgcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAgg                        >  2:5632403/1‑93 (MQ=255)
                                                                  cgcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGc                       <  2:2217143/93‑1 (MQ=255)
                                                                  cgcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGc                       >  1:687573/1‑93 (MQ=255)
                                                                  cgcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGc                       <  1:5737774/93‑1 (MQ=255)
                                                                   gcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCt                      >  1:4382734/1‑93 (MQ=255)
                                                                   gcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCt                      <  1:5662598/93‑1 (MQ=255)
                                                                     ttAACGGGCGTGCAAGAGATGACCAGTGCGCTAAAACAGGCACAAAGGCTCACAAATGTTCATGTCTATGTGACCCAGTGTAGCGCAGGCTgt                    >  1:2705757/1‑92 (MQ=255)
                                                                      tAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCg                   <  2:410598/93‑1 (MQ=255)
                                                                       aaCGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGa                  >  1:943782/1‑93 (MQ=255)
                                                                        aCGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGAc                 >  1:5047173/1‑93 (MQ=255)
                                                                        aCGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGAc                 >  1:4278613/1‑93 (MQ=255)
                                                                          gggCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActg               >  1:2804849/1‑93 (MQ=255)
                                                                          gggCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActg               <  2:4484363/93‑1 (MQ=255)
                                                                           ggCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActgg              >  2:3164773/1‑93 (MQ=255)
                                                                            gCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActggc             <  2:1625469/93‑1 (MQ=255)
                                                                             cGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActggct            <  1:5983244/93‑1 (MQ=255)
                                                                             cGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActggct            <  2:3288835/93‑1 (MQ=255)
                                                                              gTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActggctg           >  2:3732042/1‑93 (MQ=255)
                                                                              gTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActggctg           <  2:5300214/93‑1 (MQ=255)
                                                                                 aaaGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGaa        <  1:4891728/93‑1 (MQ=255)
                                                                                   agagATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGaaaa      <  2:2405008/93‑1 (MQ=255)
                                                                                    gagaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGCCTGGCTAGAAAAt     <  1:2159211/93‑1 (MQ=255)
                                                                                    gagaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTTGCTGGAAAAt     <  1:3095674/93‑1 (MQ=255)
                                                                                    gagaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAt     <  1:2411785/93‑1 (MQ=255)
                                                                                    gagaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAt     >  1:1073798/1‑93 (MQ=255)
                                                                                    gagaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAt     <  1:3381688/93‑1 (MQ=255)
                                                                                    gagaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAt     >  1:2935698/1‑93 (MQ=255)
                                                                                     agaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAtg    <  1:3998955/93‑1 (MQ=255)
                                                                                       aTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAtggt  <  2:2827893/93‑1 (MQ=255)
                                                                                       aTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAtggt  <  2:791498/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCGCAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGGCCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAATGGT  >  NC_000913/4071811‑4071990

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