Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,410,599 A→G G168G (GGT→GGC envR ← putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,410,5990AG99.6% 849.8 / ‑4.7 244G168G (GGT→GGCenvRputative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (137/106);  total (138/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

Only 100 of 246 total aligned reads displayed.

TAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAGT  >  NC_000913/3410509‑3410688
                                                                                          |                                                                                         
 aaTACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCa                                                                                        >  2:5309976/1‑93 (MQ=255)
  aTACTTTATCGACCAGAGGGGGGGTTTGTTTATAAAGATCATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCAt                                                                                       <  2:3137945/93‑1 (MQ=255)
  aTACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCAt                                                                                       >  2:323567/1‑93 (MQ=255)
   tACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATc                                                                                      >  1:111683/1‑93 (MQ=255)
   tACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATc                                                                                      >  2:250251/1‑93 (MQ=255)
   tACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATc                                                                                      >  2:2450168/1‑93 (MQ=255)
   tACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATc                                                                                      >  1:4155963/1‑93 (MQ=255)
        tATCGACCAGAGGGGGGGCTGTTTTATAAAGATAATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCaataa                                                                                 <  2:2015195/93‑1 (MQ=255)
         aTCGACCAGAGGGGGGGCTTGTTTATAAAGATCATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCaataat                                                                                <  2:2358186/93‑1 (MQ=255)
         aTCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCaataaa                                                                                >  2:4101903/1‑92 (MQ=255)
          tCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCaataata                                                                               >  2:779594/1‑93 (MQ=255)
             aCCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAatc                                                                            >  1:2701982/1‑93 (MQ=255)
                agagCGGGGGTTTTTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAatcatc                                                                         <  2:920978/93‑1 (MQ=255)
                 gagCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAAtcatca                                                                        >  1:4509709/1‑93 (MQ=255)
                   gCGGGGTTTTGTTTATAAAGATCATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATcaca                                                                      <  1:5542689/93‑1 (MQ=255)
                     gggggTTTGTTAAAAAAAATCATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCacaac                                                                    <  1:2427275/93‑1 (MQ=255)
                      ggggCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCacaaca                                                                   >  1:5680629/1‑93 (MQ=255)
                      ggggCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCacaaca                                                                   >  2:1691679/1‑93 (MQ=255)
                        ggCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATc                                                                 >  1:5053918/1‑93 (MQ=255)
                         tcTTTTTTATAAACAGCAAACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCt                                                                <  1:590910/92‑1 (MQ=255)
                          tttttttATAAAGATCATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTa                                                               <  1:4345828/89‑1 (MQ=255)
                          cTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTa                                                               >  2:5251908/1‑93 (MQ=255)
                           ttGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTaa                                                              >  2:3344583/1‑93 (MQ=255)
                             ttttAAAAAGAACATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAAt                                                            <  2:4380733/92‑1 (MQ=255)
                              tttATAAAACTCATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATc                                                           <  2:1345965/93‑1 (MQ=255)
                               ttATAAAGATCATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCg                                                          <  1:289903/93‑1 (MQ=255)
                                taaaaaGATCACACGCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGa                                                         <  1:3023548/90‑1 (MQ=255)
                                 ataAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAg                                                        >  2:2190272/1‑93 (MQ=255)
                                   aaaTATAATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGt                                                      <  1:2534860/93‑1 (MQ=255)
                                   aaaGATCATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGt                                                      <  2:5432365/93‑1 (MQ=255)
                                     aGATCATCCCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTa                                                    <  1:2770130/93‑1 (MQ=255)
                                     aGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTa                                                    >  1:1040316/1‑93 (MQ=255)
                                     aGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTa                                                    >  2:1821238/1‑93 (MQ=255)
                                     aGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTa                                                    >  2:4893887/1‑93 (MQ=255)
                                      gATCATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTAt                                                   <  1:6426549/93‑1 (MQ=255)
                                      gATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCCCTGAAGGCGCCATCAATAATCATCATCACAACATCTCAATCGTGGTTAt                                                   >  1:2171895/1‑93 (MQ=255)
                                      aaTCATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTAt                                                   <  2:661818/92‑1 (MQ=255)
                                      aaTCATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTAt                                                   <  2:4312345/92‑1 (MQ=255)
                                      aaTCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTAt                                                   <  2:2995412/92‑1 (MQ=255)
                                      aaTCAAACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTAt                                                   <  2:1198479/92‑1 (MQ=255)
                                       ataaTACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTAtt                                                  <  1:5570607/90‑1 (MQ=255)
                                        tCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTAttt                                                 >  1:756149/1‑93 (MQ=255)
                                        tCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTAttt                                                 >  1:5094024/1‑93 (MQ=255)
                                         cATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTg                                                <  1:5797750/93‑1 (MQ=255)
                                         cATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTg                                                >  1:99487/1‑93 (MQ=255)
                                         aaTACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTg                                                <  2:3641522/92‑1 (MQ=255)
                                           aaccccGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCt                                              <  2:571145/90‑1 (MQ=255)
                                           aaccccGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCt                                              <  2:2301039/90‑1 (MQ=255)
                                            aaCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTa                                             >  2:4189114/1‑93 (MQ=255)
                                              cccGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTaca                                           >  2:1619810/1‑93 (MQ=255)
                                              cccGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTaca                                           >  1:6143387/1‑93 (MQ=255)
                                               ccGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTacac                                          <  2:288323/93‑1 (MQ=255)
                                                  ccATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAAc                                       >  1:5147197/1‑93 (MQ=255)
                                                   cATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAAcc                                      >  2:3915140/1‑93 (MQ=255)
                                                     tatTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCtt                                    <  1:4022680/93‑1 (MQ=255)
                                                       ttCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCttgt                                  >  2:4608420/1‑93 (MQ=255)
                                                       ttCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCttgt                                  <  2:1967976/93‑1 (MQ=255)
                                                        tCATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCttgtt                                 >  2:1954712/1‑93 (MQ=255)
                                                         cATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTATACAACCttgttg                                >  2:5748097/1‑93 (MQ=255)
                                                         cATTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCttgttg                                <  1:2382094/93‑1 (MQ=255)
                                                          aTTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGc                               >  1:3608932/1‑93 (MQ=255)
                                                          aTTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGc                               >  2:5879963/1‑93 (MQ=255)
                                                          aTTAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGc                               <  1:4987010/93‑1 (MQ=255)
                                                           ttAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCt                              >  2:5006212/1‑93 (MQ=255)
                                                           ttAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCt                              >  2:5539813/1‑93 (MQ=255)
                                                           ttAACCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCt                              <  1:1819664/93‑1 (MQ=255)
                                                             aaCCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGa                            <  1:1439711/93‑1 (MQ=255)
                                                             aaCCAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGa                            >  1:3399419/1‑93 (MQ=255)
                                                               ccAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGaca                          >  1:839331/1‑93 (MQ=255)
                                                                cAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGacac                         >  1:3217260/1‑93 (MQ=255)
                                                                cAGTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGacac                         <  1:4826319/93‑1 (MQ=255)
                                                                  gTTTTGATCAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGc                       <  2:4682358/93‑1 (MQ=255)
                                                                  gTTTTGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGc                       >  2:4471936/1‑93 (MQ=255)
                                                                   ttttGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGcc                      <  1:6031235/93‑1 (MQ=255)
                                                                   ttttGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGcc                      >  1:4660503/1‑93 (MQ=255)
                                                                   ttttGAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGcc                      <  1:712261/93‑1 (MQ=255)
                                                                       gAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCa                  <  2:2726486/93‑1 (MQ=255)
                                                                       gAACAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCa                  <  1:3533397/93‑1 (MQ=255)
                                                                        aaCAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCaa                 >  1:2962942/1‑93 (MQ=255)
                                                                        aaCAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCaa                 <  2:621158/93‑1 (MQ=255)
                                                                          cAATTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATa               <  2:3563137/93‑1 (MQ=255)
                                                                           aaTTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATAc              >  2:3051759/1‑93 (MQ=255)
                                                                           aaTTCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTGGTTGCTGACACGCCTGCAATAc              <  2:3766985/93‑1 (MQ=255)
                                                                             ttCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACtt            <  2:5464723/93‑1 (MQ=255)
                                                                             tgccGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACtt            >  1:6370299/3‑93 (MQ=255)
                                                                              tCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTc           >  2:5065440/1‑93 (MQ=255)
                                                                              tCCGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTc           <  2:5902000/93‑1 (MQ=255)
                                                                               ccGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTcg          >  1:5403714/1‑93 (MQ=255)
                                                                               ccGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTcg          <  2:4430645/93‑1 (MQ=255)
                                                                                cGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTcgc         >  1:1921737/1‑93 (MQ=255)
                                                                                cGCTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTcgc         <  1:2117440/93‑1 (MQ=255)
                                                                                  cTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCgg       <  2:2633659/93‑1 (MQ=255)
                                                                                  cTGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCgg       >  1:5211220/1‑93 (MQ=255)
                                                                                   tGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGga      <  2:3396528/93‑1 (MQ=255)
                                                                                   tGAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGga      <  2:2618449/93‑1 (MQ=255)
                                                                                    gAAGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGgag     <  1:4312244/93‑1 (MQ=255)
                                                                                     aaGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGgaga    >  1:2855119/1‑93 (MQ=255)
                                                                                      aGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGagag   >  2:2816037/1‑93 (MQ=255)
                                                                                      aGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGagag   >  2:6024111/1‑93 (MQ=255)
                                                                                      aGGCGCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGagag   >  2:854196/1‑93 (MQ=255)
                                                                                          |                                                                                         
TAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAGT  >  NC_000913/3410509‑3410688

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