Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,245,705 T→C G52G (GGT→GGC lamB → phage lambda receptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,245,7050TC99.1% 671.4 / ‑4.8 216G52G (GGT→GGClamBphage lambda receptor protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (84/130);  total (84/132)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 224 total aligned reads displayed.

GGCAATGGCTGTTGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGTGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAGGGCGATAAGAG  >  NC_000913/4245615‑4245794
                                                                                          |                                                                                         
     tGGCTGTTGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCaaa                                                                                    <  2:2195180/93‑1 (MQ=255)
     tGGCTGTTGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCaaa                                                                                    <  1:2442002/93‑1 (MQ=255)
      ggCTGTTGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAg                                                                                   <  1:72063/93‑1 (MQ=255)
      ggCTGTTGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAg                                                                                   <  2:1143102/93‑1 (MQ=255)
      ggCTGTTGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGCAGCGGCGGGGAACAACAGTGTTTCCCCACTACCGGCGCTCAAAg                                                                                   >  1:1139036/1‑93 (MQ=255)
         tGTTGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTaa                                                                                >  2:1089646/1‑93 (MQ=255)
         tGTTGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTaa                                                                                <  1:1832109/93‑1 (MQ=255)
         tGTTGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTaa                                                                                >  2:3234194/1‑93 (MQ=255)
          gTTGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTaaa                                                                               <  1:3114108/93‑1 (MQ=255)
          gTTGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTaaa                                                                               <  2:2768611/93‑1 (MQ=255)
          gTTGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTaaa                                                                               >  2:3190085/1‑93 (MQ=255)
           ttGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAAt                                                                              <  2:3457053/93‑1 (MQ=255)
            tGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATa                                                                             <  1:4331714/93‑1 (MQ=255)
            tGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATa                                                                             <  1:1484636/93‑1 (MQ=255)
            tGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATa                                                                             <  1:1534871/93‑1 (MQ=255)
             gATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATAc                                                                            <  2:3372342/93‑1 (MQ=255)
             gATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATAc                                                                            >  2:446478/1‑93 (MQ=255)
             gATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATAc                                                                            >  2:838218/1‑93 (MQ=255)
               tttCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCg                                                                          <  2:2386552/93‑1 (MQ=255)
                ttCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGt                                                                         >  2:2686676/1‑93 (MQ=255)
                  ccACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCt                                                                       <  1:3077397/93‑1 (MQ=255)
                   cACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCtt                                                                      <  1:3353290/93‑1 (MQ=255)
                   cACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCtt                                                                      <  2:1991481/93‑1 (MQ=255)
                    aCGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTg                                                                     >  1:2246045/1‑93 (MQ=255)
                    aCGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTg                                                                     >  2:4411631/1‑93 (MQ=255)
                       gCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCa                                                                  >  2:4440585/1‑93 (MQ=255)
                        cTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCaa                                                                 <  2:3840089/93‑1 (MQ=255)
                        cTAGGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCaa                                                                 <  2:689129/93‑1 (MQ=255)
                           tGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGa                                                              >  2:3643509/1‑93 (MQ=255)
                            gCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGaa                                                             >  2:42345/1‑93 (MQ=255)
                             cACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAAt                                                            <  2:163899/93‑1 (MQ=255)
                              aCGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAAtg                                                           <  1:4341044/93‑1 (MQ=255)
                              aCGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAAtg                                                           <  1:1947086/93‑1 (MQ=255)
                                gTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAAtgtg                                                         <  2:2811217/93‑1 (MQ=255)
                                gTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGAGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTCGCAACGAAtgtg                                                         >  1:1231219/1‑93 (MQ=255)
                                 ttCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGa                                                        <  2:3863906/93‑1 (MQ=255)
                                    cGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCACCGAATGTGAAAc                                                     <  2:3596989/93‑1 (MQ=255)
                                     ggTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACt                                                    >  1:4083372/1‑93 (MQ=255)
                                     ggTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACt                                                    <  2:2079326/93‑1 (MQ=255)
                                      gTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACtt                                                   <  2:508553/93‑1 (MQ=255)
                                       tATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTa                                                  <  2:3883680/93‑1 (MQ=255)
                                       tATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTa                                                  <  1:768498/93‑1 (MQ=255)
                                        aTTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTAt                                                 <  1:1057930/93‑1 (MQ=255)
                                        aTTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTAt                                                 >  1:3542459/1‑93 (MQ=255)
                                        aTTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTAt                                                 <  2:1507393/93‑1 (MQ=255)
                                          tggttggACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTAACGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATgg                                               >  2:1241531/1‑92 (MQ=255)
                                           ggttggACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCt                                              <  2:455519/93‑1 (MQ=255)
                                           ggttggACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCt                                              >  2:3699311/1‑93 (MQ=255)
                                           ggttggACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCt                                              >  1:525753/1‑93 (MQ=255)
                                                  cAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTaa                                       <  2:4249894/93‑1 (MQ=255)
                                                  cAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTaa                                       <  1:17645/93‑1 (MQ=255)
                                                     gTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAAt                                    >  1:60646/1‑93 (MQ=255)
                                                      tAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAAtt                                   <  1:2075301/93‑1 (MQ=255)
                                                      tAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAAtt                                   <  2:3509704/93‑1 (MQ=255)
                                                        gcggcgGTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTgg                                 >  2:2626941/1‑93 (MQ=255)
                                                         cggcggTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTggg                                <  1:3452422/93‑1 (MQ=255)
                                                          ggcggTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATGCCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGt                               <  2:829841/93‑1 (MQ=255)
                                                           gcggTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTc                              >  1:4493200/1‑93 (MQ=255)
                                                             ggTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAg                            <  2:131761/93‑1 (MQ=255)
                                                             ggTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAg                            <  1:3844776/93‑1 (MQ=255)
                                                              gTGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAgg                           <  1:4018593/93‑1 (MQ=255)
                                                               tGAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGa                          <  2:2491119/93‑1 (MQ=255)
                                                                gAACAACAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGaa                         <  2:1371355/93‑1 (MQ=255)
                                                                 aacaacAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAg                        <  2:3696962/93‑1 (MQ=255)
                                                                 aacaacAGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAg                        <  2:467140/93‑1 (MQ=255)
                                                                  acagcaGTGTTTCCAGACTACCGGTGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAgt                       <  1:1174441/93‑1 (MQ=255)
                                                                  acaacaGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAgg                       >  2:1897082/1‑92 (MQ=255)
                                                                   caacaGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAgtg                      >  2:3963193/1‑93 (MQ=255)
                                                                     acaGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGtgtg                    <  1:4582088/93‑1 (MQ=255)
                                                                     acaGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGtgtg                    <  1:4638484/93‑1 (MQ=255)
                                                                     acaGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGtgtg                    >  1:1753162/1‑93 (MQ=255)
                                                                     acaGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGtgtg                    <  2:4626084/93‑1 (MQ=255)
                                                                      caGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTgg                   >  2:2882634/1‑93 (MQ=255)
                                                                       aGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGa                  >  2:3016274/1‑93 (MQ=255)
                                                                       aGTGTTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTAGGTCAGGAAGTGTGGa                  >  2:1957124/1‑93 (MQ=255)
                                                                        gtgtTTCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGaa                 >  1:3516361/1‑93 (MQ=255)
                                                                           tttCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAaga              <  2:2757161/93‑1 (MQ=255)
                                                                           tttCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAaga              <  1:4457570/93‑1 (MQ=255)
                                                                            ttCCAGACTACCGGCGCTCAAAGTAAATTCCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAagag             <  2:2158638/93‑1 (MQ=255)
                                                                            ttCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAagag             <  1:3691451/93‑1 (MQ=255)
                                                                             tCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAgg            <  2:1759503/93‑1 (MQ=255)
                                                                             tCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAgg            <  1:1189279/93‑1 (MQ=255)
                                                                             tCCAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAgg            <  1:3904478/93‑1 (MQ=255)
                                                                              ccAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAggg           <  2:3060106/93‑1 (MQ=255)
                                                                               cAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAGGGc          >  2:710751/1‑93 (MQ=255)
                                                                               cAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAGGGc          <  1:1021428/93‑1 (MQ=255)
                                                                               cAGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAGGGc          >  1:3599198/1‑93 (MQ=255)
                                                                                aGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTACAATTGGGTCAGGAAGTGTGGAAAGAggcga         >  1:3251429/1‑90 (MQ=255)
                                                                                aGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAGGGCg         >  1:2831632/1‑93 (MQ=255)
                                                                                aGACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATATGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAGGGCg         >  1:2518044/1‑93 (MQ=255)
                                                                                aGACTACCGGCGCTCAAAGGAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTTGGGTTGGCAGAGTGGAAAGAGGGCg         >  2:3777311/1‑93 (MQ=255)
                                                                                 gACTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAGGGCGa        >  2:4070322/1‑93 (MQ=255)
                                                                                  aCTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAGGGCGAt       <  2:4128139/93‑1 (MQ=255)
                                                                                   cTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAGGGCGATa      <  1:229372/93‑1 (MQ=255)
                                                                                   cTACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAGGGCGATa      >  1:2947881/1‑93 (MQ=255)
                                                                                    tACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAGGGCGATaa     >  2:4390822/1‑93 (MQ=255)
                                                                                    tACCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATCAAAATTGGGTCAGGAAGTGTGGAAAGAGGGCGATaa     <  2:4103258/93‑1 (MQ=255)
                                                                                     aCCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAGGGCGATAag    >  1:1279179/1‑93 (MQ=255)
                                                                                     aCCGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAGGGCGATAag    >  2:2308504/1‑93 (MQ=255)
                                                                                       cGGCGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAGGGCGATAagag  >  1:3953791/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGCAATGGCTGTTGATTTCCACGGCTATGCACGTTCCGGTATTGGTTGGACAGGTAGCGGCGGTGAACAACAGTGTTTCCAGACTACCGGTGCTCAAAGTAAATACCGTCTTGGCAACGAATGTGAAACTTATGCTGAATTAAAATTGGGTCAGGAAGTGTGGAAAGAGGGCGATAAGAG  >  NC_000913/4245615‑4245794

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