Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,139,669 +C coding (233/330 nt)
coding (132/159 nt)
yijI →
b4558 ←
hypothetical protein
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,139,6671.C99.7% 998.1 / ‑4.7 313coding (231/330 nt)
coding (134/159 nt)
yijI
b4558
hypothetical protein
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (136/176);  total (136/177)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 343 total aligned reads displayed.

GCTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCAT  >  NC_000913/4139575‑4139754
                                                                                             |                                                                                       
gcTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTg                                                                                          >  2:5395584/1‑93 (MQ=255)
gcTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTg                                                                                          >  2:500300/1‑93 (MQ=255)
   tCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTTCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccc                                                                                      <  1:6228674/93‑2 (MQ=255)
   tCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccc                                                                                      >  2:1388891/1‑92 (MQ=255)
   tCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccc                                                                                      >  2:6183062/1‑92 (MQ=255)
     cTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccctt                                                                                    >  1:4564496/1‑90 (MQ=255)
     cTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccctt                                                                                    >  2:198688/1‑90 (MQ=255)
     cTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGccctt                                                                                    <  2:2642989/93‑4 (MQ=255)
            tcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCa                                                                              <  2:673291/93‑1 (MQ=255)
            tcgGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCa                                                                              <  1:5605880/93‑1 (MQ=255)
               gCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTg                                                                           >  2:5817927/1‑93 (MQ=255)
                ccccGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTgc                                                                          <  1:5242257/93‑1 (MQ=255)
                 cccGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTgcg                                                                         >  2:7032037/1‑93 (MQ=255)
                   cGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCa                                                                       <  1:3377394/93‑1 (MQ=255)
                     gAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAgg                                                                     >  1:6698626/1‑93 (MQ=255)
                      aaaTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGt                                                                    <  1:3220856/93‑1 (MQ=255)
                      aaaTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGt                                                                    <  1:343798/93‑1 (MQ=255)
                      aaaTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGt                                                                    >  2:7893490/1‑93 (MQ=255)
                        aTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTa                                                                  <  1:5472537/93‑1 (MQ=255)
                        aTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTa                                                                  <  2:6224672/93‑1 (MQ=255)
                        aTCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTa                                                                  <  1:6866929/93‑1 (MQ=255)
                        aTCAGTAACTGGCAGTTATCGCCCGCCAGCCTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTa                                                                  <  1:5302459/93‑1 (MQ=255)
                         tCAGTAACTGACAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTAt                                                                 <  1:3931051/93‑1 (MQ=255)
                           aGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGc                                                               <  1:508622/93‑1 (MQ=255)
                           aGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGc                                                               <  1:1070419/93‑1 (MQ=255)
                            gTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGcc                                                              >  1:2462675/1‑93 (MQ=255)
                             tAACTGGCAGTTATCGCCCGCCAGCGCATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGcaa                                                             >  2:1448514/1‑91 (MQ=255)
                              aaCTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCAc                                                            >  2:2343228/1‑93 (MQ=255)
                                cTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTc                                                          <  1:4354988/93‑1 (MQ=255)
                                cTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTc                                                          >  1:2180364/1‑93 (MQ=255)
                                 tGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCTGTGCGCAGGTTATGCCACTCa                                                         >  1:2866335/1‑93 (MQ=255)
                                 tGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCa                                                         <  2:1093214/93‑1 (MQ=255)
                                  ggCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCaa                                                        <  2:4558925/93‑1 (MQ=255)
                                   gCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAAt                                                       <  2:7651210/93‑1 (MQ=255)
                                    cAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAAtt                                                      <  2:4385046/93‑1 (MQ=255)
                                    cAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAAtt                                                      <  1:7244639/93‑1 (MQ=255)
                                     aGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAAttt                                                     <  2:4917047/93‑1 (MQ=255)
                                       ttATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtgt                                                   >  2:5656298/1‑93 (MQ=255)
                                       ttATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtgt                                                   >  2:3182480/1‑93 (MQ=255)
                                       ttATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTtgt                                                   >  1:425733/1‑93 (MQ=255)
                                         aTCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtga                                                 >  1:3651038/1‑93 (MQ=255)
                                          tCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtgat                                                >  1:1226970/1‑93 (MQ=255)
                                          tCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtgat                                                <  2:3814227/93‑1 (MQ=255)
                                          tCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGtgat                                                <  1:1606017/93‑1 (MQ=255)
                                              ccgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaa                                            >  2:6865076/1‑93 (MQ=255)
                                              ccgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaa                                            >  2:7306978/1‑93 (MQ=255)
                                               cgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaaa                                           <  1:5751276/93‑1 (MQ=255)
                                               cgccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaaa                                           >  1:5186263/1‑93 (MQ=255)
                                                gccAGTGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaaaa                                          >  2:5200737/1‑93 (MQ=255)
                                                gccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGaaaaa                                          <  2:4690909/93‑1 (MQ=255)
                                                 ccAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAg                                         <  2:4092461/93‑1 (MQ=255)
                                                  cAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGt                                        <  1:7478496/93‑1 (MQ=255)
                                                     cGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTAtt                                     <  1:2028712/93‑1 (MQ=255)
                                                       tATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTac                                   <  2:7782653/93‑1 (MQ=255)
                                                       tATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTac                                   <  2:4263973/93‑1 (MQ=255)
                                                        aTCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTaca                                  >  1:928756/1‑93 (MQ=255)
                                                        aTCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTaca                                  >  1:595217/1‑93 (MQ=255)
                                                         tCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTacac                                 >  1:6913018/1‑93 (MQ=255)
                                                         tCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTacac                                 <  2:334116/93‑1 (MQ=255)
                                                          cGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTAcaca                                <  2:1729933/93‑1 (MQ=255)
                                                           ggTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAg                               <  2:3660408/93‑1 (MQ=255)
                                                           ggTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAg                               >  1:6959497/1‑93 (MQ=255)
                                                           ggTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAg                               >  2:7489217/1‑93 (MQ=255)
                                                            gTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGc                              <  2:252089/93‑1 (MQ=255)
                                                             tGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCg                             >  1:4152572/1‑93 (MQ=255)
                                                              gCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGt                            >  1:1028843/1‑93 (MQ=255)
                                                               ccAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGtt                           <  2:3714547/93‑1 (MQ=255)
                                                                cAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGttt                          <  1:2432204/93‑1 (MQ=255)
                                                                cAATCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGttt                          <  1:5038739/93‑1 (MQ=255)
                                                                 aaTCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTc                         >  1:2833238/1‑93 (MQ=255)
                                                                   tCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCaa                       <  1:4889219/93‑1 (MQ=255)
                                                                   tCAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCaa                       >  1:3747101/1‑93 (MQ=255)
                                                                   acAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCaa                       <  2:870258/92‑1 (MQ=255)
                                                                    cAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAAc                      <  2:6440746/93‑1 (MQ=255)
                                                                    cAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAAc                      <  1:4467932/93‑1 (MQ=255)
                                                                    cAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAAc                      <  2:4750996/93‑1 (MQ=255)
                                                                    cAGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAAc                      <  2:2788320/93‑1 (MQ=255)
                                                                     aGCGCAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACg                     <  2:167670/93‑1 (MQ=255)
                                                                      gcgcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGt                    >  1:5554441/1‑93 (MQ=255)
                                                                       cgcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTg                   <  1:2664673/93‑1 (MQ=255)
                                                                       cgcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTg                   <  2:1320762/93‑1 (MQ=255)
                                                                        gcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTgg                  >  2:7154988/1‑93 (MQ=255)
                                                                        gcAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTgg                  >  1:6472041/1‑93 (MQ=255)
                                                                         cAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGc                 >  2:5101655/1‑93 (MQ=255)
                                                                         cAAGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGc                 <  1:497567/93‑1 (MQ=255)
                                                                           aGGGCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCtt               <  2:7604441/93‑1 (MQ=255)
                                                                             ggCGCTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGc             >  1:1469660/1‑93 (MQ=255)
                                                                              gcgcTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGcc            <  1:1823607/93‑1 (MQ=255)
                                                                              gcgcTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGcc            <  2:5890520/93‑1 (MQ=255)
                                                                               cgcTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCg           >  2:6488046/1‑93 (MQ=255)
                                                                               cgcTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCg           <  2:4162407/93‑1 (MQ=255)
                                                                               cgcTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCg           <  1:4616614/93‑1 (MQ=255)
                                                                                gcTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCgg          <  2:5610238/93‑1 (MQ=255)
                                                                                gcTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCgg          >  1:4111200/1‑93 (MQ=255)
                                                                                 cTTTTGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGa         <  1:7773881/93‑1 (MQ=255)
                                                                                  ttttGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGAc        >  2:629300/1‑93 (MQ=255)
                                                                                  ttttGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGAc        >  2:2118428/1‑93 (MQ=255)
                                                                                    ttGGCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGAcgc      <  1:1439244/93‑1 (MQ=255)
                                                                                      ggCGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACgcgc    >  2:3471481/1‑93 (MQ=255)
                                                                                        cGTTGCCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCAt  <  1:4725146/93‑1 (MQ=255)
                                                                                             |                                                                                       
GCTTCCTGTTCGTCGGCCCCGGAAATCAGTAACTGGCAGTTATCGCCCGCCAGCGTATCGGTGCCAATCAGCGCAAGGGCGCTTTTGGCGTTGCCTTGCGGTCAGTGCGCAGGTTATGCCACTCAATTTGTGATGAAAAAGTATTACACAGCGTTTCAACGTGGCTTGCCGGACGCGCAT  >  NC_000913/4139575‑4139754

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