Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 823,777 (A)5→6 intergenic (‑57/‑76) ybhP ← / → ybhQ hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913823,7721.A98.8% 575.9 / ‑0.8 164intergenic (‑52/‑81)ybhP/ybhQhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (85/77);  total (86/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.78e-01

Only 100 of 187 total aligned reads displayed.

GTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCAC‑‑AAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGT  >  NC_000913/823680‑823859
                                                                                             |                                                                                         
gTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTcac                                                                                            <  2:1823842/93‑1 (MQ=255)
gTGAGCACCTTGAACGAAAATTGTTGTGTATGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTcac                                                                                            <  1:2456709/93‑1 (MQ=255)
 tGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTcac‑‑a                                                                                         <  1:590118/93‑1 (MQ=255)
      aCCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCAC‑‑aaaaaa                                                                                   <  2:279478/93‑2 (MQ=255)
      aCCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCAC‑‑aaaaaa                                                                                   <  2:2702902/93‑2 (MQ=255)
      aCCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCAC‑‑aaaaaa                                                                                   >  1:2363308/1‑92 (MQ=255)
      aCCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCAC‑‑aaaaaa                                                                                   <  2:228158/93‑2 (MQ=255)
      aCCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCAC‑‑aaaaaa                                                                                   >  2:4727440/1‑92 (MQ=255)
       ccTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCAC‑‑aaaaaag                                                                                  >  2:1758105/1‑91 (MQ=255)
        cTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCAC‑‑aaaaaagg                                                                                 <  2:1520272/93‑4 (MQ=255)
         ttGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCAC‑‑aaaaaaggt                                                                                <  2:4479059/93‑5 (MQ=255)
           gAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGc                                                                               <  1:1447788/93‑1 (MQ=255)
           gAACGAAAATTGTTGTGTGTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGc                                                                               <  1:4618542/93‑1 (MQ=255)
            aaCGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCg                                                                              <  1:729947/93‑1 (MQ=255)
             aCGAAAATTGTTGTGTTTTATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCGCATTGAGATAGTGTAGTCTGCGTCACAAAAAAGGTGCgg                                                                             >  1:2939025/1‑93 (MQ=255)
             aCGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCgg                                                                             >  2:3265328/1‑93 (MQ=255)
                aaaaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCt                                                                          <  2:3718211/93‑1 (MQ=255)
                 aaaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCAC‑‑AAAAAGGTGCGGTCTTa                                                                        <  2:4410539/93‑1 (MQ=255)
                  aaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTa                                                                        >  1:1662070/1‑93 (MQ=255)
                  aaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTa                                                                        <  2:3967077/93‑1 (MQ=255)
                  aaTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGCGTAGTCAGCGTCACAAAAAAGGTGCGGTCTTa                                                                        >  2:565892/1‑93 (MQ=255)
                   aTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTAc                                                                       >  1:2732493/1‑93 (MQ=255)
                   aTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTAc                                                                       >  1:3933521/1‑93 (MQ=255)
                   aTTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTAc                                                                       >  1:134152/1‑93 (MQ=255)
                    ttgttgTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACg                                                                      <  1:1282836/93‑1 (MQ=255)
                     tgttgtGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACgg                                                                     >  2:3213570/1‑93 (MQ=255)
                       ttgtGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGaa                                                                   >  1:3365935/1‑93 (MQ=255)
                        tgtgTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAAt                                                                  <  1:1480719/93‑1 (MQ=255)
                          tgtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAAttt                                                                <  2:427835/93‑1 (MQ=255)
                           gtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAAGGTGCGGTCTTACGGAAttt                                                                <  1:1884048/93‑1 (MQ=255)
                           gtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAAtttt                                                               >  2:2400634/1‑93 (MQ=255)
                           gtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAAtttt                                                               >  1:2583546/1‑93 (MQ=255)
                           gtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCAC‑‑AAAACAGGTGCGGGCTTCCGGAAttgt                                                               >  1:2583547/1‑91 (MQ=255)
                               gATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCTGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGt                                                           >  2:516904/1‑93 (MQ=255)
                                aTCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTa                                                          >  2:812384/1‑93 (MQ=255)
                                 tCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTaa                                                         <  1:3176000/93‑1 (MQ=255)
                                  cGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTaaa                                                        >  2:4735590/1‑93 (MQ=255)
                                  cGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTaaa                                                        <  2:2496456/93‑1 (MQ=255)
                                   gggCATAATTTTCCTGTCTTTGCATCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAg                                                       <  2:2593089/93‑1 (MQ=255)
                                    ggCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGt                                                      >  1:1211185/1‑93 (MQ=255)
                                       aTAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCg                                                   >  1:1056863/1‑93 (MQ=255)
                                       aTAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCg                                                   >  2:226299/1‑93 (MQ=255)
                                       aTAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCg                                                   >  1:2115254/1‑93 (MQ=255)
                                       aTAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCg                                                   >  1:31259/1‑93 (MQ=255)
                                           ttttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTAc                                               >  1:3696407/1‑93 (MQ=255)
                                           ttttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTAc                                               <  1:320886/93‑1 (MQ=255)
                                           ttttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTAc                                               >  2:3788926/1‑93 (MQ=255)
                                           ttttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTAc                                               <  1:43819/93‑1 (MQ=255)
                                             ttCGTGTTTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTActc                                             <  1:4175373/93‑1 (MQ=255)
                                              tCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTActct                                            <  1:3388126/93‑1 (MQ=255)
                                              tCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTActct                                            >  1:4297358/1‑93 (MQ=255)
                                                cTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTgg                                          >  1:3094689/1‑92 (MQ=255)
                                                cTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGa                                          <  2:2522789/93‑1 (MQ=255)
                                                cTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGa                                          >  1:988919/1‑93 (MQ=255)
                                                 tGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGag                                         <  2:175272/93‑1 (MQ=255)
                                                 tGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGag                                         <  1:1118779/93‑1 (MQ=255)
                                                 tGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGag                                         >  2:3156139/1‑93 (MQ=255)
                                                  gTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagt                                        <  1:642533/93‑1 (MQ=255)
                                                  gTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagt                                        <  1:638514/93‑1 (MQ=255)
                                                    cTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaa                                      <  1:3630510/93‑1 (MQ=255)
                                                     tttGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaag                                     >  1:3518409/1‑93 (MQ=255)
                                                      ttGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaagt                                    <  2:4370292/93‑1 (MQ=255)
                                                      ttGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCAC‑‑AAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaagta                                   >  1:4105250/1‑93 (MQ=255)
                                                       tGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaagta                                   >  1:4716986/1‑93 (MQ=255)
                                                       tGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaagta                                   <  2:702439/93‑1 (MQ=255)
                                                        gCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTag                                  >  1:3527225/1‑93 (MQ=255)
                                                        gCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTag                                  >  1:1667309/1‑93 (MQ=255)
                                                         ccTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTaga                                 >  2:1449995/1‑93 (MQ=255)
                                                           tctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAgaga                               >  2:661141/1‑93 (MQ=255)
                                                           tctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAgaga                               >  2:1636271/1‑93 (MQ=255)
                                                           tctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAgaga                               <  1:3359146/93‑1 (MQ=255)
                                                             ttATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATa                             >  1:2784604/1‑93 (MQ=255)
                                                                tctcATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGATTAAGTAGAGATAAAt                          >  1:4208244/1‑93 (MQ=255)
                                                                tctcATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAt                          <  2:362045/93‑1 (MQ=255)
                                                                  tcATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttc                        >  2:3506040/1‑93 (MQ=255)
                                                                  tcATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttc                        >  2:865469/1‑93 (MQ=255)
                                                                  tcATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttc                        >  1:91192/1‑93 (MQ=255)
                                                                  tcATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttc                        <  2:1395878/93‑1 (MQ=255)
                                                                   cATTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttct                       >  2:1285863/1‑93 (MQ=255)
                                                                    aTTGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttctt                      <  1:2227385/93‑1 (MQ=255)
                                                                     ttGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttcttc                     <  2:3104991/93‑1 (MQ=255)
                                                                     ttGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttcttc                     >  1:4558463/1‑93 (MQ=255)
                                                                     ttGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttcttc                     >  2:1149699/1‑93 (MQ=255)
                                                                     ttGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttcttc                     <  1:2636254/93‑1 (MQ=255)
                                                                     ttGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttcttc                     >  1:2308212/1‑93 (MQ=255)
                                                                     ttGAAATAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttcttc                     <  2:1591439/93‑1 (MQ=255)
                                                                         aaTAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGga                 <  2:1267047/93‑1 (MQ=255)
                                                                         aaTAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGga                 <  2:1864272/93‑1 (MQ=255)
                                                                          aTAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGCGTAAGTAGATATAAATTCTTCAGgag                >  1:3781414/1‑93 (MQ=255)
                                                                          aTAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGgag                >  1:705740/1‑93 (MQ=255)
                                                                          aTAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGgag                >  2:454493/1‑93 (MQ=255)
                                                                           tAGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGgaga               <  1:3857112/93‑1 (MQ=255)
                                                                            aGTGTAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGgagag              <  2:3984293/93‑1 (MQ=255)
                                                                             gtgtAGTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGgagaga             >  2:2965709/1‑93 (MQ=255)
                                                                                  gTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCa        >  1:4064056/1‑93 (MQ=255)
                                                                                  gTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCa        >  1:522524/1‑93 (MQ=255)
                                                                                  gTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCa        >  2:4393189/1‑93 (MQ=255)
                                                                                  gTCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCa        <  2:3683843/93‑1 (MQ=255)
                                                                                   tCGGCGTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCAt       >  2:140302/1‑93 (MQ=255)
                                                                                        gTCACAAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGt  >  2:220746/1‑93 (MQ=255)
                                                                                             |                                                                                         
GTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCAC‑‑AAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGT  >  NC_000913/823680‑823859

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