Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 985,663 A→G G181G (GGT→GGC ompF ← outer membrane protein 1a

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913985,6630AG98.9% 573.9 / ‑5.5 185G181G (GGT→GGCompFouter membrane protein 1a
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (109/74);  total (110/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

CCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAA  >  NC_000913/985574‑985751
                                                                                         |                                                                                        
ccAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGccc                                                                                       >  2:4550591/1‑93 (MQ=255)
 cAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCa                                                                                      >  2:2615955/1‑93 (MQ=255)
 cAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCa                                                                                      >  1:3629593/1‑93 (MQ=255)
   aCGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAgg                                                                                    >  2:269215/1‑93 (MQ=255)
   aCGATACCAAAGCCTTCGTATTCGTAACTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAgg                                                                                    >  1:4177243/1‑93 (MQ=255)
    cGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGt                                                                                   <  1:1490590/93‑1 (MQ=255)
    cGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGt                                                                                   >  2:501051/1‑93 (MQ=255)
     gATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTa                                                                                  >  1:1176938/1‑93 (MQ=255)
     gATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTa                                                                                  <  1:1894311/93‑1 (MQ=255)
      aTACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGCCACGCTCGTTTTTGCCCAgggcc                                                                                 >  2:3935943/1‑90 (MQ=255)
       tACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACt                                                                                >  2:3575854/1‑93 (MQ=255)
          cAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGaa                                                                             >  1:2324293/1‑93 (MQ=255)
            aaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACa                                                                           >  1:2318005/1‑93 (MQ=255)
            aaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACa                                                                           <  1:3505891/93‑1 (MQ=255)
            aaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGGGTCACGCTCGTTTTTGCCCAGGTACTGAACa                                                                           >  1:796191/1‑93 (MQ=255)
                cTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGa                                                                       >  2:2034401/1‑93 (MQ=255)
                cTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGa                                                                       >  1:2362709/1‑93 (MQ=255)
                cTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGCACAGCGa                                                                       >  1:4532372/1‑93 (MQ=255)
                 ttcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGaa                                                                      >  1:3043873/1‑93 (MQ=255)
                    gtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCGGGTACTGAACAGCGAAGtt                                                                   >  1:2539362/1‑93 (MQ=255)
                     tattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTc                                                                  >  1:4283598/1‑93 (MQ=255)
                     tattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTc                                                                  >  2:569933/1‑93 (MQ=255)
                     tattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTc                                                                  >  1:3162068/1‑93 (MQ=255)
                     tattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTc                                                                  >  1:3948289/1‑93 (MQ=255)
                     tattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACCGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAAAAGCGACGTTc                                                                  >  2:17508/1‑93 (MQ=255)
                       ttcgtaGCTGATAGACCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAg                                                                <  1:2268889/93‑1 (MQ=255)
                       ttcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAg                                                                >  2:2797679/1‑93 (MQ=255)
                         cgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGc                                                              <  1:1265675/93‑1 (MQ=255)
                           taGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCa                                                            >  1:252772/1‑93 (MQ=255)
                           taGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCa                                                            >  2:4040414/1‑93 (MQ=255)
                           taGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCa                                                            >  2:4545613/1‑93 (MQ=255)
                            aGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCAt                                                           >  1:1966121/1‑93 (MQ=255)
                             gCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATc                                                          >  2:696575/1‑93 (MQ=255)
                             gCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATc                                                          >  1:2094736/1‑93 (MQ=255)
                              cTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCa                                                         >  1:1567866/1‑93 (MQ=255)
                               tGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCaa                                                        >  2:375362/1‑93 (MQ=255)
                                  tAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCa                                                     >  1:836279/1‑93 (MQ=255)
                                   aGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCAAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAg                                                    >  1:435836/1‑93 (MQ=255)
                                    gAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGa                                                   >  2:3508395/1‑93 (MQ=255)
                                    gAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGa                                                   >  1:3324610/1‑93 (MQ=255)
                                    gAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGa                                                   >  1:786168/1‑93 (MQ=255)
                                     aaCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAc                                                  >  2:1214016/1‑93 (MQ=255)
                                     aaCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAc                                                  >  2:4083277/1‑93 (MQ=255)
                                      aCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAcc                                                 <  2:4495554/93‑1 (MQ=255)
                                         gCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCaaa                                              >  2:1489606/1‑93 (MQ=255)
                                         gCCAACACCGTCGCCGTTAGAACGGCGTGCAGCGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCaaa                                              >  1:3010306/1‑93 (MQ=255)
                                          ccAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaag                                             <  2:2199178/93‑1 (MQ=255)
                                           cAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaaga                                            <  1:1285151/93‑1 (MQ=255)
                                            aaCACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAaagaa                                           >  2:3410364/1‑93 (MQ=255)
                                            aaCACCGTCGCAGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGAGTAGTTCCGGCCATCACCCAGCCCaaacag                                           >  2:3074164/1‑90 (MQ=255)
                                              ccccGGCGCCGTTAGAACGGCGTGCAGGGGCCCGCGTGTTTTTGCCCAGGTTCTGAACAGCGAAGGTCAGGCCATCAACCAGACCAAAGAAGt                                         <  1:1017807/91‑1 (MQ=21)
                                              cacCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGt                                         <  1:1017808/93‑1 (MQ=255)
                                              cacCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGt                                         <  1:1100922/93‑1 (MQ=255)
                                              cacCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGt                                         <  1:1545710/93‑1 (MQ=255)
                                               acCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGtt                                        <  1:4041805/93‑1 (MQ=255)
                                               acCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGtt                                        <  2:1680664/93‑1 (MQ=255)
                                                ccGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTg                                       <  1:4177548/93‑1 (MQ=255)
                                                ccGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTg                                       >  2:136189/1‑93 (MQ=255)
                                                ccGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTg                                       <  1:3580614/93‑1 (MQ=255)
                                                ccGTCGCCGTTAGAACGGCGTGCAGTGGCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTg                                       >  2:1588722/1‑93 (MQ=255)
                                                ccGTCGCCGTTAGAACGGCGCGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTg                                       >  2:1588721/1‑93 (MQ=255)
                                                  gTCGCCGTTAGAACGGCGTGCAGTGTCACGCTAGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGa                                     <  1:3697299/93‑1 (MQ=255)
                                                  ctCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGa                                     <  1:2649171/92‑1 (MQ=255)
                                                    cGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGt                                   <  2:3261022/93‑1 (MQ=255)
                                                     gCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGtt                                  <  2:4179583/93‑1 (MQ=255)
                                                      ccGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTa                                 >  2:3204585/1‑93 (MQ=255)
                                                       cGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTAc                                <  2:3476528/93‑1 (MQ=255)
                                                       cGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTAc                                <  2:1003379/93‑1 (MQ=255)
                                                        gTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACg                               >  1:631602/1‑93 (MQ=255)
                                                         ttAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGa                              >  1:313967/1‑93 (MQ=255)
                                                         ttAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGa                              >  1:3199277/1‑93 (MQ=255)
                                                          tAGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGAt                             >  2:1994608/1‑93 (MQ=255)
                                                           aGAACGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATa                            <  1:2437983/93‑1 (MQ=255)
                                                              aCGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGt                         <  2:3556337/93‑1 (MQ=255)
                                                               cGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTATTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTa                        <  2:1900025/93‑1 (MQ=255)
                                                               cGGCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTa                        >  2:3202182/1‑93 (MQ=255)
                                                                ggCGTGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAg                       >  1:490440/1‑93 (MQ=255)
                                                                    tGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAc                   <  1:3653046/93‑1 (MQ=255)
                                                                    tGCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAc                   <  2:3911654/93‑1 (MQ=255)
                                                                    ggCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAc                   <  1:1653811/92‑1 (MQ=255)
                                                                     gCAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAcg                  >  2:3979096/1‑93 (MQ=255)
                                                                      cAGTGTCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAcgc                 <  1:4429163/93‑1 (MQ=255)
                                                                        gtgtCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAcgccg               >  1:28116/1‑93 (MQ=255)
                                                                         tgtCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAcgccgc              >  1:2035270/1‑93 (MQ=255)
                                                                         tgtCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAcgccgc              <  2:3673464/93‑1 (MQ=255)
                                                                          gtCCCGCTCGTTTTTGCCCAGGTACTGAAGAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACgccgcc             <  2:196952/93‑1 (MQ=255)
                                                                          gtCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACgccgcc             >  2:959839/1‑93 (MQ=255)
                                                                          gtCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACgccgcc             >  2:3832278/1‑93 (MQ=255)
                                                                          gtCACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACgccgcc             <  1:890295/93‑1 (MQ=255)
                                                                            cACGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCaa           >  1:4237494/1‑93 (MQ=255)
                                                                              cGCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAaca         <  2:1239918/93‑1 (MQ=255)
                                                                               gCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAacac        <  1:2188292/93‑1 (MQ=255)
                                                                               gCTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAacac        >  1:4345193/1‑93 (MQ=255)
                                                                                cTCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACg       >  2:424702/1‑93 (MQ=255)
                                                                                 tCGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGa      >  2:2726543/1‑93 (MQ=255)
                                                                                  cGTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGAc     <  2:4049034/93‑1 (MQ=255)
                                                                                   gTTTTTGCCCAGGTACTGAACATCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGAcc    >  1:850128/1‑93 (MQ=255)
                                                                                   gTTTTTGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGAcc    <  1:346558/93‑1 (MQ=255)
                                                                                    tttttGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGCAGCAACGCCGCCAACACGACCa   >  1:2845445/1‑93 (MQ=255)
                                                                                     ttttGCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCaa  <  1:3856935/93‑1 (MQ=255)
                                                                                         |                                                                                        
CCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAA  >  NC_000913/985574‑985751

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