Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,199,211 T→C V273V (GTA→GTG intE ← prophage e14 integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,199,2110TC99.6% 913.6 / ‑5.0 270V273V (GTA→GTGintEprophage e14 integrase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (132/136);  total (133/137)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

Only 100 of 287 total aligned reads displayed.

GGCTTTAACCTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCTACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCTTTTCCTGAATGACATG  >  NC_000913/1199121‑1199300
                                                                                          |                                                                                         
ggCTTTAACCTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGccac                                                                                         >  2:6696856/1‑93 (MQ=255)
ggCTTTAACCTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGccac                                                                                         >  2:4040314/1‑93 (MQ=255)
ggCTTTAACCTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGccac                                                                                         <  2:3817519/93‑1 (MQ=255)
 gCTTTAACCTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGccact                                                                                        <  1:7150877/93‑1 (MQ=255)
 gCTTTAACCTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGccact                                                                                        <  2:4718636/93‑1 (MQ=255)
  cTTTAACCTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCATCGGGccacta                                                                                       <  1:792435/93‑1 (MQ=255)
  cTTTAACCTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGccacta                                                                                       <  1:174751/93‑1 (MQ=255)
   tttAACCTGCGCTGCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATTCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCcactac                                                                                      <  1:1959423/93‑1 (MQ=255)
   tttAACCTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCcactac                                                                                      >  2:6047844/1‑93 (MQ=255)
     tAACCTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACAt                                                                                    <  2:7587591/93‑1 (MQ=255)
     tAACCTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACAt                                                                                    <  2:1251336/93‑1 (MQ=255)
      aaCCTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATc                                                                                   <  1:3110823/93‑1 (MQ=255)
       aCCTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATcg                                                                                  >  2:2399695/1‑93 (MQ=255)
       aCCTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATcg                                                                                  <  1:6080677/93‑1 (MQ=255)
        ccTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATcgc                                                                                 >  1:4006062/1‑93 (MQ=255)
          tGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATcgcga                                                                               >  1:6586427/1‑92 (MQ=255)
          tGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGt                                                                               >  1:3199361/1‑93 (MQ=255)
          tGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGt                                                                               >  1:3387123/1‑93 (MQ=255)
          tGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGt                                                                               >  1:810872/1‑93 (MQ=255)
          tGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGt                                                                               >  1:1538937/1‑93 (MQ=255)
            cgcTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTaa                                                                             <  1:3804622/93‑1 (MQ=255)
             gcTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTTCATCGCGTaaa                                                                            <  2:6808268/93‑1 (MQ=255)
              cTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAAc                                                                           >  2:7127107/1‑93 (MQ=255)
               tCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACt                                                                          <  2:877705/93‑1 (MQ=255)
               tCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACt                                                                          <  1:7273261/93‑1 (MQ=255)
               tCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACt                                                                          >  1:5273715/1‑93 (MQ=255)
                  gcgcTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCa                                                                       <  2:3404388/93‑1 (MQ=255)
                  gcgcTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCa                                                                       <  2:5054945/93‑1 (MQ=255)
                  gcgcTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCa                                                                       <  1:4754116/93‑1 (MQ=255)
                  gcgcTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAAATCCa                                                                       >  1:11219/1‑93 (MQ=255)
                    gcTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGt                                                                     >  1:654860/1‑93 (MQ=255)
                     ctctGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGtt                                                                    >  1:5534468/1‑93 (MQ=255)
                       ctGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTaa                                                                  <  1:7577674/93‑1 (MQ=255)
                       ctGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTaa                                                                  >  2:5422588/1‑93 (MQ=255)
                        tGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAAt                                                                 >  2:51852/1‑93 (MQ=255)
                        tGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAAt                                                                 >  1:6213341/1‑93 (MQ=255)
                          cTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCg                                                               <  1:3014046/93‑1 (MQ=255)
                           ttGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGc                                                              <  1:6307274/93‑1 (MQ=255)
                            tgtgAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCa                                                             >  1:1352147/1‑93 (MQ=255)
                              tgAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCAtt                                                           <  1:7071167/93‑1 (MQ=255)
                              tgAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCAtt                                                           <  1:3505217/93‑1 (MQ=255)
                                aaGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGa                                                         <  2:2542941/93‑1 (MQ=255)
                                 aGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAg                                                        <  1:6266313/93‑1 (MQ=255)
                                 aGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAg                                                        >  2:4177730/1‑93 (MQ=255)
                                  gTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAgg                                                       <  2:595927/93‑1 (MQ=255)
                                  gTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAgg                                                       >  2:8199513/1‑93 (MQ=255)
                                   ttGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGc                                                      <  1:6759195/93‑1 (MQ=255)
                                    tGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCg                                                     <  2:3331436/93‑1 (MQ=255)
                                     gAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGa                                                    <  1:2799007/93‑1 (MQ=255)
                                      aaCGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGaa                                                   <  1:4973400/93‑1 (MQ=255)
                                      aaCGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGaa                                                   <  1:4802999/93‑1 (MQ=255)
                                       acaaaaaaaaTGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCCGGCCACTACATCGTGTAAACTCCAGTTAATCGCATTGAGGCGAag                                                  <  2:7977556/90‑1 (MQ=255)
                                       aCGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGGGTAAAATCCAGTTAATAGCATTGAGGCCaaa                                                  >  2:4225476/1‑92 (MQ=255)
                                        cGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAaga                                                 <  1:8259553/93‑1 (MQ=255)
                                        cGAAAAAAATGCACCAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAaga                                                 <  1:643800/93‑1 (MQ=255)
                                        agAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAaga                                                 <  1:5832206/92‑1 (MQ=255)
                                         gAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAagag                                                >  2:6883759/1‑93 (MQ=255)
                                             aaaaTGCACTAAGTATGCACTGACTGCATAGTCACTGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAag                                            <  1:6896578/93‑1 (MQ=255)
                                              aaaTGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAaga                                           <  1:4153581/93‑1 (MQ=255)
                                               aaTGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAagag                                          <  1:443229/93‑1 (MQ=255)
                                               aaTGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAagag                                          >  2:1476732/1‑93 (MQ=255)
                                                aTGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAgg                                         >  1:7357016/1‑93 (MQ=255)
                                                aTGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAgg                                         <  2:8457237/93‑1 (MQ=255)
                                                 tGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGa                                        <  2:7392232/93‑1 (MQ=255)
                                                 tGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGa                                        <  1:4306426/93‑1 (MQ=255)
                                                  gCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGaa                                       <  1:6999629/93‑1 (MQ=255)
                                                  gCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGaa                                       >  2:4054854/1‑93 (MQ=255)
                                                   cACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAGCTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAAt                                      <  2:3496536/93‑1 (MQ=255)
                                                    aCTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAAtt                                     >  2:7691425/1‑93 (MQ=255)
                                                    aCTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAAtt                                     >  2:8337550/1‑93 (MQ=255)
                                                    aCTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCAACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAAtt                                     >  2:443018/1‑93 (MQ=255)
                                                     cTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTg                                    >  1:8823166/1‑93 (MQ=255)
                                                     cTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTg                                    >  2:6602113/1‑93 (MQ=255)
                                                      tAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGc                                   <  1:2924085/93‑1 (MQ=255)
                                                           aTGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGAttt                              <  1:6637076/93‑1 (MQ=255)
                                                            tGCACTGACTGCATAGTCACGGCAGCGGGCCACTACATAGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGGATTGCAAtttt                             <  2:7798034/93‑1 (MQ=255)
                                                               actgactgCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCt                          >  1:6652072/1‑93 (MQ=255)
                                                                ctgactgCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTc                         >  2:4245120/1‑93 (MQ=255)
                                                                ctgactgCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTc                         <  2:1797377/93‑1 (MQ=255)
                                                                 tgactgCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGACGAGAAAGAGGAATTGCGATTTTGCTcc                        >  2:6408046/1‑93 (MQ=255)
                                                                 tgactgCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTcc                        <  1:2405847/93‑1 (MQ=255)
                                                                  gactgCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTccc                       <  1:1602443/93‑1 (MQ=255)
                                                                  gactgCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTccc                       <  1:2194655/93‑1 (MQ=255)
                                                                   actgCATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCg                      <  1:5140587/93‑1 (MQ=255)
                                                                       cATAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCt                  >  1:7516260/1‑93 (MQ=255)
                                                                        aTAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCtt                 <  2:5935052/93‑1 (MQ=255)
                                                                        aTAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCtt                 <  2:3154541/93‑1 (MQ=255)
                                                                        aTAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCtt                 >  2:2247885/1‑93 (MQ=255)
                                                                        aTAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCtt                 <  1:2612337/93‑1 (MQ=255)
                                                                        aTAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCtt                 >  2:4025716/1‑93 (MQ=255)
                                                                         tAGTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCttt                <  1:2476471/93‑1 (MQ=255)
                                                                           gTCACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCTTTTc              >  2:3911596/1‑93 (MQ=255)
                                                                             cACGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCTTTTCCt            <  1:3571465/93‑1 (MQ=255)
                                                                              aCGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCTTTTCCTg           >  2:1870536/1‑93 (MQ=255)
                                                                               cGGCAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCTTTTCCTGa          <  2:6307462/93‑1 (MQ=255)
                                                                                  cAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCTTTTCCTGAATg       >  2:4827306/1‑93 (MQ=255)
                                                                                  cAGCGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCTTTTCCTGAATg       <  1:8686042/93‑1 (MQ=255)
                                                                                     cGGGCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCTTTTCCTGAATGACa    >  1:8581682/1‑93 (MQ=255)
                                                                                      gggCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCTTTTCCTGAATGACAt   >  1:7693174/1‑93 (MQ=255)
                                                                                       ggCCACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCTTTTCCTGAATGACATg  <  2:2983921/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGCTTTAACCTGCGCTCCGCGCTCTGCTTGTGAAGTTGAACGAAAAAAATGCACTAAGTATGCACTGACTGCATAGTCACGGCAGCGGGCTACTACATCGCGTAAACTCCAGTTAATCGCATTGAGGCGAAGAGAAAGAGGAATTGCGATTTTGCTCCCGGTCTTTTCCTGAATGACATG  >  NC_000913/1199121‑1199300

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