Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,439,845 T→A 16.2% intergenic (+52/+27) ytfK → / ← ytfL DUF1107 family protein/UPF0053 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,439,8450TA16.2% 540.5 / 85.6 217intergenic (+52/+27)ytfK/ytfLDUF1107 family protein/UPF0053 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (76/104);  new base A (19/16);  total (96/121)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

Only 100 of 223 total aligned reads displayed.

TCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATC  >  NC_000913/4439746‑4439944
                                                                                                   |                                                                                                   
tCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAt                                                                                                     >  2:3721708/1‑100 (MQ=255)
tCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAt                                                                                                     <  1:4106336/100‑1 (MQ=255)
 ccGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAtt                                                                                                    >  2:2464222/1‑100 (MQ=255)
  cGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTc                                                                                                   <  2:2877479/100‑1 (MQ=255)
  cGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTc                                                                                                   >  1:3375431/1‑100 (MQ=255)
  cGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATc                                                                                                   <  2:1077764/100‑1 (MQ=255)
     aGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAgg                                                                                                <  1:3703367/100‑1 (MQ=255)
     aGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAgg                                                                                                <  1:221797/100‑1 (MQ=255)
      gTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGa                                                                                               <  1:3805716/100‑1 (MQ=255)
        tAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGc                                                                                             >  1:2514548/1‑100 (MQ=255)
         aaCCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGcc                                                                                            <  1:1121366/100‑1 (MQ=255)
          aCCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAATTGCTATGCAGTAAAAAGTGCTATGCAGTAATACTACGGCTCCTGAATCAGGAGCCg                                                                                           <  2:1474721/100‑1 (MQ=255)
          aCCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCg                                                                                           <  2:163283/100‑1 (MQ=255)
           ccGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGCCGt                                                                                          <  1:880318/100‑1 (MQ=255)
            cGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGtt                                                                                         <  1:590204/100‑1 (MQ=255)
            cGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGtt                                                                                         <  2:1372184/100‑1 (MQ=255)
            cGTCAGGTTATGCGTCTGCAAACAGAGATGGCTAAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGtt                                                                                         >  2:2526031/1‑100 (MQ=255)
              tCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGa                                                                                       <  1:3046086/100‑1 (MQ=255)
               cAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGAt                                                                                      <  1:1460799/100‑1 (MQ=255)
                aGGATATACGTCAGCACACAGACCTGGCTTAACCAAAGTACCATGCCGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGAt                                                                                      <  1:2717011/99‑1 (MQ=255)
                 ggTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGt                                                                                    >  1:2417325/1‑100 (MQ=255)
                    tATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTc                                                                                 <  1:3751169/100‑1 (MQ=255)
                    tATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTc                                                                                 <  1:962183/100‑1 (MQ=255)
                    tATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGCCGtcttattact                                                                                 <  2:3177461/100‑10 (MQ=38)
                      tGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTg                                                                               >  1:3831739/1‑100 (MQ=255)
                        cGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTggg                                                                             <  1:2662401/100‑1 (MQ=255)
                        cGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTggg                                                                             >  1:165514/1‑100 (MQ=255)
                        cGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTggg                                                                             <  2:3048590/100‑1 (MQ=255)
                         gTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTgggg                                                                            >  2:196958/1‑100 (MQ=255)
                          tCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTggggg                                                                           >  1:3385271/1‑100 (MQ=255)
                          tCTGCAAACAGAAATGGCTTAACCAACGTGCGATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTgggg                                                                            <  2:166131/99‑1 (MQ=255)
                            tGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGtt                                                                         >  1:2327332/1‑100 (MQ=255)
                            tGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGtt                                                                         >  2:1826403/1‑100 (MQ=255)
                            tGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGtt                                                                         >  2:4005905/1‑100 (MQ=255)
                              cAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGGAATACGACGGCTCATGATTCAGAAGCCGTTGATGTTTCTGGGGATTa                                                                        >  2:1228882/1‑99 (MQ=255)
                               aaaCAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAcg                                                                      <  1:3219887/100‑1 (MQ=255)
                               aaaCAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAc                                                                       <  2:479421/99‑1 (MQ=255)
                                aaCAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATCCAGGAGCCGTTGATGTTTCTGGGGGTTAcgc                                                                     <  1:2820854/100‑1 (MQ=255)
                                      gaTGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCt                                                               >  1:1633594/1‑100 (MQ=255)
                                          gCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTttct                                                           <  1:3960454/100‑1 (MQ=255)
                                          gCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTttct                                                           <  2:3367838/100‑1 (MQ=255)
                                           cTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTttctt                                                          >  1:822739/1‑100 (MQ=255)
                                             tAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGTTCCTGATTCAGGAGCCTTTGATGTTTCTGGGGGTTACGCGACGCTTtcttct                                                        <  2:4179346/100‑1 (MQ=255)
                                              aaCCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTcttctt                                                       >  2:2473372/1‑100 (MQ=255)
                                              aaCCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTcttctt                                                       <  1:198245/100‑1 (MQ=255)
                                              aaCCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTcttctt                                                       >  1:3102597/1‑100 (MQ=255)
                                               aCCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtcttt                                                      <  1:1608293/100‑1 (MQ=255)
                                                ccAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtcttta                                                     >  2:450467/1‑100 (MQ=255)
                                                ccAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtcttta                                                     <  1:1001375/100‑1 (MQ=255)
                                                 cAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGACTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttat                                                    <  1:148629/100‑1 (MQ=255)
                                                  aaagtgctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatc                                                   <  1:2396448/100‑1 (MQ=255)
                                                       gctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAg                                              >  2:1539459/1‑100 (MQ=255)
                                                         tatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCg                                            <  1:2039180/100‑1 (MQ=255)
                                                         tatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCg                                            >  2:4140420/1‑100 (MQ=255)
                                                         tatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCg                                            <  2:771034/100‑1 (MQ=255)
                                                         tatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCg                                            >  2:657364/1‑100 (MQ=255)
                                                         tatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCg                                            >  2:3810741/1‑100 (MQ=255)
                                                             cagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAg                                        <  1:2529634/100‑1 (MQ=255)
                                                              agtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAgg                                       <  1:191325/100‑1 (MQ=255)
                                                              agtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAg                                        >  1:3718357/1‑99 (MQ=255)
                                                               gtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGc                                      >  1:3907928/1‑100 (MQ=255)
                                                                 aaaaagtgctatgcagtaaTAATACGGCTCCTGAGTCAGGAGCCGTTTATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAg                                    >  1:1822775/1‑100 (MQ=255)
                                                                   aaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGTAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTGGCGTCCGGCAg                                    >  2:4039664/1‑98 (MQ=255)
                                                                    aagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGttt                                 >  1:3528065/1‑100 (MQ=255)
                                                                    aagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGttt                                 <  2:2325849/100‑1 (MQ=255)
                                                                     agtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGtttt                                <  2:2984556/100‑1 (MQ=255)
                                                                       tgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTgg                              <  1:1753192/100‑1 (MQ=255)
                                                                        gctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTg                               <  2:2768811/98‑1 (MQ=255)
                                                                        gctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGc                             <  1:3170827/100‑1 (MQ=255)
                                                                        gctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGc                             >  1:3899717/1‑100 (MQ=255)
                                                                          tatgtagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTTATTTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCg                            <  2:815627/99‑1 (MQ=255)
                                                                          tatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCg                            >  1:2158307/1‑99 (MQ=255)
                                                                          tatgcagtaaTAAGACGGCCCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGGTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTCTTGGCGa                           >  2:1023582/1‑100 (MQ=255)
                                                                            tgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTCCGCGACGCTTTATTATATATATTTAGCGTCAGGCAGTTTTGGCGaa                          >  2:1810250/1‑99 (MQ=255)
                                                                                gtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGc                     <  2:2171437/100‑1 (MQ=255)
                                                                                 taataaGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCCTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCg                    >  2:2691327/1‑100 (MQ=255)
                                                                                   ataaGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGt                  >  1:2212241/1‑100 (MQ=255)
                                                                                   ataaGACGCCTCCTGATGCAGGAGTCGTGGATGTTTCTGGGGGGTACGCGACGTTTTCTTCGTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGt                  <  2:3935169/100‑1 (MQ=255)
                                                                                    atcaaCGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGTGTTACGTGACGCTTTCTTCTTTATCTTTTGCGTCAGGCAGTTTTGGCGAAAGGGCGGTg                 >  2:3531314/5‑100 (MQ=38)
                                                                                    atcaaCGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTg                 >  2:2559109/5‑100 (MQ=255)
                                                                                     aaGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTgg                <  1:1193521/100‑1 (MQ=255)
                                                                                      aGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGc               >  1:3731561/1‑100 (MQ=255)
                                                                                       aaCGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGcc              >  1:1245675/2‑100 (MQ=255)
                                                                                       aaCGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGcc              >  2:1188855/2‑100 (MQ=255)
                                                                                        aCGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCt             >  1:777764/1‑100 (MQ=255)
                                                                                        aCGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGAATTCATCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCt             >  2:2103269/1‑100 (MQ=255)
                                                                                        aCGGCTATTGAACCAAGAGCCGTTGATGTTTCTGGGGGTTTCGCGACGCTCTCTTCTTTATTTTTAGCGTCAGGAAGTTTTGGCGCAAGGCCGCTGGCCt             >  2:1454839/1‑100 (MQ=255)
                                                                                         cGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCtt            <  1:2721701/100‑1 (MQ=255)
                                                                                         cGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCtt            >  1:1352878/1‑100 (MQ=255)
                                                                                         cGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCtt            >  2:3977990/1‑100 (MQ=255)
                                                                                          ggCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTg           >  1:670554/1‑100 (MQ=255)
                                                                                           gCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGc          <  2:1112790/100‑1 (MQ=255)
                                                                                           atctcTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTgn          <  1:2013326/96‑2 (MQ=255)
                                                                                             tCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTg        <  2:3731540/100‑1 (MQ=255)
                                                                                              ccTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGt       >  2:3253688/1‑100 (MQ=255)
                                                                                               ntGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGcagat      >  2:2013323/2‑96 (MQ=255)
                                                                                               cTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTc      <  1:1594575/100‑1 (MQ=255)
                                                                                               cTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTc      <  2:3831714/100‑1 (MQ=255)
                                                                                                tGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTTGCCTTGCTGtag     >  2:3895570/1‑98 (MQ=255)
                                                                                                   ttCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATc  <  2:742767/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGATC  >  NC_000913/4439746‑4439944

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