Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,490,055 T→G 100% E3D (GAA→GAC yjgR ← DUF853 family protein with NTPase fold

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,490,0550TG96.2% 568.4 / 6.6 165E3D (GAA→GACyjgRDUF853 family protein with NTPase fold
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base G (76/82);  total (78/87)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

TACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTG  >  NC_000913/4489958‑4490150
                                                                                                 |                                                                                               
 aCCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCa                                                                                              >  2:400368/1‑100 (MQ=255)
  cccGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCAc                                                                                             <  1:1154849/100‑1 (MQ=255)
  cccGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCAc                                                                                             <  1:1016921/100‑1 (MQ=255)
   ccGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACt                                                                                            >  2:2252707/1‑100 (MQ=255)
       ccccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATg                                                                                        <  2:2439021/100‑1 (MQ=255)
       ccccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATg                                                                                        >  2:1870686/1‑100 (MQ=255)
        cccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATga                                                                                       <  1:2441491/100‑1 (MQ=255)
        cccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATga                                                                                       <  1:3750752/100‑1 (MQ=255)
        cccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATga                                                                                       <  1:1187241/100‑1 (MQ=255)
         ccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                      >  1:172549/1‑100 (MQ=255)
         ccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                      <  1:4174903/100‑1 (MQ=255)
         ccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                      <  2:1011081/100‑1 (MQ=255)
          cGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATgaga                                                                                     <  2:4074875/100‑1 (MQ=255)
            ttGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGata                                                                                   <  2:3654063/100‑1 (MQ=255)
            ttGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGata                                                                                   <  1:320724/100‑1 (MQ=255)
             tGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGatat                                                                                  >  2:4172105/1‑100 (MQ=255)
              gcgcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtc                                                                                 >  1:444147/1‑100 (MQ=255)
               cgcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcc                                                                                >  1:3097554/1‑100 (MQ=255)
               cgcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcc                                                                                <  2:2508743/100‑1 (MQ=255)
               cgcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcc                                                                                >  2:15543/1‑100 (MQ=255)
                gcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcct                                                                               >  2:1720341/1‑100 (MQ=255)
                  cGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcctcc                                                                             <  1:1468267/100‑1 (MQ=255)
                    gTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAg                                                                           <  1:1706603/100‑1 (MQ=255)
                     tAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGt                                                                          >  2:1447836/1‑100 (MQ=255)
                      aaTCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGtt                                                                         >  2:792222/1‑100 (MQ=255)
                        tCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTc                                                                       >  1:634999/1‑100 (MQ=255)
                        tCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCACGAGATATCCTCCAGTTTc                                                                       <  1:4005105/100‑1 (MQ=255)
                         cAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcc                                                                      <  2:2613836/100‑1 (MQ=255)
                         cAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcc                                                                      <  1:3967677/100‑1 (MQ=255)
                         cAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcc                                                                      <  2:3065153/100‑1 (MQ=255)
                         cAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTCTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcc                                                                      <  1:2734045/100‑1 (MQ=255)
                          aGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcct                                                                     >  2:2159588/1‑100 (MQ=255)
                          aGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcct                                                                     <  1:3286784/100‑1 (MQ=255)
                               aTGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcctg                                                                <  2:3271590/100‑1 (MQ=255)
                                 ggtcgggaGCCATGGCCGGCAGGAAAAGCAGGGGTGGTGCGGGCGTGCGGGGAATGAAGAGGGGGTCACGCATGAGATATCCTCCAGTTTcctgcc                                                                  <  1:4075239/89‑1 (MQ=255)
                                  aCGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAg                                                             <  2:927530/100‑1 (MQ=255)
                                  aCGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCCTGCGGGCAATTAACAGGGGGTCACTCATGAGATCTACTCCAGTTTCCTGCCTGGAg                                                             <  1:2489829/100‑1 (MQ=255)
                                   cGGTTAGCCATTCCCGGTAGTAAAAACAGTGCTGTGTCCGGCGTGCGGGCCATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGt                                                            <  1:2747478/100‑1 (MQ=255)
                                   cGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATTCTCCAGTTTCCTGCCTGGAGt                                                            <  1:3049680/100‑1 (MQ=255)
                                      ttAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGCCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGtata                                                         >  1:2602388/1‑100 (MQ=255)
                                         gCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                      >  2:1831985/1‑100 (MQ=255)
                                         gCCATTCCCGGCAGTAAAAACAGTTCTGGGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                      <  1:250112/100‑1 (MQ=255)
                                           cATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTCACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCa                                                     <  1:2143893/99‑1 (MQ=255)
                                            aTTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAAc                                                   >  2:1034930/1‑100 (MQ=255)
                                             ttCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACg                                                  <  2:749385/100‑1 (MQ=255)
                                               cccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTc                                                >  1:116858/1‑100 (MQ=255)
                                               cccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGt                                                 >  1:3684963/1‑99 (MQ=255)
                                               cccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTc                                                >  1:4103295/1‑100 (MQ=255)
                                               cccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTc                                                >  1:158240/1‑100 (MQ=255)
                                               cccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTc                                                >  1:1386377/1‑100 (MQ=255)
                                                  ggCAGTAAAAACAGTTCTGTGTCCGGCTTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCgg                                             <  1:205272/100‑1 (MQ=255)
                                                   gCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGc                                            <  2:146379/100‑1 (MQ=255)
                                                     aGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCaa                                          <  1:3276059/100‑1 (MQ=255)
                                                      gTAAAAACAGTTCAGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCaaa                                         >  2:2904033/1‑100 (MQ=255)
                                                       tAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAag                                        <  2:3708064/100‑1 (MQ=255)
                                                        aaaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAaga                                       >  2:463468/1‑100 (MQ=255)
                                                        aaaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAaga                                       >  2:1630302/1‑100 (MQ=255)
                                                        aaaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAaga                                       <  1:581735/100‑1 (MQ=255)
                                                        aaaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGCGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAaga                                       <  1:437791/100‑1 (MQ=255)
                                                         aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCATTTTCTTACTTGGGGATAAGCCACGTTCGGCgaa                                         >  1:2304221/1‑96 (MQ=255)
                                                         aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                      <  2:3008581/100‑1 (MQ=255)
                                                          tccgatctCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAgaga                                     <  2:1994266/93‑1 (MQ=255)
                                                          aaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAGCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAgaga                                     >  2:4054371/1‑100 (MQ=255)
                                                          aaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAgaga                                     >  2:3114405/1‑100 (MQ=255)
                                                           aaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGaa                                    >  2:2989608/1‑100 (MQ=255)
                                                           aaCAGTTCTGTGTCCGGCGTGCGGGCAATGAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGaa                                    <  2:4159963/100‑1 (MQ=255)
                                                             cAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATa                                  <  2:116858/100‑1 (MQ=255)
                                                               gTTCTGTGTCTGGCGTGCGGGCAATTACCAGGGGGTTACGCATGAGGTATAAGCCGGTTTCCTGGGCGGAGAATAGGCAACGTGCGGCGAAGCGAAt                                   >  1:2366864/1‑97 (MQ=255)
                                                               gTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACg                                >  2:1609108/1‑100 (MQ=255)
                                                               gTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACGCAGGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGGAAAGTCCGGCAAAGAGAATagg                                >  1:1085084/1‑98 (MQ=255)
                                                                ttCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACgt                               <  2:1323065/100‑1 (MQ=255)
                                                                 tCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTTACTCATGAGATATCCTACAGTTTCCTGCATGGCGTATAGGCAACGTCCGGCAAAGAGAATACgtg                              <  1:3840894/100‑1 (MQ=255)
                                                                 tCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACgtg                              <  2:2380507/100‑1 (MQ=255)
                                                                 tCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACgtg                              <  2:557230/100‑1 (MQ=255)
                                                                 tCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGAAAAGAGAATACgt                               >  1:1017750/1‑99 (MQ=255)
                                                                  cTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGaagatcgga                             >  1:1994262/1‑93 (MQ=255)
                                                                  cTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACgtgt                             <  2:1136849/100‑1 (MQ=255)
                                                                      gtCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATg                         >  2:817856/1‑100 (MQ=255)
                                                                          ggCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAaga                     <  1:1074781/100‑1 (MQ=255)
                                                                          ggCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAaga                     >  1:1749064/1‑100 (MQ=255)
                                                                            cGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAgaga                   <  1:4054351/100‑1 (MQ=255)
                                                                             gTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaa                  <  2:2742218/100‑1 (MQ=255)
                                                                              tGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaaa                 <  2:623702/100‑1 (MQ=255)
                                                                              tGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaaa                 >  1:1851263/1‑100 (MQ=255)
                                                                                  ggCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATaa             >  1:1831160/1‑100 (MQ=255)
                                                                                  ggCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATaa             >  1:3579631/1‑100 (MQ=255)
                                                                                  ggCAATTAACAGGGGGCCCCTCATGAGATATCCTCCAGTTACCTGCCGGGAGTCTAGGCAACGTCCGGCAAAGCGAATACGAGTAATGAAGAGAAAATa              >  1:3871872/1‑99 (MQ=255)
                                                                                    cAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           <  2:3811064/100‑1 (MQ=255)
                                                                                    cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           >  2:3299770/1‑100 (MQ=255)
                                                                                    cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           <  2:2958673/100‑1 (MQ=255)
                                                                                    cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           >  1:2382763/1‑100 (MQ=255)
                                                                                    cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           <  2:3767565/100‑1 (MQ=255)
                                                                                      aTTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTa         >  1:1359101/1‑100 (MQ=255)
                                                                                      aTTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTa         >  1:348977/1‑100 (MQ=255)
                                                                                       ttAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTaa        >  2:751148/1‑100 (MQ=255)
                                                                                       ttAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTaa        >  1:296070/1‑100 (MQ=255)
                                                                                       ttAACAGGGGGTCACTCATGAGATATCCTCCAGTCTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGCATACGTGTAATGAAGAGAAAATACCCTaa        >  2:355334/1‑100 (MQ=255)
                                                                                          aCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGGCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGc     >  1:2285481/1‑100 (MQ=255)
                                                                                           cAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCg    >  1:1892244/1‑100 (MQ=255)
                                                                                             ggggGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCgtg  <  2:3214312/100‑1 (MQ=255)
                                                                                                 |                                                                                               
TACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTG  >  NC_000913/4489958‑4490150

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