Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,490,055 T→G 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,0550TG97.7% 441.2 / ‑4.4 129E3D (GAA→GAC) yjgRDUF853 family protein with NTPase fold
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (66/60);  total (68/61)
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 130 total aligned reads displayed.

CCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTG  >  NC_000913/4489960‑4490150
                                                                                               |                                                                                               
cccGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCAc                                                                                             >  1:3332143/1‑100 (MQ=255)
cccGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCAc                                                                                             <  2:1525401/100‑1 (MQ=255)
cccGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCAc                                                                                             >  1:1486743/1‑100 (MQ=255)
   gTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCa                                                                                          <  2:2578231/100‑1 (MQ=255)
     ccccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATg                                                                                        >  1:1431446/1‑100 (MQ=255)
      cccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGAGTTCGGGCAATTAACAGGGGGTCACACATga                                                                                       >  2:880954/1‑100 (MQ=255)
       ccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                      <  1:2478025/100‑1 (MQ=255)
       ccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                      <  1:689438/100‑1 (MQ=255)
       ccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                      >  2:449870/1‑100 (MQ=255)
       ccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                      <  1:39551/100‑1 (MQ=255)
        cGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgaga                                                                                     >  2:1337524/1‑100 (MQ=255)
         gTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGat                                                                                    >  2:1070166/1‑100 (MQ=255)
             cgcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcc                                                                                <  2:1564278/100‑1 (MQ=255)
              gcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAGCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcct                                                                               >  1:314581/1‑100 (MQ=255)
                cGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcctcc                                                                             <  2:351141/100‑1 (MQ=255)
                  gTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAg                                                                           <  2:30039/100‑1 (MQ=255)
                   tAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGt                                                                          >  2:3122496/1‑100 (MQ=255)
                       cAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcc                                                                      <  2:2002023/100‑1 (MQ=255)
                        aGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcct                                                                     >  2:2886207/1‑100 (MQ=255)
                        aGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcct                                                                     >  1:860717/1‑100 (MQ=255)
                         gCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctn                                                                    <  2:1710335/100‑2 (MQ=255)
                         gCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctg                                                                    <  1:450046/100‑1 (MQ=255)
                           ccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCGATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcc                                                                  >  1:2230420/1‑100 (MQ=255)
                           ccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcc                                                                  <  2:1122550/100‑1 (MQ=255)
                             aTGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcctg                                                                >  2:3212279/1‑100 (MQ=255)
                              tGACGGTTAGCCATTCCCGGCAGTAAAGACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTgg                                                               <  1:3174976/100‑1 (MQ=255)
                              tGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTgg                                                               >  1:1124439/1‑100 (MQ=255)
                                aCGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAg                                                             >  2:3254376/1‑100 (MQ=255)
                                aCGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAg                                                             <  1:3442176/100‑1 (MQ=255)
                                aCGGTTAGCCATTCCCGGCAGTAAAAACAGCTCCGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAg                                                             >  1:2283535/1‑100 (MQ=255)
                                 cGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGt                                                            >  2:1856165/1‑100 (MQ=255)
                                  ggTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGta                                                           <  2:2674509/100‑1 (MQ=255)
                                    ttAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGtata                                                         <  2:724405/100‑1 (MQ=255)
                                    ttAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGtata                                                         <  2:2704960/100‑1 (MQ=255)
                                     tAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAg                                                        >  2:319347/1‑100 (MQ=255)
                                      aGCCATTTCCGGCAGTAAAAACAGTGCTGTGCCCGGCGTGAGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAgg                                                       <  1:1888559/100‑1 (MQ=255)
                                      aGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGTGGTCAGTCATGAGATATCGTCCAGTTTCCTGCCTGGAGTATAgg                                                       <  1:2778581/100‑1 (MQ=255)
                                      aGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAgg                                                       <  2:3181020/100‑1 (MQ=255)
                                       gCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                      <  2:1863757/100‑1 (MQ=255)
                                        ccATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCGCTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCa                                                     >  2:2944392/1‑100 (MQ=255)
                                        ccATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCa                                                     <  2:2778703/100‑1 (MQ=255)
                                         cATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCaa                                                    <  2:1686060/100‑1 (MQ=255)
                                         cATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGATTCCTGCCTGGAGTATAGGCaa                                                    >  2:2903614/1‑100 (MQ=255)
                                          aTTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAAc                                                   <  1:621924/100‑1 (MQ=255)
                                          aTTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAAc                                                   >  2:1483993/1‑100 (MQ=255)
                                           ttCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACg                                                  <  1:156884/100‑1 (MQ=255)
                                           ttCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTACCTGGAGTATAGGCAACg                                                  <  1:1695111/100‑1 (MQ=255)
                                             cccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTc                                                >  2:2572388/1‑100 (MQ=255)
                                              ccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTcc                                               >  1:497597/1‑100 (MQ=255)
                                                  cAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCa                                           >  2:3114650/1‑100 (MQ=255)
                                                  cAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCa                                           <  1:2376607/100‑1 (MQ=255)
                                                   aGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTATCCTGCCTGCAGTATAGGCAACGTCCGGCa                                           >  1:3012482/1‑99 (MQ=255)
                                                   aGTAAAAACAGTTCTGTGTACGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCaa                                          >  2:451969/1‑100 (MQ=255)
                                                    gTAAAAACAGTTCAGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGAAATCCTCCAGTTTCCTACCTGGAGAATAGGCAACGACCGGCaaa                                         >  1:1651600/1‑100 (MQ=255)
                                                     tAAGAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAag                                        <  1:589885/100‑1 (MQ=255)
                                                     tAAAAACAGTTCTGTGTCCTGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAag                                        <  2:120059/100‑1 (MQ=255)
                                                     tAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAag                                        >  1:3010871/1‑100 (MQ=255)
                                                      aaaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAaga                                       >  1:1216551/1‑100 (MQ=255)
                                                       aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                      >  2:34525/1‑100 (MQ=255)
                                                       aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                      <  2:897531/100‑1 (MQ=255)
                                                       aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                      <  2:1547613/100‑1 (MQ=255)
                                                       aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                      >  1:3443174/1‑100 (MQ=255)
                                                        aaaaAGTTCTGTGGTCGGCGTGCGGGCAATTGACAGGGGGGCACTGATGAGATATACTCCACTTTGCTGCCGGGAGTATGGGCAACGTCCGGGACagag                                      <  1:1579098/99‑1 (MQ=255)
                                                        aaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAgaga                                     <  2:1638675/100‑1 (MQ=255)
                                                           cAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATa                                  <  1:1858274/100‑1 (MQ=255)
                                                           cAGTTCTGTGCCCGGCGTGCGGGCAACTAACAGGGGGTCACTCATGAGATATCCTCCAGTATCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATa                                  <  1:1319651/100‑1 (MQ=255)
                                                            aGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATAc                                 <  2:3127089/100‑1 (MQ=255)
                                                             gTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCGTGGAGTATAGGCAACGTCCGGCAAAGAGCAt                                   >  1:3139540/1‑97 (MQ=255)
                                                               tCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACgtg                              <  2:2061448/100‑1 (MQ=255)
                                                               tCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACgtg                              <  1:34525/100‑1 (MQ=255)
                                                               tCTGTGTCCGACGTGCGAACAATTAACAGGGGGTCACTAAAGACATATCCTACAGTTTACTGCCTGGAATATAGGCAACGTCCGGCCAAGAGAATAc                                 >  1:461070/1‑97 (MQ=255)
                                                                cTGTGTCCGGCGTGCGGGGAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACgtg                              <  2:2748628/99‑1 (MQ=255)
                                                                   tgtCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAAt                          <  2:1303807/100‑1 (MQ=255)
                                                                    gtCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATg                         >  2:1191926/1‑100 (MQ=255)
                                                                    gtCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATg                         >  1:3083437/1‑100 (MQ=255)
                                                                       cGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAag                      <  1:301158/100‑1 (MQ=255)
                                                                         gCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCACAGAGAATACGTGTAATGAagag                    >  1:2762367/1‑100 (MQ=255)
                                                                         gCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAagag                    <  2:459793/100‑1 (MQ=255)
                                                                          cGTGCGGGCAATTAACAGGGGGTCACTCATGAGACATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAgaga                   <  1:1834047/100‑1 (MQ=255)
                                                                           gTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaa                  <  2:2428273/100‑1 (MQ=255)
                                                                            tGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaaa                 <  2:3270552/100‑1 (MQ=255)
                                                                              cGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAAt               >  2:2417684/1‑100 (MQ=255)
                                                                                ggCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATaa             <  2:2642195/100‑1 (MQ=255)
                                                                                ggCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATaa             >  2:2394498/1‑100 (MQ=255)
                                                                                ggCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATaa             >  1:2102259/1‑100 (MQ=255)
                                                                                ggCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGAAATGAAGAGAAAATaa             >  1:974784/1‑100 (MQ=255)
                                                                                  cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATaa             >  1:3151001/1‑98 (MQ=255)
                                                                                  cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           >  2:3156121/1‑100 (MQ=255)
                                                                                  cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           <  1:1483971/100‑1 (MQ=255)
                                                                                    aTTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTa         >  1:2379900/1‑100 (MQ=255)
                                                                                     ttAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTaa        >  1:1155553/1‑100 (MQ=255)
                                                                                     ttAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTaa        >  1:13047/1‑100 (MQ=255)
                                                                                       aaCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAg      >  2:862965/1‑100 (MQ=255)
                                                                                       aaCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAg      >  2:1936006/1‑100 (MQ=255)
                                                                                       aaCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAAGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAg      <  1:2543870/100‑1 (MQ=255)
                                                                                        aCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGc     <  1:1115925/100‑1 (MQ=255)
                                                                                         cAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCg    <  1:1070153/100‑1 (MQ=255)
                                                                                         cAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCg    <  2:445146/100‑1 (MQ=255)
                                                                                           ggggGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCgtg  >  1:3258940/1‑100 (MQ=255)
                                                                                            ggggTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCgtg  >  1:2653611/1‑99 (MQ=255)
                                                                                               |                                                                                               
CCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTG  >  NC_000913/4489960‑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: ‑