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,0550TG97.4% 532.2 / 2.5 156E3D (GAA→GACyjgRDUF853 family protein with NTPase fold
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base G (65/87);  total (68/88)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

Only 100 of 156 total aligned reads displayed.

TACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTG  >  NC_000913/4489958‑4490150
                                                                                                 |                                                                                               
tACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTc                                                                                               >  2:1504091/1‑100 (MQ=255)
    cttCCCCGTTTCGCCGGTAATCACCCCATGACCGTTAGCCATTTCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCAc                                                                                             <  2:2797173/96‑1 (MQ=255)
    cGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTc                                                                                           <  2:1257860/100‑1 (MQ=255)
    cGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTc                                                                                           >  2:2662500/1‑100 (MQ=255)
     gTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCa                                                                                          <  1:1251021/100‑1 (MQ=255)
     gTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTATGTCCGGCGGGCGGGCAATTAACAGGGGGTCACTCa                                                                                          >  2:2374396/1‑100 (MQ=255)
      tCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCAt                                                                                         <  2:2145049/100‑1 (MQ=255)
       ccccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATg                                                                                        >  1:2250882/1‑100 (MQ=255)
       ccccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATg                                                                                        <  2:128490/100‑1 (MQ=255)
       ccccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATg                                                                                        >  2:2791819/1‑100 (MQ=255)
       ccccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATg                                                                                        <  2:1703413/100‑1 (MQ=255)
        cccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATga                                                                                       <  2:2686924/100‑1 (MQ=255)
        cccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATga                                                                                       >  1:1826300/1‑100 (MQ=255)
         ctgTTGCGCCGGTAATCAGCCCATGAAGGTTAGCCATTCCCGGCAGTAAAAAAAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATg                                                                                        >  1:306906/3‑98 (MQ=255)
          cGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgaga                                                                                     >  2:233248/1‑100 (MQ=255)
          cGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgaga                                                                                     >  1:428457/1‑100 (MQ=255)
            ttGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGata                                                                                   <  2:2612061/100‑1 (MQ=255)
            ttGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGata                                                                                   <  1:897993/100‑1 (MQ=255)
             tGCGCCGGTAATCAGCCCATGACGGTGAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGatat                                                                                  <  1:350043/100‑1 (MQ=255)
               cgcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcc                                                                                >  2:3134863/1‑100 (MQ=255)
                gcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcct                                                                               >  1:497341/1‑100 (MQ=255)
                 ccGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcctc                                                                              >  2:1432724/1‑100 (MQ=255)
                 ccGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcctc                                                                              >  2:2744154/1‑100 (MQ=255)
                    gTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAg                                                                           <  1:1885755/100‑1 (MQ=255)
                       aTCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGttt                                                                        >  2:1734935/1‑100 (MQ=255)
                       aTCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGttt                                                                        >  1:689902/1‑100 (MQ=255)
                       aTCAGCCCATGACGGTTAGCCATTCCAGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATTCTCCAGtt                                                                         >  1:1242550/1‑99 (MQ=255)
                        tCAGCCCATGTCGGTTCGCCATTCCCGGCAGTAAAAGCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGt                                                                          <  1:1164805/97‑1 (MQ=255)
                        tCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTc                                                                       >  2:1176801/1‑100 (MQ=255)
                          aGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGCCACTCATGAGATATCCTCCAGTTTcct                                                                     >  1:3108349/1‑100 (MQ=255)
                            cccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgc                                                                   <  2:62020/100‑1 (MQ=255)
                            cccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgc                                                                   >  1:2125266/1‑100 (MQ=255)
                             ccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcc                                                                  <  1:2737154/100‑1 (MQ=255)
                             ccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcc                                                                  >  2:1686752/1‑100 (MQ=255)
                              cATGACGGTTAGCCATGCCCGGCAGTAAAAACGGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctg                                                                    <  1:1556946/97‑1 (MQ=255)
                               aTGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcctg                                                                <  1:2570813/100‑1 (MQ=255)
                               aTGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcctg                                                                <  2:1937716/100‑1 (MQ=255)
                               aTGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcctg                                                                <  2:1042485/100‑1 (MQ=255)
                                tGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcctg                                                                <  1:2668369/99‑1 (MQ=255)
                                  aCGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAg                                                             <  2:1735517/100‑1 (MQ=255)
                                   cGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGt                                                            >  2:121969/1‑100 (MQ=255)
                                   cGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGt                                                            >  2:1046460/1‑100 (MQ=255)
                                      ttAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGtata                                                         <  2:1853924/100‑1 (MQ=255)
                                       tAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAg                                                        <  1:2019243/100‑1 (MQ=255)
                                       tAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAg                                                        <  1:3295127/100‑1 (MQ=255)
                                        aGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAgg                                                       >  1:2077249/1‑100 (MQ=255)
                                        aGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAgg                                                       <  1:3078602/100‑1 (MQ=255)
                                         gCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                      <  1:714696/100‑1 (MQ=255)
                                          ccATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCa                                                     <  1:930044/100‑1 (MQ=255)
                                          ccATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCa                                                     <  2:1077995/100‑1 (MQ=255)
                                          ccATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCa                                                     <  1:930058/100‑1 (MQ=255)
                                           cATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCaa                                                    >  1:2542256/1‑100 (MQ=255)
                                             ttCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACg                                                  <  2:1632371/100‑1 (MQ=255)
                                                ccGGCAGTAAAAACAGTTCTGTGTCCGGCGCGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGCATATGCAACGTcc                                               >  1:431942/1‑100 (MQ=255)
                                                 cGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCg                                              >  2:2039609/1‑100 (MQ=255)
                                                  ggCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCgg                                             <  1:3314455/100‑1 (MQ=255)
                                                  ggCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCgg                                             >  1:2433163/1‑100 (MQ=255)
                                                  gacAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCgg                                             >  2:2285640/3‑100 (MQ=255)
                                                   gCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGc                                            <  1:1322653/100‑1 (MQ=255)
                                                     aGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCa                                           >  1:645785/1‑99 (MQ=255)
                                                      gTAAAAAGAGTTCTGTGTCCGGCGTGCGGTCACTTAACAGGGGGCCACACAAGAGATATCCTCCAGTTTCCTGCCTGGAGTATCGGCAACTGCCGGCa                                           >  1:23153/1‑98 (MQ=255)
                                                      gTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGc                                            <  1:577374/97‑1 (MQ=255)
                                                      gTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCaaa                                         <  1:1034095/100‑1 (MQ=255)
                                                      gTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCaaa                                         >  2:1526249/1‑100 (MQ=255)
                                                       tAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAag                                        <  2:2318502/100‑1 (MQ=255)
                                                         aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                      <  2:1405431/100‑1 (MQ=255)
                                                         aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                      <  1:513040/100‑1 (MQ=255)
                                                         aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                      <  1:1736043/100‑1 (MQ=255)
                                                          aaccAGTTCTGTGTCCGGCGTGCGGACAATTAACAGGGGGTCACTCATGAGATAGCCTCCAGTTGCCTGCCTGGAGTATAGACAACGTCCGGCAAagag                                      <  1:566630/96‑1 (MQ=255)
                                                          aaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAgaga                                     <  2:3357883/100‑1 (MQ=255)
                                                               gTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACg                                <  1:2703086/100‑1 (MQ=255)
                                                               gTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACg                                >  1:2481222/1‑100 (MQ=255)
                                                               gTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACg                                >  1:1420033/1‑100 (MQ=255)
                                                               gTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGACTGGAGTATAGGCAACGTCAGGCAAAGAGAATatg                                >  1:2430536/1‑98 (MQ=255)
                                                                 tCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACgtg                              <  2:30127/100‑1 (MQ=255)
                                                                 tCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACgtg                              <  2:243731/100‑1 (MQ=255)
                                                                   tgtgTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTa                            >  2:2073186/1‑100 (MQ=255)
                                                                   tgtgTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTa                            >  1:1543101/1‑100 (MQ=255)
                                                                     tgtCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAAt                          <  2:2863836/100‑1 (MQ=255)
                                                                        ccGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGaa                       <  1:3032653/100‑1 (MQ=255)
                                                                        ccGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGaa                       >  1:1150657/1‑100 (MQ=255)
                                                                         cGGCGTGCGGGCAATTAACAGGTGGTCACTCCTGAGACATCCTCCAGCTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAag                      <  1:2285624/100‑1 (MQ=255)
                                                                         cGGCGTGCGGGCAATTAACAGGGGGTCACCCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGTCAACGTCCGGCAAAGAGAATACGTGTAATGAag                      <  1:2036926/100‑1 (MQ=255)
                                                                         cGGCGTGCGGGCAATTAACAGGGGGTCACACATGAGATATCCTCCCGTTTCCTGCCTGTACTATCGGCAACGTCCGGCAAAGAGTATATGCGTaa                           >  1:1382247/1‑95 (MQ=255)
                                                                            cGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAgaga                   <  1:2343226/100‑1 (MQ=255)
                                                                             gTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaa                  <  1:1085352/100‑1 (MQ=255)
                                                                             gTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaa                  <  2:3118166/100‑1 (MQ=255)
                                                                              tGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaaa                 <  2:1644594/100‑1 (MQ=255)
                                                                                cGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGACAACGTCCGCCAAAGAGAATACGTGTAATGAAGAGaaaa                >  1:1647193/1‑99 (MQ=255)
                                                                                   gCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAc            <  2:3071515/100‑1 (MQ=255)
                                                                                    cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           >  1:3026239/1‑100 (MQ=255)
                                                                                    cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           <  1:88091/100‑1 (MQ=255)
                                                                                         aaCAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAg      >  2:453221/1‑100 (MQ=255)
                                                                                         aaCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAATGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAg      >  2:1617037/1‑100 (MQ=255)
                                                                                         aaCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAg      <  1:3278386/100‑1 (MQ=255)
                                                                                           cAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCg    <  1:2711999/100‑1 (MQ=255)
                                                                                            aGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCgt   <  2:2498248/100‑1 (MQ=255)
                                                                                             ggggTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCgtg  >  1:1410413/1‑100 (MQ=255)
                                                                                             ggggGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCgtg  >  2:303930/1‑100 (MQ=255)
                                                                                             ggggGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCgtg  >  1:1324861/1‑100 (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: