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.4% 653.6 / ‑0.5 187E3D (GAA→GACyjgRDUF853 family protein with NTPase fold
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (94/90);  total (95/92)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

Only 100 of 187 total aligned reads displayed.

TACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTG  >  NC_000913/4489958‑4490150
                                                                                                 |                                                                                               
ttcccGTCCTCGTTGCTCCGGTAATTGACCCATGTCAGTTAGCCCTTCCCGGCAGTAAAAACAGTTCTGGGTCCTGCGTGCGGGCAATTAACAGGGGGTc                                                                                               <  2:3639059/98‑1 (MQ=255)
tACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTc                                                                                               <  1:100233/100‑1 (MQ=255)
 aCCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCa                                                                                              >  1:840234/1‑100 (MQ=255)
  cccGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCAc                                                                                             >  1:2125258/1‑100 (MQ=255)
   ccGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACt                                                                                            >  1:1925121/1‑100 (MQ=255)
    cGTCCCCGTTGCGCCGGTAATCAGCCCATGCCGGTTCGCCCTTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGTGGTCACTc                                                                                           <  2:2850324/100‑1 (MQ=255)
    cGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTc                                                                                           >  1:2891290/1‑100 (MQ=255)
        cccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATga                                                                                       >  1:1323101/1‑100 (MQ=255)
        cccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATga                                                                                       <  1:594036/100‑1 (MQ=255)
        cccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATga                                                                                       <  1:1977145/100‑1 (MQ=255)
        cccGGTGCGCCGGTATTCAGCCCATGGCGGTTCGCCGTTGCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGGCGGCTATTAACCGGAGGTCGTTCATg                                                                                        <  2:1897752/99‑1 (MQ=255)
         ccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                      >  1:1383862/1‑100 (MQ=255)
         ccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                      <  2:2738304/100‑1 (MQ=255)
         ccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                      <  1:3553377/100‑1 (MQ=255)
         ccGTTGCACCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgag                                                                                      >  1:3852448/1‑100 (MQ=255)
            ttGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGata                                                                                   <  1:336038/100‑1 (MQ=255)
             tGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGatat                                                                                  >  2:1136702/1‑100 (MQ=255)
             tGCGCCGGTAATCAGCCCATGACGGTTAGCCATCCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGata                                                                                   <  1:2049375/99‑1 (MQ=255)
              gcgcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtc                                                                                 <  1:1130669/100‑1 (MQ=255)
               cgcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcc                                                                                >  2:873045/1‑100 (MQ=255)
               cgcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcc                                                                                >  1:74935/1‑100 (MQ=255)
               cgcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAACTAACAGGGGGTCACTCATGAGATAtcc                                                                                <  1:3745755/100‑1 (MQ=255)
                gcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGCGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcct                                                                               >  2:1559157/1‑100 (MQ=255)
                  cGGTAATCAGCCCATGACGGGTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcctcc                                                                             <  2:1133050/100‑1 (MQ=255)
                   ggTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCa                                                                            <  1:2764920/100‑1 (MQ=255)
                   ggTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCa                                                                            >  2:632385/1‑100 (MQ=255)
                      aaTCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGtt                                                                         >  2:2127083/1‑100 (MQ=255)
                       aTCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGttt                                                                        >  1:106912/1‑100 (MQ=255)
                         cAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcc                                                                      <  1:3597368/100‑1 (MQ=255)
                         cAGCACATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGGCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcc                                                                      <  2:2275774/100‑1 (MQ=255)
                          agtccATGACGGTTAGCCATCCCCGGCAGTAAAAACACTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcc                                                                      <  2:2244414/96‑1 (MQ=255)
                             ccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcc                                                                  <  1:244319/100‑1 (MQ=255)
                             ccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcc                                                                  <  1:3563703/100‑1 (MQ=255)
                             ccATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcc                                                                  <  1:3322570/100‑1 (MQ=255)
                              cATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcct                                                                 >  2:3254610/1‑100 (MQ=255)
                               aTGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctgcctg                                                                <  1:1000944/100‑1 (MQ=255)
                                  aCGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAg                                                             <  2:1164805/100‑1 (MQ=255)
                                   cGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCCATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGt                                                            >  2:708968/1‑100 (MQ=255)
                                   cGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGt                                                            <  1:2065564/100‑1 (MQ=255)
                                   cGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGt                                                            >  1:2932410/1‑100 (MQ=255)
                                   cGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGt                                                            >  2:1248201/1‑100 (MQ=255)
                                    ggTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGta                                                           <  1:1506198/100‑1 (MQ=255)
                                       tAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAg                                                        <  2:3727254/100‑1 (MQ=255)
                                        aGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAgg                                                       <  2:2455483/100‑1 (MQ=255)
                                        aGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAgg                                                       <  2:179812/100‑1 (MQ=255)
                                        aGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAgg                                                       <  1:3941005/100‑1 (MQ=255)
                                        aGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAgg                                                       >  1:2724155/1‑100 (MQ=255)
                                        aGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAgg                                                       <  1:2039589/100‑1 (MQ=255)
                                        aGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAgg                                                       >  2:693111/1‑100 (MQ=255)
                                        aGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAgg                                                       <  2:3094059/100‑1 (MQ=255)
                                         gCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                      >  2:2317285/1‑100 (MQ=255)
                                         gCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                      >  2:959583/1‑100 (MQ=255)
                                           cATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGGGATATCCTCCAGTTTCCTGCCTGGAGTCTAGGCaa                                                    >  1:1791065/1‑100 (MQ=255)
                                           cATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCaa                                                    >  2:4105765/1‑100 (MQ=255)
                                               ttgggCAGTAAAAACGTTGCTGTGTCCGGCGTGCGGGCAATTAGCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTc                                                >  2:3874788/4‑100 (MQ=255)
                                               cccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCaa                                                    <  2:2378488/96‑1 (MQ=255)
                                               cccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTc                                                >  1:1137892/1‑100 (MQ=255)
                                               cccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTc                                                >  2:2325354/1‑100 (MQ=255)
                                               cccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGCCACTCATGAGATCTCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTc                                                >  2:1493573/1‑100 (MQ=255)
                                                ccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTcc                                               <  1:3880190/100‑1 (MQ=255)
                                                 cGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTGAAGAGCTATCCTCGAGTTTCCTGCCTGGAGTATAGGCAACGTTCg                                              >  2:2610778/1‑100 (MQ=255)
                                                 cGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCg                                              >  1:2672531/1‑100 (MQ=255)
                                                      gTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCaaa                                         >  1:481424/1‑100 (MQ=255)
                                                       tAAAAACAGACCTGTGTCCGGCGCGCGGGCAACAAACAAGGCGACACCCATGAGATATCCTCCAGCTCCCTGACTGGAGAACAGACAACCTTCGGCAAag                                        >  2:79592/1‑100 (MQ=255)
                                                        aaaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAaga                                       >  1:1283266/1‑100 (MQ=255)
                                                         aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                      >  1:1963938/1‑100 (MQ=255)
                                                          aaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAgaga                                     >  2:4008468/1‑100 (MQ=255)
                                                             cAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAAt                                   <  1:1248209/99‑1 (MQ=255)
                                                               gTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACg                                >  1:2597357/1‑100 (MQ=255)
                                                               gTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACg                                >  2:3754701/1‑100 (MQ=255)
                                                                  cTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACgtgt                             >  1:2674571/1‑100 (MQ=255)
                                                                      gtCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATg                         >  1:67383/1‑100 (MQ=255)
                                                                       tCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGa                        >  1:3146475/1‑100 (MQ=255)
                                                                        ccGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGaa                       >  2:3413086/1‑100 (MQ=255)
                                                                         cGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAag                      <  1:3799758/100‑1 (MQ=255)
                                                                         cGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAag                      <  2:840228/100‑1 (MQ=255)
                                                                          ggCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAaga                     >  1:3415434/1‑100 (MQ=255)
                                                                               gCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaaaa                <  1:98801/100‑1 (MQ=255)
                                                                                cGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAAt               >  1:3393566/1‑100 (MQ=255)
                                                                                 gggCACTTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATa              <  1:2251337/100‑1 (MQ=255)
                                                                                 gggCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATa              >  2:3466724/1‑100 (MQ=255)
                                                                                  ggCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATaa             <  2:289840/100‑1 (MQ=255)
                                                                                  ggCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATaa             >  1:1180972/1‑100 (MQ=255)
                                                                                   gCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCGACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAc            <  2:2070706/100‑1 (MQ=255)
                                                                                   gCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCGGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAc            <  2:218891/100‑1 (MQ=255)
                                                                                    cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCTGCAAAGAGAATACGATTAATGAAGAGaaaa                <  2:2033166/95‑1 (MQ=255)
                                                                                    cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           >  1:2100486/1‑100 (MQ=255)
                                                                                    cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           <  1:3140522/100‑1 (MQ=255)
                                                                                    cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           >  1:4124936/1‑100 (MQ=255)
                                                                                    cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           <  1:3103871/100‑1 (MQ=255)
                                                                                       ttAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTaa        >  2:1972001/1‑100 (MQ=255)
                                                                                        tAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTaaa       >  2:3159997/1‑100 (MQ=255)
                                                                                         aaCAGGGGGTGTCTCATGAGATATCCGCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAGCCTAAAg      <  2:559055/100‑1 (MQ=255)
                                                                                         aaCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAg      >  2:3641152/1‑100 (MQ=255)
                                                                                         aaCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAg      >  1:3613088/1‑100 (MQ=255)
                                                                                         aaCAGGGGGTCACTCATGAGATATCCTCCAGGTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTaa        >  2:1375769/1‑98 (MQ=255)
                                                                                           cAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCg    <  1:516513/100‑1 (MQ=255)
                                                                                           cAGGGGGTCCCTCATGAGATATCCTCCACTTTCCTGCCTGGAGTCGAGGCAACGGCCGGCCAGCAGAGCACGCGCAGTGAAGCGAAACGCACCTAAAGc     >  2:3622898/1‑99 (MQ=255)
                                                                                             ggggGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGCATACGTGTAATGAAGAGAAAATAACCTAAAGCgtg  >  2:3820357/1‑100 (MQ=255)
                                                                                             ggggGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCgtg  >  2:1858425/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 < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: