Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,517,027 C→T A231T (GCA→ACA)  fecR ← anti‑sigma transmembrane signal transducer for ferric citrate transport; periplasmic FecA‑bound ferric citrate sensor and cytoplasmic FecI ECF sigma factor activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,517,0270CT99.0% 717.7 / ‑2.9 208A231T (GCA→ACA) fecRanti‑sigma transmembrane signal transducer for ferric citrate transport; periplasmic FecA‑bound ferric citrate sensor and cytoplasmic FecI ECF sigma factor activator
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (93/112);  total (93/115)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

Only 100 of 210 total aligned reads displayed.

AACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCT  >  NC_000913/4516928‑4517122
                                                                                                   |                                                                                               
   gggTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGcc                                                                                              >  1:3573284/1‑100 (MQ=255)
    ggTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCa                                                                                             >  1:3410834/1‑100 (MQ=255)
     gTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCaa                                                                                            >  1:3089266/1‑100 (MQ=255)
       tAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAAc                                                                                          <  2:3274052/100‑1 (MQ=255)
       tAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAAc                                                                                          <  2:1380806/100‑1 (MQ=255)
       tAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAAc                                                                                          <  1:4600731/100‑1 (MQ=255)
       tAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAAc                                                                                          <  2:4158219/100‑1 (MQ=255)
       tAGCGTGGCTATCACCTCACCCAGCGGTTTATAGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAAc                                                                                          <  2:1565104/100‑1 (MQ=255)
        aGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACt                                                                                         >  1:3506789/1‑100 (MQ=255)
        aGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACt                                                                                         <  1:1666211/100‑1 (MQ=255)
         gCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTc                                                                                        >  2:2332386/1‑100 (MQ=255)
          cGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCa                                                                                       >  1:597905/1‑100 (MQ=255)
          cGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCa                                                                                       <  1:3868852/100‑1 (MQ=255)
          cGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCa                                                                                       <  2:1089767/100‑1 (MQ=255)
           gTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCag                                                                                      <  2:4234293/100‑1 (MQ=255)
            tGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGACGTCCTTCGTCCAGATTTTACTCTCGCCATCCACCTGTTTCACTGTGCCATACTCaga                                                                                     <  2:2145964/100‑1 (MQ=255)
              gCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAgg                                                                                   <  1:689255/100‑1 (MQ=255)
              gCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTAAGGATGTCCCTCGTCCAGCTTGTACTCTCGTCATCCATCGGTCTCACTGTGCCAAACTCAGAgg                                                                                   >  2:3583638/1‑100 (MQ=255)
                tATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGgcg                                                                                 >  2:2573869/1‑100 (MQ=255)
                 aTCACCTCACCCAGCGGTTTATCGCTGAAGCTCACGATGTCCTTTGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGgcgc                                                                                <  2:5242101/100‑1 (MQ=255)
                   cacctcaccCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTg                                                                              <  1:148572/100‑1 (MQ=255)
                     cctcaccCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGaa                                                                            >  1:1054856/1‑100 (MQ=255)
                          cccAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCa                                                                       <  2:4667005/100‑1 (MQ=255)
                          cccAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCa                                                                       <  2:2835144/100‑1 (MQ=255)
                          cccAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCa                                                                       >  2:3980126/1‑100 (MQ=255)
                          cccAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCa                                                                       >  2:1588363/1‑100 (MQ=255)
                           ccAGCGGTTTATCGCTGAAGCTCAGGCTGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAg                                                                      >  2:1270793/1‑100 (MQ=255)
                           ccAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAg                                                                      <  2:1170534/100‑1 (MQ=255)
                            cAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAgg                                                                     >  1:5493121/1‑100 (MQ=255)
                             aGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGc                                                                    >  2:1459968/1‑100 (MQ=255)
                              gCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCCTCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCt                                                                   <  1:554551/100‑1 (MQ=255)
                              gCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCt                                                                   >  2:4343966/1‑100 (MQ=255)
                              gCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCt                                                                   >  2:3186633/1‑100 (MQ=255)
                                 gTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTTTCGTCATCCAGCGGTTTCACTGTGCCTAACTCAGAGGCGCTGCACTCCCGGCTTTc                                                                <  2:4534634/100‑1 (MQ=255)
                                 gTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTGTCGTCATCCAGCGGTTTCACTGTGCCTAACTCAGAGGCGCTGTACTTCTGGCTTTc                                                                <  2:4535768/100‑1 (MQ=255)
                                  tttATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCa                                                               <  1:3948089/100‑1 (MQ=255)
                                   ttATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCAc                                                              <  1:3307993/100‑1 (MQ=255)
                                    tATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCAcc                                                             <  1:2932386/100‑1 (MQ=255)
                                    tATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCAcc                                                             <  1:3649317/100‑1 (MQ=255)
                                    tATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCACCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCAcc                                                             <  1:2118864/100‑1 (MQ=255)
                                     aTCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCa                                                            <  2:368744/100‑1 (MQ=255)
                                      tCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAg                                                           <  1:2450268/100‑1 (MQ=255)
                                      tCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAg                                                           <  2:5365587/100‑1 (MQ=255)
                                       cGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGc                                                          <  1:1078050/100‑1 (MQ=255)
                                        gCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCg                                                         <  1:4257456/100‑1 (MQ=255)
                                         cTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGt                                                        >  1:2091104/1‑100 (MQ=255)
                                           gAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGttn                                                      <  1:505757/100‑2 (MQ=255)
                                            aaGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCa                                                     <  2:5393512/100‑1 (MQ=255)
                                             aGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCAc                                                    <  2:723187/100‑1 (MQ=255)
                                             aGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCAc                                                    <  2:3540251/100‑1 (MQ=255)
                                              gCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACg                                                   <  1:781052/100‑1 (MQ=255)
                                               cTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGa                                                  <  1:2356188/100‑1 (MQ=255)
                                               cTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGa                                                  >  2:2084719/1‑100 (MQ=255)
                                                   ggATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCg                                              <  1:5194060/100‑1 (MQ=255)
                                                      tGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGttt                                           >  1:5500887/1‑100 (MQ=255)
                                                      tGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGttt                                           >  1:1962774/1‑100 (MQ=255)
                                                        tCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGAGttttt                                         >  2:1602563/1‑100 (MQ=255)
                                                         ccTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTg                                        >  1:4976360/1‑100 (MQ=255)
                                                               tCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCggg                                  <  1:866721/100‑1 (MQ=255)
                                                               tCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCggg                                  <  1:2826604/100‑1 (MQ=255)
                                                               tCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCggg                                  <  1:3362825/100‑1 (MQ=255)
                                                               tCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCggg                                  <  2:2233021/100‑1 (MQ=255)
                                                               tCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCggg                                  <  1:3932954/100‑1 (MQ=255)
                                                                ccAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAACCTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCgggg                                 >  2:1898378/1‑100 (MQ=255)
                                                                ccAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCgggg                                 >  2:4890221/1‑100 (MQ=255)
                                                                ccAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCgggg                                 <  1:4213635/100‑1 (MQ=255)
                                                                 cAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCggggg                                <  1:4636345/100‑1 (MQ=255)
                                                                    cTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCAc                             >  1:399644/1‑100 (MQ=255)
                                                                    cTTGTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCAc                             >  2:4385528/1‑100 (MQ=255)
                                                                       gTACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTgg                          <  2:3928303/100‑1 (MQ=255)
                                                                        tACTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGGACTGCAGGCTTTCACCAGCGTTCACGATGCGGTTTTGCGCGGGGGCACTGGc                         <  2:3724387/100‑1 (MQ=255)
                                                                         aCTCTCGTCATCCAGCGGTTTCACTGTGCCAACCTCAGAGGCGCTGAGCTGCAGGCTTTCACCTGCGTTCACGAGGCGTTTTTGCGCGGGGGCACTGGCg                        <  2:2629607/100‑1 (MQ=255)
                                                                         aCTCTCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCg                        >  1:3248031/1‑100 (MQ=255)
                                                                          ctctCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCga                       >  1:4684013/1‑100 (MQ=255)
                                                                          ctctCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCga                       >  1:1762897/1‑100 (MQ=255)
                                                                          ctctCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCga                       >  1:2037296/1‑100 (MQ=255)
                                                                          ctctCGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCga                       >  1:3198814/1‑100 (MQ=255)
                                                                           tctcGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCgag                      >  2:5310155/1‑100 (MQ=255)
                                                                            ctcGTCATCCAGCGGTTTCACTGTGCCAAACTCAGCGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCgaga                     <  2:2763081/100‑1 (MQ=255)
                                                                            ctcGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCgaga                     <  2:5419721/100‑1 (MQ=255)
                                                                              cGTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAg                   <  2:3441985/100‑1 (MQ=255)
                                                                              cGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAg                   <  2:3819266/100‑1 (MQ=255)
                                                                               gTCATCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGc                  >  2:3082901/1‑100 (MQ=255)
                                                                                 cagccAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGCCTTTCACCAGCGTTCACGCTGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCAc                <  2:5081775/97‑1 (MQ=255)
                                                                                  aTCCAGCGGTTTCGCTGTGCCAAACTCGGAGGCGCTGAACTGCAGGCTTTCACCAGCGTGCACGATGCGTTTTTGCGCGGGGGCACTGCCGAGAAGCACt               >  1:4912633/1‑100 (MQ=255)
                                                                                  aTCCAGCGGTTTCACTGTGCCAAACTCTGAGTCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACt               <  1:4483618/100‑1 (MQ=255)
                                                                                   tCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACtt              <  2:6498/100‑1 (MQ=255)
                                                                                   tCCAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACtt              >  2:3545854/1‑100 (MQ=255)
                                                                                    ccAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTc             >  1:2596899/1‑100 (MQ=255)
                                                                                     cAGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTcc            >  2:1961850/1‑100 (MQ=255)
                                                                                      aGCGGTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCca           >  1:3003160/1‑100 (MQ=255)
                                                                                        cGGGTTCACTGTGCCAATCTCAGAGGCGCTGAACTGCACGCTTTCACCAGCGTTCTCGGTGCTTTTTTGCGCGGGGGCACGGGCGAGAAGCACGTTcacg         <  2:103401/100‑2 (MQ=255)
                                                                                         ggTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAg        >  1:4729656/1‑100 (MQ=255)
                                                                                         ggTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAg        <  2:3864546/100‑1 (MQ=255)
                                                                                          gTTTCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGc       >  2:178872/1‑100 (MQ=255)
                                                                                            ttCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGt     <  1:1127481/100‑1 (MQ=255)
                                                                                            ttCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGt     <  2:204356/100‑1 (MQ=255)
                                                                                             tCACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtg    >  1:722533/1‑100 (MQ=255)
                                                                                              cACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtgc   <  1:1849266/100‑1 (MQ=255)
                                                                                              cACTGTGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtgc   <  2:3474274/100‑1 (MQ=255)
                                                                                                   |                                                                                               
AACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAGCGGTTTCACTGCGCCAAACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCT  >  NC_000913/4516928‑4517122

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