Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 269,425 A→G intergenic (‑19/‑41) perR ← / → yi91a putative transcriptional regulator LYSR‑type/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913269,4250AG84.4% 691.9 / 123.2 301intergenic (‑19/‑41)perR/yi91aputative transcriptional regulator LYSR‑type/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/22);  new base G (121/133);  total (146/155)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 314 total aligned reads displayed.

AAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATT  >  NC_000913/269335‑269517
                                                                                          |                                                                                            
aaaaGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTggg                                                                                            <  1:108176/93‑1 (MQ=255)
aaaaGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTggg                                                                                            <  1:6198235/93‑1 (MQ=255)
aaaaGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTggg                                                                                            >  2:4420680/1‑93 (MQ=255)
 aaaGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTgggg                                                                                           >  1:6630291/1‑93 (MQ=255)
 aaaGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTgggg                                                                                           >  2:3943614/1‑93 (MQ=255)
 aaaGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTgggg                                                                                           >  1:2479431/1‑93 (MQ=255)
 aaaGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTgggg                                                                                           >  2:5896267/1‑93 (MQ=255)
  aaGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGa                                                                                          >  1:6944819/1‑93 (MQ=255)
  aaGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGa                                                                                          <  1:474394/93‑1 (MQ=255)
     cTCCGGTACGATCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGa                                                                                       >  2:4654229/1‑93 (MQ=255)
     cTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGa                                                                                       <  1:6105522/93‑1 (MQ=255)
        cGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGt                                                                                    <  1:2377837/93‑1 (MQ=255)
        cGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGt                                                                                    <  2:2626352/93‑1 (MQ=255)
         ggTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTg                                                                                   <  1:5810921/93‑1 (MQ=255)
         ggTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTg                                                                                   >  1:509819/1‑93 (MQ=255)
         ggTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTg                                                                                   <  2:3012588/93‑1 (MQ=255)
         ggTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTg                                                                                   <  2:3814367/93‑1 (MQ=255)
          gTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTgg                                                                                  <  1:5702956/93‑1 (MQ=255)
               aCCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTAc                                                                             <  2:5401750/93‑1 (MQ=255)
               aCCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTAc                                                                             <  2:2430137/93‑1 (MQ=255)
                    ccGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAAt                                                                        <  2:4353286/93‑1 (MQ=255)
                    ccGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAAt                                                                        <  2:2249753/93‑1 (MQ=255)
                     cGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAAtt                                                                       <  2:7458279/93‑1 (MQ=255)
                     cGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAAtt                                                                       <  2:2473914/93‑1 (MQ=255)
                      gCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGACTCCTTTTCACTAGGGATGAAGTGGTTTACTGAAttt                                                                      >  1:2894406/1‑93 (MQ=255)
                        ttCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTgg                                                                    >  1:1162943/1‑93 (MQ=255)
                        ttCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTgg                                                                    <  2:1062869/93‑1 (MQ=255)
                         tCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGc                                                                   >  1:1003477/1‑93 (MQ=255)
                          cAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGcc                                                                  <  1:7312323/93‑1 (MQ=255)
                           aaaGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCa                                                                 >  2:2127144/1‑93 (MQ=255)
                              ggCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCt                                                              <  2:5176456/93‑1 (MQ=255)
                              ggCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCt                                                              <  1:319840/93‑1 (MQ=255)
                                 gcgTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGaa                                                           >  2:4188257/1‑93 (MQ=255)
                                   gTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACa                                                         >  1:2313105/1‑93 (MQ=255)
                                     aaCAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACaga                                                       <  1:3947378/93‑1 (MQ=255)
                                     aaCAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACaga                                                       >  1:4813065/1‑93 (MQ=255)
                                      aCAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACagag                                                      >  2:3136086/1‑93 (MQ=255)
                                       cAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAgg                                                     <  2:4821848/93‑1 (MQ=255)
                                       cAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAgg                                                     <  2:5469021/93‑1 (MQ=255)
                                        aGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGt                                                    <  2:7329663/93‑1 (MQ=255)
                                         gATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTg                                                   >  1:293906/1‑93 (MQ=255)
                                         gATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTg                                                   >  2:6563730/1‑93 (MQ=255)
                                           tttAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGat                                                 >  2:1215598/1‑93 (MQ=255)
                                            ttAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGata                                                <  2:5691006/93‑1 (MQ=255)
                                             tAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGatat                                               <  2:85424/93‑1 (MQ=255)
                                              agagGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATg                                              >  1:5964929/1‑93 (MQ=255)
                                                agGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCt                                            <  1:6031916/93‑1 (MQ=255)
                                                 ggTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTc                                           >  2:714988/1‑93 (MQ=255)
                                                  gTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCa                                          <  2:3341916/93‑1 (MQ=255)
                                                   tGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCAc                                         >  1:1290545/1‑93 (MQ=255)
                                                     cTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCt                                       >  2:5247244/1‑93 (MQ=255)
                                                     cTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCt                                       <  2:6909987/93‑1 (MQ=255)
                                                     cTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCt                                       >  2:6463355/1‑93 (MQ=255)
                                                      ttttGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTc                                      <  1:331213/93‑1 (MQ=255)
                                                       tttGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCa                                     <  2:4514281/93‑1 (MQ=255)
                                                       tttGCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCa                                     <  1:1617339/93‑1 (MQ=255)
                                                         tGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGa                                   <  1:3403792/93‑1 (MQ=255)
                                                          gCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGaa                                  <  2:5102277/93‑1 (MQ=255)
                                                          gCTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGaa                                  >  1:6747806/1‑93 (MQ=255)
                                                           cTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGaac                                 <  2:3519335/93‑1 (MQ=255)
                                                           cTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGaac                                 >  1:2074274/1‑93 (MQ=255)
                                                           cTAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGaac                                 <  1:4007287/93‑1 (MQ=255)
                                                            tAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGaaca                                >  2:6377509/1‑93 (MQ=255)
                                                            tAAGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGaaca                                >  2:6023172/1‑93 (MQ=255)
                                                             aaGAGCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGaacaa                               <  1:7160415/93‑1 (MQ=255)
                                                              agagCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGaacaac                              <  2:1161541/93‑1 (MQ=255)
                                                               gagCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAacaaca                             <  2:6306311/93‑1 (MQ=255)
                                                               gagCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAacaaca                             <  1:24383/93‑1 (MQ=255)
                                                               gaacTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAacaaca                             <  1:4235989/90‑1 (MQ=255)
                                                                agCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAacac                            <  1:4447091/93‑1 (MQ=255)
                                                                agCTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAacac                            <  1:1505419/93‑1 (MQ=255)
                                                                  cTTCATTCATGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAg                          >  2:6761073/1‑93 (MQ=255)
                                                                   ttcattcaTGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAgg                         <  2:1627998/93‑1 (MQ=255)
                                                                   ttcattcaTGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAgg                         >  2:1260737/1‑93 (MQ=255)
                                                                    tcattcatGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACGGAGGTGATATGCTCACCTCAGAACAACACAGgt                        >  2:3698370/1‑93 (MQ=255)
                                                                    tcattcatGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGgt                        <  1:4557230/93‑1 (MQ=255)
                                                                    tcattcatGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGgt                        <  2:3578807/93‑1 (MQ=255)
                                                                     cattcatGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGgtg                       >  1:1338886/1‑93 (MQ=255)
                                                                      attcatGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGgtgt                      <  1:3505333/93‑1 (MQ=255)
                                                                      attcatGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGgtgt                      >  1:3123559/1‑93 (MQ=255)
                                                                      attcatGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGgtgt                      >  1:2797255/1‑93 (MQ=255)
                                                                      attcatGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGgtgt                      >  1:3837323/1‑93 (MQ=255)
                                                                      attcatGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGgtgt                      >  1:4288432/1‑93 (MQ=255)
                                                                       ttcatGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTc                     <  1:5412307/93‑1 (MQ=255)
                                                                        tcatGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCa                    >  2:4980334/1‑93 (MQ=255)
                                                                        tcatGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCa                    <  2:4618509/93‑1 (MQ=255)
                                                                        tcatGAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCa                    >  2:7028388/1‑93 (MQ=255)
                                                                            gAATCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAAt                >  1:636411/1‑93 (MQ=255)
                                                                            gAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAAt                <  1:451168/93‑1 (MQ=255)
                                                                             aaTCCTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATg               >  2:4267371/1‑93 (MQ=255)
                                                                              aTCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGa              <  2:1426075/93‑1 (MQ=255)
                                                                               tCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGaa             >  1:1380098/1‑93 (MQ=255)
                                                                                ccTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGaaa            <  1:6842871/93‑1 (MQ=255)
                                                                                 cTTTTCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGaaaa           >  1:83172/1‑93 (MQ=255)
                                                                                  ttttCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGaaaaa          >  1:2635947/1‑93 (MQ=255)
                                                                                   tttCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGaaaaaa         >  2:3915681/1‑93 (MQ=255)
                                                                                    ttCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGaaaaaaa        <  1:4331734/93‑1 (MQ=255)
                                                                                     tCACTGGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAg       <  1:2159599/93‑1 (MQ=255)
                                                                                         tAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAAt   <  2:2790706/93‑1 (MQ=255)
                                                                                          aGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAAtt  <  1:5292457/93‑1 (MQ=255)
                                                                                          |                                                                                            
AAAAGCTCCGGTACGACCAGCCGCTTCAAAGGCGCGTAACAGATTTAGAGGTGCTTTTGCTAAGAGCTTCATTCATGAATCCTTTTCACTAGGGATGAAGTGGTTTACTGAATTTGGCCACCTGAACAGAGGTGATATGCTCACCTCAGAACAACACAGGTGTCATAATGAAAAAAAGAAATT  >  NC_000913/269335‑269517

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