Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,358,838 G→A A10T (GCT→ACT)  gltF → GltF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,358,8380GA99.0% 342.2 / ‑4.9 97A10T (GCT→ACT) gltFGltF
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (95/1);  total (96/1)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

GCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAGCTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTAT  >  NC_000913/3358766‑3358927
                                                                        |                                                                                         
gCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACtt                                                                       >  1:3476183/1‑93 (MQ=255)
           taAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTgcagc                                                            >  1:1783331/1‑93 (MQ=255)
               aaataaatTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTc                                                        >  2:394029/1‑93 (MQ=255)
               aaataaatTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTc                                                        >  2:2045571/1‑93 (MQ=255)
                 ataaatTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAg                                                      >  2:2587212/1‑93 (MQ=255)
                 ataaatTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAg                                                      >  2:3311019/1‑93 (MQ=255)
                 ataaatTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAg                                                      >  2:2587213/1‑93 (MQ=255)
                  taaatTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGt                                                     >  1:88258/1‑93 (MQ=255)
                     atTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGcc                                                  >  1:3246047/1‑93 (MQ=255)
                     atTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGcc                                                  >  1:3374846/1‑93 (MQ=255)
                      tttGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCa                                                 >  1:781035/1‑93 (MQ=255)
                      tttGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCa                                                 >  1:1582215/1‑93 (MQ=255)
                      tttGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCa                                                 >  1:3552598/1‑93 (MQ=255)
                      tttGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCa                                                 >  2:4201711/1‑93 (MQ=255)
                       ttGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCAt                                                >  1:1755005/1‑93 (MQ=255)
                       ttGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCAt                                                >  2:2990793/1‑93 (MQ=255)
                         gTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATgg                                              >  1:4578652/1‑93 (MQ=255)
                         gTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATgg                                              >  2:2209008/1‑93 (MQ=255)
                          ttattaAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGc                                             >  2:2152863/1‑93 (MQ=255)
                           tattaAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCa                                            >  2:2187071/1‑93 (MQ=255)
                            attaAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCaa                                           >  1:4138240/1‑93 (MQ=255)
                            attaAGGATATGTTCATATGTTTTTCAAAAAGAACCTCAAAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCac                                           >  2:1140927/1‑92 (MQ=255)
                             ttaAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAAc                                          >  1:3615060/1‑93 (MQ=255)
                             ttaAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAAc                                          >  1:4444179/1‑93 (MQ=255)
                             ttaAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAAc                                          >  1:1666725/1‑93 (MQ=255)
                              taAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAGCTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACa                                         >  1:2352243/1‑93 (MQ=255)
                               aaGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAg                                        >  2:3679280/1‑93 (MQ=255)
                               aaGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAg                                        >  2:827211/1‑93 (MQ=255)
                               aaGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAg                                        >  2:827210/1‑93 (MQ=255)
                               aaGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAg                                        >  2:716315/1‑93 (MQ=255)
                                aGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGa                                       >  1:1484023/1‑93 (MQ=255)
                                aGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGa                                       >  1:1935227/1‑93 (MQ=255)
                                 ggATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGAt                                      >  1:2783114/1‑93 (MQ=255)
                                 ggATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGAt                                      >  2:2361635/1‑93 (MQ=255)
                                 ggATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGAt                                      >  2:4403760/1‑93 (MQ=255)
                                 ggATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGAt                                      >  2:4351667/1‑93 (MQ=255)
                                 ggATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGAt                                      >  2:293430/1‑93 (MQ=255)
                                  gATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGAtt                                     >  1:218960/1‑93 (MQ=255)
                                   atatGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTc                                    >  2:4621992/1‑93 (MQ=255)
                                   atatGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTc                                    >  2:3693716/1‑93 (MQ=255)
                                    tatGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCa                                   >  1:2093734/1‑93 (MQ=255)
                                    tatGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCa                                   >  2:4053343/1‑93 (MQ=255)
                                    tatGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCa                                   >  2:3107170/1‑93 (MQ=255)
                                    tatGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCa                                   >  1:850094/1‑93 (MQ=255)
                                    tatGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCa                                   >  2:1517555/1‑93 (MQ=255)
                                     atGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCaa                                  >  1:2428361/1‑93 (MQ=255)
                                      tGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAGCAGATTCAAc                                 >  1:1290462/1‑93 (MQ=255)
                                      tGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAAc                                 >  2:3495004/1‑93 (MQ=255)
                                      tGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAAc                                 >  2:3530262/1‑93 (MQ=255)
                                      tGTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAAc                                 >  1:2908741/1‑93 (MQ=255)
                                       gTTCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACt                                >  1:131388/1‑93 (MQ=255)
                                        ttCATATGTTTTTCAAAAAGAACCTCACAGCAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTg                               >  2:3361113/1‑93 (MQ=255)
                                         tCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGa                              >  2:1782851/1‑93 (MQ=255)
                                         tCATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGa                              >  2:1782720/1‑93 (MQ=255)
                                          cATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGAt                             >  1:1615519/1‑93 (MQ=255)
                                          cATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGAt                             >  1:2972329/1‑93 (MQ=255)
                                          cATATGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGCAAACAGATTAAACtaaa                             >  1:1615518/1‑90 (MQ=255)
                                           atatGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATa                            >  2:4506516/1‑93 (MQ=255)
                                           atatGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATa                            >  1:3430360/1‑93 (MQ=255)
                                            tatGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATAc                           >  2:1956942/1‑93 (MQ=255)
                                            tatGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATAc                           >  1:3161488/1‑93 (MQ=255)
                                            tatGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATAc                           >  1:4384187/1‑93 (MQ=255)
                                             atGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACa                          >  2:2064948/1‑93 (MQ=255)
                                              tGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAg                         >  2:2664403/1‑93 (MQ=255)
                                              tGTTTTTCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAg                         >  1:4459121/1‑93 (MQ=255)
                                                tttttCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGat                       >  1:1393801/1‑92 (MQ=255)
                                                tttttCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGaa                       >  2:2156199/1‑93 (MQ=255)
                                                 ttttCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCACCAGTTTAAACTGATAAAGAAt                      >  1:723484/1‑93 (MQ=255)
                                                  tttCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAAtt                     >  2:2485178/1‑93 (MQ=255)
                                                  tttCAAAAAGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAAtt                     >  2:717187/1‑93 (MQ=255)
                                                      aaaaaGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACt                 >  1:2347243/1‑93 (MQ=255)
                                                       aaaaGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTa                >  2:4324599/1‑93 (MQ=255)
                                                       aaaaGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTa                >  1:886446/1‑93 (MQ=255)
                                                        aaaGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTat               >  2:583146/1‑93 (MQ=255)
                                                        aaaGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTat               >  2:583145/1‑93 (MQ=255)
                                                         aaGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTatc              >  1:1914078/1‑93 (MQ=255)
                                                          aGAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTatca             >  2:3400710/1‑93 (MQ=255)
                                                           gAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTatcat            >  1:2242281/1‑93 (MQ=255)
                                                           gAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTatcat            >  2:1067865/1‑93 (MQ=255)
                                                           gAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTatcat            >  1:1732601/1‑93 (MQ=255)
                                                           gAACCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTatcat            >  2:2134625/1‑93 (MQ=255)
                                                            aaCCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTatcatc           <  1:3548871/93‑1 (MQ=255)
                                                             aCCTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCg          >  1:2021865/1‑93 (MQ=255)
                                                              ccTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCgg         >  1:4723989/1‑93 (MQ=255)
                                                               cTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGc        >  1:2551078/1‑93 (MQ=255)
                                                               cTCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGc        >  1:1655232/1‑93 (MQ=255)
                                                                tCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCg       >  2:4326951/1‑93 (MQ=255)
                                                                tCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCg       >  1:2894241/1‑93 (MQ=255)
                                                                tCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCg       >  2:4717423/1‑93 (MQ=255)
                                                                tCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCg       >  2:2896258/1‑93 (MQ=255)
                                                                tCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCg       >  2:1192302/1‑93 (MQ=255)
                                                                tCACAACAACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCg       >  2:1897710/1‑93 (MQ=255)
                                                                  acaacaACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAg     >  1:3202088/1‑93 (MQ=255)
                                                                  acaacaACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAg     >  2:1946609/1‑93 (MQ=255)
                                                                   caacaaCTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGt    >  2:370219/1‑93 (MQ=255)
                                                                   caacaaCTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGt    >  1:3701786/1‑93 (MQ=255)
                                                                     acaACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGtat  >  2:4015110/1‑93 (MQ=255)
                                                                     acaACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGtat  >  2:336132/1‑93 (MQ=255)
                                                                     acaACTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGtat  >  2:2790781/1‑93 (MQ=255)
                                                                        |                                                                                         
GCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAACCTCACAACAGCTGCTATTTGTGCAGCACTTTCTGTTGCAGCATTCAGTGCCATGGCAACAGATTCAACTGATACAGAATTAACTATCATCGGCGAGTAT  >  NC_000913/3358766‑3358927

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