Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 317,571 T→C 18.0% intergenic (+4/+155) ykgQ → / ← rclC pseudogene, putative dehydrogenase/reactive chlorine species (RCS) stress resistance inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913317,5710TC18.0% 412.5 / 79.2 188intergenic (+4/+155)ykgQ/rclCpseudogene, putative dehydrogenase/reactive chlorine species (RCS) stress resistance inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (62/91);  new base C (16/18);  total (79/109)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 192 total aligned reads displayed.

TTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATC  >  NC_000913/317473‑317670
                                                                                                  |                                                                                                   
   gTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATc                                                                                                 <  1:3136259/100‑1 (MQ=255)
     tCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCAt                                                                                               >  1:113065/1‑100 (MQ=255)
       gCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGt                                                                                             >  2:1564257/1‑100 (MQ=255)
       gCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGt                                                                                             >  2:2489788/1‑100 (MQ=255)
        cTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTa                                                                                            >  2:2423189/1‑100 (MQ=255)
        cTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTa                                                                                            <  1:3177496/100‑1 (MQ=255)
         ttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAt                                                                                           >  1:1785761/1‑100 (MQ=255)
          tGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtt                                                                                          <  2:3412060/100‑1 (MQ=255)
           gATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAttt                                                                                         >  1:1726272/1‑100 (MQ=255)
            aTTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtttt                                                                                        <  1:2673684/100‑1 (MQ=255)
            aTTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtttt                                                                                        <  2:3091671/100‑1 (MQ=255)
            aTTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttt                                                                                        <  2:3861919/100‑1 (MQ=255)
             tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTa                                                                                       <  1:1397161/100‑1 (MQ=255)
             tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTa                                                                                       >  1:26242/1‑100 (MQ=255)
             tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTa                                                                                       <  1:1964063/100‑1 (MQ=255)
              ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTa                                                                                       <  2:2823579/99‑1 (MQ=255)
              ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAc                                                                                      <  1:2293326/100‑1 (MQ=255)
              ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAc                                                                                      <  1:2143411/100‑1 (MQ=255)
              ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAc                                                                                      >  1:2135720/1‑100 (MQ=255)
               tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttac                                                                                     <  2:4057700/100‑4 (MQ=255)
               tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATACTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttac                                                                                     <  2:402407/100‑4 (MQ=38)
                gggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttaca                                                                                    <  1:3751056/100‑5 (MQ=38)
                   ttATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttacagca                                                                                 <  2:1677842/100‑8 (MQ=37)
                     aTGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgac                                                                               >  2:29671/1‑100 (MQ=255)
                     aTGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgac                                                                               >  1:2780574/1‑100 (MQ=255)
                      tGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacg                                                                              <  1:373535/100‑1 (MQ=255)
                      tGTTCGGGAGTCGCTACAACACGTTAATCCTGCATTGGTTCGTGGGGGTCAGGATGTAGAAATACATGATGGGTATTAATGATGTATTTTACCTaa                                                                                  >  2:4114342/1‑96 (MQ=255)
                         tGGGGAGTCGGTAAAACATTTTAATCCTGCATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATAAGGTATTTTACCTAAGACGACa                                                                           <  2:1882951/100‑1 (MQ=255)
                         tGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACa                                                                           <  2:419514/100‑1 (MQ=255)
                          ggggAGTCGGTAAAACATTTTAATCCTTAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACaa                                                                          <  2:2909858/100‑1 (MQ=255)
                          ggggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACaa                                                                          <  2:1249279/100‑1 (MQ=255)
                          ggggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACaa                                                                          >  2:2872170/1‑100 (MQ=255)
                              aGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAc                                                                      <  2:73264/100‑1 (MQ=255)
                               gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                                     <  1:3388649/100‑1 (MQ=255)
                                  ggTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCttt                                                                  >  1:746466/1‑100 (MQ=255)
                                   gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtttt                                                                 <  2:2700249/100‑1 (MQ=255)
                                   gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtttt                                                                 <  1:853299/100‑1 (MQ=255)
                                   gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAAAAAAAAAAAAAAAAACCttt                                                                  <  2:116044/99‑1 (MQ=255)
                                    tAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTa                                                                <  1:3564829/100‑1 (MQ=255)
                                      aaaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTaa                                                               >  1:4051292/1‑99 (MQ=255)
                                      aaaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAAc                                                              >  2:4073588/1‑100 (MQ=255)
                                      aaaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAACCATGTATTTTACCTAAGACGACAAAAACCTTTTAAc                                                              <  1:2957673/100‑1 (MQ=255)
                                       aaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACa                                                             >  1:1123766/1‑100 (MQ=255)
                                          aTTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAggg                                                          >  1:211716/1‑100 (MQ=255)
                                            tttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAg                                                        <  2:1381588/100‑1 (MQ=255)
                                             ttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGt                                                       >  1:3296011/1‑100 (MQ=255)
                                             ttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGt                                                       <  1:3208310/100‑1 (MQ=255)
                                                   cTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATAtt                                                 >  1:276633/1‑100 (MQ=255)
                                                       aTTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATc                                             <  2:2216148/100‑1 (MQ=255)
                                                       aTTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATc                                             >  1:2388442/1‑100 (MQ=255)
                                                       aTTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATc                                             <  2:2150639/100‑1 (MQ=255)
                                                        tttGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCt                                            <  2:211716/100‑1 (MQ=255)
                                                         ttGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCtt                                           <  2:4148876/100‑1 (MQ=255)
                                                          tGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTc                                          <  2:573860/100‑1 (MQ=255)
                                                            ttCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCAc                                        <  1:1511752/100‑1 (MQ=255)
                                                             tCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCAct                                       >  1:3366600/1‑100 (MQ=255)
                                                              cGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActc                                      <  1:2661402/100‑1 (MQ=255)
                                                              cGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActc                                      >  2:864069/1‑100 (MQ=255)
                                                               gTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActct                                     <  1:2813727/100‑1 (MQ=255)
                                                                tGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTg                                    >  1:1021110/1‑100 (MQ=255)
                                                                 ggTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGt                                   >  1:999596/1‑100 (MQ=255)
                                                                  gTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTc                                  <  2:922859/100‑1 (MQ=255)
                                                                   ttGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCa                                 <  1:3630786/100‑1 (MQ=255)
                                                                   ttGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCa                                 <  1:3801209/100‑1 (MQ=255)
                                                                    tGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCat                                <  1:728182/100‑1 (MQ=255)
                                                                       gtcttaggtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtc                             >  1:278049/10‑100 (MQ=25)
                                                                          gCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTcc                          <  1:1827909/100‑1 (MQ=255)
                                                                          gCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTcc                          <  2:2325397/100‑1 (MQ=255)
                                                                          gCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTcc                          <  2:2030811/100‑1 (MQ=255)
                                                                           cTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCg                         >  1:3920170/1‑100 (MQ=255)
                                                                           cTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTAACTCTGTCATAtctac                          >  2:3915990/1‑97 (MQ=255)
                                                                            tGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCgg                        >  2:983931/1‑100 (MQ=255)
                                                                            tGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCgg                        <  2:3617346/100‑1 (MQ=255)
                                                                            tGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCgg                        >  1:2463440/1‑100 (MQ=255)
                                                                             gTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                       >  2:2692586/1‑100 (MQ=255)
                                                                             gTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                       >  2:768535/1‑100 (MQ=255)
                                                                              gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                       >  1:1266988/3‑99 (MQ=255)
                                                                              gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                      >  1:4169300/3‑100 (MQ=255)
                                                                              gtaaaaTACATGATTAATATCAATCATGTATTTTACCAAAGACGACAAAAACCTTTTAACAGGGTGTGATATTGATCTTCACTCTGTCATATCTCCg                         >  2:2248830/3‑97 (MQ=37)
                                                                               taaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTaa                     >  2:3469446/2‑100 (MQ=255)
                                                                               aaaaaTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTaa                     >  2:3584667/1‑100 (MQ=255)
                                                                                aaaaTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAATGATATTGATCTTTACTCTCTCATATCTCGGGTAat                    >  2:2917035/1‑100 (MQ=255)
                                                                                aaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAat                    >  1:1298576/1‑100 (MQ=255)
                                                                                aaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAat                    >  1:1322453/1‑100 (MQ=255)
                                                                                aaaaTACATGATAGATATTAATCATGAATTTAACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAat                    <  1:1184031/100‑1 (MQ=255)
                                                                                  aaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAatat                  >  2:1771267/1‑100 (MQ=255)
                                                                                  aaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAatat                  >  2:3322879/1‑100 (MQ=255)
                                                                                   aTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATG                 <  1:2512254‑M1/100‑1 (MQ=255)
                                                                                   aTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATG                 >  1:2987223‑M1/1‑100 (MQ=255)
                                                                                        tgattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTC            <  1:290440‑M1/100‑1 (MQ=255)
                                                                                         gattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGTGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGT             >  2:2954454‑M1/1‑98 (MQ=255)
                                                                                         gattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCA           <  2:1785757‑M1/100‑1 (MQ=255)
                                                                                           ttgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGG         <  2:446457‑M1/100‑1 (MQ=255)
                                                                                            tgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGC        <  2:1942282‑M1/100‑1 (MQ=255)
                                                                                               tatTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtc     <  1:1117503/100‑1 (MQ=255)
                                                                                               tatTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtc     <  2:656637/100‑1 (MQ=255)
                                                                                                atTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtca    <  2:1311746/100‑1 (MQ=255)
                                                                                                 ttAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtcat   >  2:2551495/1‑100 (MQ=255)
                                                                                                  tAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtcatc  <  1:3481517/100‑1 (MQ=255)
                                                                                                  tAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtcat   >  2:2652162/1‑99 (MQ=255)
                                                                                                  |                                                                                                   
TTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATC  >  NC_000913/317473‑317670

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