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,0550TG98.7% 388.5 / ‑4.1 115E3D (GAA→GACyjgRDUF853 family protein with NTPase fold
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (57/56);  total (57/58)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 116 total aligned reads displayed.

ACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTGT  >  NC_000913/4489959‑4490151
                                                                                                |                                                                                                
aCCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCa                                                                                               >  2:2734137/1‑100 (MQ=255)
 cccGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCAc                                                                                              >  1:1641972/1‑100 (MQ=255)
 cccGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCAc                                                                                              >  1:1718895/1‑100 (MQ=255)
    gTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCa                                                                                           <  2:2704822/100‑1 (MQ=255)
      ccccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCTGGCAATTAACAGGGGGTCACTCATg                                                                                         <  2:692876/100‑1 (MQ=255)
      ccccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATg                                                                                         <  2:1377090/100‑1 (MQ=255)
      ccccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATg                                                                                         <  1:3016363/100‑1 (MQ=255)
      caccGTTGCGCCGGTAATCAGTCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATg                                                                                         <  1:2624369/98‑1 (MQ=255)
       cccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATga                                                                                        <  2:460168/100‑1 (MQ=255)
        ccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                       <  2:1054656/100‑1 (MQ=255)
        ccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGAGCAATTAACAGGGGGTCACTCATaag                                                                                       >  1:640835/1‑100 (MQ=255)
          gTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGat                                                                                     >  1:832198/1‑100 (MQ=255)
           ttGCGTCTGTAATCAGCCCTTGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCTGGCAATTTACAGGGGGTCACTCATGTGat                                                                                     <  2:222627/99‑1 (MQ=255)
           ttGCGCCTTTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgtg                                                                                       <  1:2055427/97‑3 (MQ=255)
              cgcCGGTAATCAGCTCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcc                                                                                 <  1:1340541/100‑1 (MQ=255)
              cgcCGGTAATCAGCCCATGACGGTTAGCCATTGCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcc                                                                                 >  1:2889123/1‑100 (MQ=255)
              cgcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcc                                                                                 >  1:1160111/1‑100 (MQ=255)
               gcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAGAAGTACTGTGTCCGGCGGGCGGGCAATTAACGGCGGGTCACTGGTGGGATAtcc                                                                                 >  2:2490048/1‑99 (MQ=255)
               gcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGACACTCATGAGATAtcct                                                                                >  2:708627/1‑100 (MQ=255)
                 cGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGata                                                                                    >  2:2950314/1‑94 (MQ=255)
                        cAGCCCATGACGGTTAGCTATTCCCGGTAGTAAAAATATTTTTGTTTCCTGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcc                                                                       <  1:68593/100‑1 (MQ=255)
                        cAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTcc                                                                       <  1:1428506/100‑1 (MQ=255)
                        cAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcc                                                                       <  2:1678407/100‑1 (MQ=255)
                        cAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcc                                                                       >  1:2235857/1‑100 (MQ=255)
                        cAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcc                                                                       <  1:2008353/100‑1 (MQ=255)
                        cAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGGCACTCATGAGATATCCTCCAGTTTcc                                                                       >  2:2509044/1‑100 (MQ=255)
                         aGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcct                                                                      >  1:2118576/1‑100 (MQ=255)
                         aGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGAAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcct                                                                      >  1:1795365/1‑100 (MQ=255)
                          gCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctg                                                                     <  2:2722968/100‑1 (MQ=255)
                          gCCCATGACGGTAAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACGCATGAGATATCCTCCAGTTTcctg                                                                     >  2:3048150/1‑100 (MQ=255)
                           cccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgc                                                                    <  2:2466482/100‑1 (MQ=255)
                            ccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcc                                                                   <  2:640831/100‑1 (MQ=255)
                            ccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcc                                                                   <  2:660848/100‑1 (MQ=255)
                              aTGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcctg                                                                 <  1:1238859/100‑1 (MQ=255)
                              aTGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcctg                                                                 <  2:1549873/100‑1 (MQ=255)
                                 aCGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAg                                                              <  1:883952/100‑1 (MQ=255)
                                 aCGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAg                                                              <  1:2790271/100‑1 (MQ=255)
                                 aCGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAg                                                              <  1:2790236/100‑1 (MQ=255)
                                 aCGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAg                                                              >  2:2489546/1‑100 (MQ=255)
                                   ggTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGta                                                            >  1:2696483/1‑100 (MQ=255)
                                    gTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGta                                                            <  1:883858/99‑1 (MQ=255)
                                     ttAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGtata                                                          >  1:665012/1‑100 (MQ=255)
                                      tAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGTAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAg                                                         <  2:1503655/100‑1 (MQ=255)
                                      tAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAg                                                         <  1:1634543/100‑1 (MQ=255)
                                      tAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAg                                                         <  1:2357508/100‑1 (MQ=255)
                                      tAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAg                                                         <  1:2159377/100‑1 (MQ=255)
                                      tAGCCATTCCCGACAGTAAAAACAGTACTGTGTCCGGCGTGCGGGCAATTAACAGGGAGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAg                                                         >  2:413685/1‑100 (MQ=255)
                                        gCCATTCCCGGCAGTAAAAACAGTTCTGTGTGCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                       >  2:1178435/1‑100 (MQ=255)
                                        gCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                       <  1:1193650/100‑1 (MQ=255)
                                        gCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                       <  2:324086/100‑1 (MQ=255)
                                         ccATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCa                                                      >  1:1625865/1‑100 (MQ=255)
                                         ccATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCa                                                      >  1:2403448/1‑100 (MQ=255)
                                            ttCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCGCCAGTTTCCTGCCTGGAGTATAGGCAACg                                                   <  1:2914783/100‑1 (MQ=255)
                                             tcgcGGCAGTAAAGACTGTTCTGTGGCCGGCGTGCGGGCAATTAACAGGGGGGCACTCAGGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCa                                                      <  2:2982698/93‑1 (MQ=255)
                                              cccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTc                                                 >  2:392282/1‑100 (MQ=255)
                                              accGGCAGTAAAAACAGTCCTATGTCCCGCATGCGGGCAAAAAACACGGGGTCACTCATGAGATATCCTCCATTTTCCTACATGGAGtat                                                           >  2:2303930/2‑90 (MQ=255)
                                               ccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTcc                                                >  2:1452778/1‑100 (MQ=255)
                                                 ggtaGTAAAAACAGGTCTGTGTCGAGCGTGCGTGCAATTAACAGGGGTTCACTCATGCGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCgg                                              <  2:2943190/97‑1 (MQ=255)
                                                  gCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGc                                             <  2:2502505/100‑1 (MQ=255)
                                                  gCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGc                                             >  1:1978702/1‑100 (MQ=255)
                                                     gTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCaaa                                          >  1:1942759/1‑100 (MQ=255)
                                                     gTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCaaa                                          <  1:708631/100‑1 (MQ=255)
                                                        aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGTGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                       <  1:2525735/100‑1 (MQ=255)
                                                        aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                       <  1:2096620/100‑1 (MQ=255)
                                                        aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                       <  2:2387789/100‑1 (MQ=255)
                                                        aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                       <  2:2387809/100‑1 (MQ=255)
                                                        aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                       <  1:2489575/100‑1 (MQ=255)
                                                        aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                       <  1:60787/100‑1 (MQ=255)
                                                         aaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAgaga                                      >  1:2628316/1‑100 (MQ=255)
                                                          aaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGaa                                     >  2:2481487/1‑100 (MQ=255)
                                                            cAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATa                                   <  1:2851709/100‑1 (MQ=255)
                                                            cAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATa                                   <  2:2591242/100‑1 (MQ=255)
                                                              gTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTACTGCAAAGAGAATa                                   >  2:540690/1‑98 (MQ=255)
                                                              atcttgtgTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACg                                 <  1:879138/96‑1 (MQ=255)
                                                                tCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACgtg                               <  2:255019/100‑1 (MQ=255)
                                                                 atgtgTCCGTCGTGCGGTCAATTAACAGGGAGTCACTCATGAGATATCCACCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACgt                                >  2:670132/2‑98 (MQ=255)
                                                                  tgtgTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATGCgtgt                              >  2:1849531/1‑99 (MQ=255)
                                                                  tgtgTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACgaga                              >  2:879132/1‑96 (MQ=255)
                                                                  tgtgTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTa                             >  1:959521/1‑100 (MQ=255)
                                                                  tgtgTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTa                             >  1:2035650/1‑100 (MQ=255)
                                                                       ccGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGAATAGGCAACGTACGGCAAAGAGAATACGTGTAATGaa                        >  1:1946075/1‑100 (MQ=255)
                                                                        cGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATTTCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAag                       <  1:1011835/100‑1 (MQ=255)
                                                                        cGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAag                       >  1:2304814/1‑100 (MQ=255)
                                                                             tGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaaa                  >  2:178595/1‑100 (MQ=255)
                                                                              gCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaaaa                 >  2:1995612/1‑100 (MQ=255)
                                                                              gCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaaaa                 >  2:3027873/1‑100 (MQ=255)
                                                                               cGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAAt                >  2:2086393/1‑100 (MQ=255)
                                                                                gggCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATa               >  2:2337199/1‑100 (MQ=255)
                                                                                 ggCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATaa              >  2:1907089/1‑100 (MQ=255)
                                                                                  gCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAc             >  1:1428625/1‑100 (MQ=255)
                                                                                   cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc            <  1:2155624/100‑1 (MQ=255)
                                                                                   cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc            <  2:2333313/100‑1 (MQ=255)
                                                                                      ttAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTaa         <  2:1718875/100‑1 (MQ=255)
                                                                                        aaCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAg       >  2:582038/1‑100 (MQ=255)
                                                                                        aaCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAg       >  2:1553547/1‑100 (MQ=255)
                                                                                        aaCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAg       >  2:2782386/1‑100 (MQ=255)
                                                                                         aCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCAAGAGTATAGGCAAAATCCGCCAAAGAAAATACGGGTAATGAAGAGAAAAAAAAATaaa        >  1:2305315/1‑98 (MQ=255)
                                                                                          cAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGCGAAAATAACCTAAAGCg     >  2:396259/1‑100 (MQ=255)
                                                                                          cAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGAAAAGTACGGAAAAGAGAATACGTGTAATGAAGAGAAAATAAAATAAAg       >  1:2635369/1‑98 (MQ=255)
                                                                                             ggggTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCgtgt  >  2:2867307/1‑100 (MQ=255)
                                                                                                |                                                                                                
ACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTGT  >  NC_000913/4489959‑4490151

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