Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 969,158 C→T R183C (CGC→TGC)  ycaQ → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913969,1580CT87.5% 766.0 / 64.4 266R183C (CGC→TGC) ycaQhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/14);  new base T (125/107);  total (145/121)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-02

Only 100 of 273 total aligned reads displayed.

GCGGCTGGTGGGAATGGAAGCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGA  >  NC_000913/969066‑969249
                                                                                            |                                                                                           
  ggCTGGTGGGAATGGAAGCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAgtgc                                                                                           >  1:6664099/1‑93 (MQ=255)
   gCTGGTGGGAATGGAAGCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGACGGTGATTGAACGGCGCAACTTCCAGtgcg                                                                                          <  1:5312760/93‑1 (MQ=255)
     tggtggGAATGGAAGCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGtt                                                                                        <  1:4271101/93‑1 (MQ=255)
     tggtggGAATGGAAGCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGtt                                                                                        <  2:4243580/93‑1 (MQ=255)
     tggtggGAATGGAAGCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGtt                                                                                        >  1:5031404/1‑93 (MQ=255)
          ggAATGGAAGCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGa                                                                                   >  1:2024599/1‑93 (MQ=255)
          ggAATGGAAGCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGa                                                                                   >  2:5714526/1‑93 (MQ=255)
              tGGAAGCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTa                                                                               >  2:3923746/1‑93 (MQ=255)
               ggAAGCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTaa                                                                              <  2:2645277/93‑1 (MQ=255)
                gAAGCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAAc                                                                             <  1:2168004/93‑1 (MQ=255)
                  aGCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAAccc                                                                           <  2:65804/93‑1 (MQ=255)
                  aGCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAAccc                                                                           <  2:2396399/93‑1 (MQ=255)
                  aGCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAAccc                                                                           <  1:3359985/93‑1 (MQ=255)
                   gCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCa                                                                          >  1:6191610/1‑93 (MQ=255)
                   gCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCa                                                                          >  2:559281/1‑93 (MQ=255)
                   gCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCa                                                                          <  1:34111/93‑1 (MQ=255)
                    ccGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCAc                                                                         >  1:1682029/1‑93 (MQ=255)
                      gCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCg                                                                       <  2:1835713/93‑1 (MQ=255)
                       cATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCgt                                                                      >  2:4406657/1‑93 (MQ=255)
                       cATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCgt                                                                      <  2:2517624/93‑1 (MQ=255)
                       cATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCgt                                                                      <  1:765509/93‑1 (MQ=255)
                        aTAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCgtg                                                                     >  2:6809882/1‑93 (MQ=255)
                          aaaCGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTc                                                                   <  1:4274765/93‑1 (MQ=255)
                           aaCGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCa                                                                  >  2:3531759/1‑93 (MQ=255)
                             cGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATg                                                                <  1:4455030/93‑1 (MQ=255)
                               gCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGCCATGcc                                                              <  2:6988173/93‑1 (MQ=255)
                                 aTCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTg                                                            >  1:4735806/1‑93 (MQ=255)
                                  tCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGa                                                           >  1:6991131/1‑93 (MQ=255)
                                  tCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGa                                                           >  1:4903522/1‑93 (MQ=255)
                                  tCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGa                                                           >  2:5110413/1‑93 (MQ=255)
                                  tCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGa                                                           >  2:5907875/1‑93 (MQ=255)
                                  tCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGa                                                           >  2:22820/1‑93 (MQ=255)
                                   cTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGAc                                                          >  2:963752/1‑93 (MQ=255)
                                   cTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGAc                                                          <  2:7225964/93‑1 (MQ=255)
                                    tGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACt                                                         >  1:3445971/1‑93 (MQ=255)
                                    tGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACt                                                         <  1:2599215/93‑1 (MQ=255)
                                    tGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACt                                                         >  2:6074887/1‑93 (MQ=255)
                                      gAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTgg                                                       <  1:4078141/93‑1 (MQ=255)
                                      gAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTgg                                                       <  2:2857377/93‑1 (MQ=255)
                                      gAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTgg                                                       <  2:1675949/93‑1 (MQ=255)
                                         ggTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGgat                                                    <  2:2470666/93‑1 (MQ=255)
                                         ggTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGgat                                                    >  1:5649689/1‑93 (MQ=255)
                                         ggTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGgat                                                    <  1:1558873/93‑1 (MQ=255)
                                           tttatttaCTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTATCCCACCGTGTCATGCCTGACTGGgatga                                                  <  1:7200424/93‑1 (MQ=255)
                                           tttatttaCTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGgatga                                                  <  1:1037163/93‑1 (MQ=255)
                                           tttatttaCTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGgatga                                                  <  1:1379409/93‑1 (MQ=255)
                                            ttatttaCTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGgatgat                                                 <  2:4857925/93‑1 (MQ=255)
                                             tatttaCTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGatgatg                                                <  1:6726825/93‑1 (MQ=255)
                                             tatttaCTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGatgatg                                                <  1:3462056/93‑1 (MQ=255)
                                             tatttaCTGCCGGAAAGGTGATGGGGATTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGatgatg                                                >  1:3637925/1‑93 (MQ=255)
                                               tttaCTGCCGGAAAGTTGATGGTGATTGAACGGCGCAACTTCCAGGGCGTTTATGATTTAACCCACCGTGGCATGCCTGACTGGGATGATGAg                                              >  1:6816096/1‑93 (MQ=255)
                                               tttaCTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAg                                              >  2:7288868/1‑93 (MQ=255)
                                               tttaCTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAg                                              <  1:3017201/93‑1 (MQ=255)
                                                 taCTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAgcg                                            >  2:782527/1‑93 (MQ=255)
                                                    tGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGa                                         >  1:2464697/1‑93 (MQ=255)
                                                      ccGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGAtc                                       <  2:4090984/93‑1 (MQ=255)
                                                      ccGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGAtc                                       <  2:2089324/93‑1 (MQ=255)
                                                        ggAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGAtctc                                     <  1:386505/93‑1 (MQ=255)
                                                        ggAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGAtctc                                     >  2:2179220/1‑93 (MQ=255)
                                                        ggAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGAtctc                                     >  2:1992612/1‑93 (MQ=255)
                                                        ggAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGAtctc                                     <  2:1050395/93‑1 (MQ=255)
                                                         gAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCg                                    <  1:5400347/93‑1 (MQ=255)
                                                          aaaGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGt                                   <  2:4704660/93‑1 (MQ=255)
                                                          aaaGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGt                                   >  1:4622855/1‑93 (MQ=255)
                                                           aaGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGtt                                  <  1:908082/93‑1 (MQ=255)
                                                           aaGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGtt                                  <  1:5117853/93‑1 (MQ=255)
                                                           aaGGTGATGGTGATTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGtt                                  >  1:3153029/1‑93 (MQ=255)
                                                            aGGTGATGGTGATTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGttt                                 <  2:5457990/93‑1 (MQ=255)
                                                             ggtgatggtgatTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTc                                <  2:3322706/93‑1 (MQ=255)
                                                              gtgatggtgatTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCg                               >  1:927938/1‑93 (MQ=255)
                                                                  tggtgatTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCaaa                           >  1:6282767/1‑93 (MQ=255)
                                                                   ggtgatTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAAc                          <  1:1034137/93‑1 (MQ=255)
                                                                   ggtgatTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAAc                          <  1:4511447/93‑1 (MQ=255)
                                                                   ggtgatTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAAc                          >  1:4641711/1‑93 (MQ=255)
                                                                    gtgatTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACa                         <  2:5964739/93‑1 (MQ=255)
                                                                     tgatTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAg                        >  1:2357162/1‑93 (MQ=255)
                                                                     tgatTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAg                        >  1:7384105/1‑93 (MQ=255)
                                                                      gatTGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGa                       <  2:1243171/93‑1 (MQ=255)
                                                                        ttGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAg                     <  1:2194569/93‑1 (MQ=255)
                                                                        ttGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAg                     >  2:4945775/1‑93 (MQ=255)
                                                                         tGAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGc                    <  1:1126548/93‑1 (MQ=255)
                                                                          gAACGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCa                   <  1:5704646/93‑1 (MQ=255)
                                                                           aaCGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAg                  >  1:6501100/1‑93 (MQ=255)
                                                                            aCGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                 <  1:7524764/93‑1 (MQ=255)
                                                                            aCGGCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGa                 <  2:4786045/93‑1 (MQ=255)
                                                                              ggCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGaaa               <  1:2787746/93‑1 (MQ=255)
                                                                              ggCGCAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGaaa               >  1:6742451/1‑93 (MQ=255)
                                                                               gcgcAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAat              >  2:3237123/1‑93 (MQ=255)
                                                                               gcgcAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAat              >  1:5968676/1‑93 (MQ=255)
                                                                                 gcAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatca            <  2:1118582/93‑1 (MQ=255)
                                                                                  cAACTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcat           <  1:2966468/93‑1 (MQ=255)
                                                                                   aaCTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcatc          >  1:3893937/1‑93 (MQ=255)
                                                                                   aaCTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcatc          >  1:7327413/1‑93 (MQ=255)
                                                                                   aaCTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAatcatc          >  2:4116289/1‑93 (MQ=255)
                                                                                    aCTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAAtcatca         >  2:4502406/1‑93 (MQ=255)
                                                                                     cTTCCAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATcatcat        <  2:632496/93‑1 (MQ=255)
                                                                                     cTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATcatcat        <  2:6496545/93‑1 (MQ=255)
                                                                                      ttCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATg       >  2:5231056/1‑93 (MQ=255)
                                                                                        ccAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCt     <  1:302330/93‑1 (MQ=255)
                                                                                        ccAGTGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCt     >  2:1576561/1‑93 (MQ=255)
                                                                                            |                                                                                           
GCGGCTGGTGGGAATGGAAGCCGCATAAACGGCATCTGGAAGGTTTATTTACTGCCGGAAAGGTGATGGTGATTGAACGGCGCAACTTCCAGCGCGTTTATGATTTAACCCACCGTGTCATGCCTGACTGGGATGATGAGCGCGATCTCGTTTCGCAAACAGAAGCAGAAATCATCATGCTGGA  >  NC_000913/969066‑969249

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