Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,561,503 A→G intergenic (+366/‑94) b2447 → / → b2448 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,561,5030AG99.4% 525.7 / ‑4.7 170intergenic (+366/‑94)b2447/b2448hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (107/62);  total (108/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

Only 100 of 175 total aligned reads displayed.

CACCAAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTG  >  NC_000913/2561413‑2561592
                                                                                          |                                                                                         
caccaAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGt                                                                                         <  2:4333197/93‑1 (MQ=255)
 accaAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGtt                                                                                        <  2:963794/93‑1 (MQ=255)
 accaAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGtt                                                                                        >  1:2061238/1‑93 (MQ=255)
 accaAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGtt                                                                                        >  1:2944961/1‑93 (MQ=255)
  ccaAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTc                                                                                       >  2:270703/1‑93 (MQ=255)
   caAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTcc                                                                                      >  2:2272186/1‑93 (MQ=255)
    aaaGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTccc                                                                                     <  2:3643641/93‑1 (MQ=255)
     aaGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCca                                                                                    >  1:4049443/1‑93 (MQ=255)
     aaGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCca                                                                                    >  2:3878275/1‑93 (MQ=255)
      aGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTGCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaa                                                                                   >  1:2684698/1‑93 (MQ=255)
       gCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaac                                                                                  <  2:2188268/93‑1 (MQ=255)
         caacaaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaa                                                                                <  1:3986639/93‑1 (MQ=255)
         caacaaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaa                                                                                <  1:517572/93‑1 (MQ=255)
           acaaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaaca                                                                              >  2:3697698/1‑93 (MQ=255)
            catgATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaacaa                                                                             <  2:3393523/90‑1 (MQ=255)
             aaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACaaa                                                                            <  1:3912150/93‑1 (MQ=255)
               gATGGGAGGTGAATTTTGTGATGGAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAATGAAAGGCGGTTCCCAACAACAAAAt                                                                          >  2:3936657/1‑93 (MQ=255)
                 tGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATaa                                                                        >  1:1739229/1‑93 (MQ=255)
                  gggAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATaaa                                                                       <  1:872106/93‑1 (MQ=255)
                  gggAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATaaa                                                                       <  2:3830630/93‑1 (MQ=255)
                  gggAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATaaa                                                                       >  2:2527856/1‑93 (MQ=255)
                   ggAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAg                                                                      <  1:3083912/93‑1 (MQ=255)
                    gAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGa                                                                     <  2:1747192/93‑1 (MQ=255)
                     aGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGAt                                                                    <  1:301278/93‑1 (MQ=255)
                      ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATg                                                                   >  2:3550324/1‑93 (MQ=255)
                       gTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATgg                                                                  <  1:4092040/93‑1 (MQ=255)
                         gAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGa                                                                >  1:2659608/1‑93 (MQ=255)
                          aaTTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAg                                                               <  2:1097573/93‑1 (MQ=255)
                           aTTTTGTGGTGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGCgg                                                              >  2:1917974/1‑93 (MQ=255)
                           aTTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAgg                                                              >  1:2783763/1‑93 (MQ=255)
                             tttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTa                                                            >  1:968831/1‑93 (MQ=255)
                             tttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTa                                                            >  2:3281668/1‑93 (MQ=255)
                             tttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTa                                                            >  1:1539670/1‑93 (MQ=255)
                               tgtgATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTaaa                                                          >  1:4365607/1‑93 (MQ=255)
                                 tgATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaa                                                        <  2:2588548/93‑1 (MQ=255)
                                  gATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAt                                                       >  2:2910051/1‑93 (MQ=255)
                                   aTGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATg                                                      >  1:3949465/1‑93 (MQ=255)
                                    tGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGa                                                     >  2:312676/1‑93 (MQ=255)
                                     gCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGaa                                                    >  2:4110628/1‑93 (MQ=255)
                                      cAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGaaa                                                   <  1:3804452/93‑1 (MQ=255)
                                       aGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAAc                                                  >  2:2938500/1‑93 (MQ=255)
                                       aGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAAc                                                  <  1:3440766/93‑1 (MQ=255)
                                       aGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAAc                                                  >  2:4146050/1‑93 (MQ=255)
                                       aGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAAc                                                  >  2:582699/1‑93 (MQ=255)
                                       aGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAAc                                                  >  1:1245379/1‑93 (MQ=255)
                                        gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACg                                                 >  1:1194253/1‑93 (MQ=255)
                                        gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACg                                                 >  1:250047/1‑93 (MQ=255)
                                         aTCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACgg                                                >  1:464343/1‑93 (MQ=255)
                                          tCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGt                                               >  2:2796246/1‑93 (MQ=255)
                                            gCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTgg                                             >  2:2419607/1‑93 (MQ=255)
                                              tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGa                                           >  1:2918411/1‑93 (MQ=255)
                                              tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGa                                           >  2:3052740/1‑93 (MQ=255)
                                              tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGa                                           >  1:13128/1‑93 (MQ=255)
                                              tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGa                                           <  1:60833/93‑1 (MQ=255)
                                              tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGGAAAAATGAAACGGTGGGa                                           >  2:1551743/1‑93 (MQ=255)
                                                  tCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTg                                       >  2:1259651/1‑93 (MQ=255)
                                                  tCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTg                                       <  2:16045/93‑1 (MQ=255)
                                                   cccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTgg                                      >  2:2990365/1‑93 (MQ=255)
                                                     cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGTAAGTGGCt                                    <  2:4573921/93‑1 (MQ=255)
                                                      ttCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCtt                                   >  2:2779130/1‑93 (MQ=255)
                                                      ttCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCtt                                   >  1:532915/1‑93 (MQ=255)
                                                        cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTa                                 <  1:3284608/93‑1 (MQ=255)
                                                         aGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTAc                                >  2:4159455/1‑93 (MQ=255)
                                                         aGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTAc                                <  2:1647201/93‑1 (MQ=255)
                                                          ggCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACt                               <  1:2966340/93‑1 (MQ=255)
                                                           gCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTg                              >  1:2862133/1‑93 (MQ=255)
                                                           gCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTg                              >  2:1537028/1‑93 (MQ=255)
                                                             aTTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGt                            >  1:1645992/1‑93 (MQ=255)
                                                              tttttGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGATGTGGCTTTACTGGTg                           >  1:4523594/1‑93 (MQ=255)
                                                              tttttGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                           >  2:4564249/1‑93 (MQ=255)
                                                               ttttGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGa                          >  1:919412/1‑93 (MQ=255)
                                                                 ttGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATg                        <  2:3913505/93‑1 (MQ=255)
                                                                  tGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATgg                       <  2:2026048/93‑1 (MQ=255)
                                                                   gTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATggg                      >  1:1745050/1‑92 (MQ=255)
                                                                   gTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGt                      >  1:2064600/1‑93 (MQ=255)
                                                                    tttAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTa                     >  1:689039/1‑93 (MQ=255)
                                                                      tAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAAt                   <  1:3832071/93‑1 (MQ=255)
                                                                      tAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAAt                   >  2:1156677/1‑93 (MQ=255)
                                                                       aaTTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATg                  >  2:388128/1‑93 (MQ=255)
                                                                       aaTTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATg                  >  1:1190263/1‑93 (MQ=255)
                                                                        aTTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGt                 >  1:230259/1‑93 (MQ=255)
                                                                         ttCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                >  2:1703766/1‑93 (MQ=255)
                                                                         ttCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                >  1:3757799/1‑93 (MQ=255)
                                                                         ttCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                >  2:243533/1‑93 (MQ=255)
                                                                          tcaatcaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTg               >  1:3122520/1‑93 (MQ=255)
                                                                           caatcaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTgg              >  1:3727906/1‑93 (MQ=255)
                                                                             atcaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtt            <  1:435165/93‑1 (MQ=255)
                                                                                aaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGc         >  2:264398/1‑93 (MQ=255)
                                                                                 aCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCt        >  2:1011223/1‑93 (MQ=255)
                                                                                 aCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGCGGGAAGTGGCTTCACTGGCGATGGTAATGTTGGTTTGCt        >  2:2048454/1‑93 (MQ=255)
                                                                                   gAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGt      >  1:4470283/1‑93 (MQ=255)
                                                                                     aaGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt    <  1:4115409/93‑1 (MQ=255)
                                                                                     aaGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt    >  2:3996042/1‑93 (MQ=255)
                                                                                     aaGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt    <  2:3037007/93‑1 (MQ=255)
                                                                                     aaGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt    >  2:1366515/1‑93 (MQ=255)
                                                                                      aGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAtt   <  2:3244236/93‑1 (MQ=255)
                                                                                       ggCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTg  <  1:47665/93‑1 (MQ=255)
                                                                                       ggCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTg  >  2:4434501/1‑93 (MQ=255)
                                                                                       ggCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTg  >  2:1762625/1‑93 (MQ=255)
                                                                                       ggCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTg  >  2:3351942/1‑93 (MQ=255)
                                                                                          |                                                                                         
CACCAAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTG  >  NC_000913/2561413‑2561592

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