Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,650,355 T→C pseudogene (781/1158 nt) intQ → pseudogene, Qin prophage; predicted defective integrase;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,650,3550TC99.9% 2869.5 / NA 686pseudogene (781/1158 nt)intQpseudogene, Qin prophage; predicted defective integrase;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (265/419);  total (266/420)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 695 total aligned reads displayed.

CCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGT  >  NC_000913/1650257‑1650452
                                                                                                  |                                                                                                 
 nnntCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCa                                                                                                <  2:9864753/98‑1 (MQ=255)
 cAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCa                                                                                                <  1:8583222/101‑1 (MQ=255)
 cAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCa                                                                                                >  2:566613/1‑101 (MQ=255)
 cAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCa                                                                                                >  2:3530431/1‑101 (MQ=255)
 cAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCa                                                                                                >  2:1706509/1‑101 (MQ=255)
   nnngggTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCata                                                                                              <  2:11015730/98‑3 (MQ=255)
   aTCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCaga                                                                                              <  1:2620176/101‑1 (MQ=255)
    tCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagag                                                                                             <  1:5999919/101‑1 (MQ=255)
     nnngTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagaga                                                                                            <  2:14045105/98‑1 (MQ=255)
     cGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagaga                                                                                            <  1:14685698/101‑1 (MQ=255)
      gggTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagagag                                                                                           >  1:14149811/1‑101 (MQ=255)
      gggTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagagag                                                                                           >  2:3738837/1‑101 (MQ=255)
       ggTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGt                                                                                          <  1:2895430/101‑1 (MQ=255)
        gTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTa                                                                                         <  1:4262893/101‑1 (MQ=255)
         nnnttcatCTTATTAAGCCAGCAATCGACGCCCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTAt                                                                                        <  2:8875728/98‑1 (MQ=255)
          nnntcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATg                                                                                       <  2:8246355/98‑1 (MQ=255)
           nnncatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATgg                                                                                      <  2:14674479/98‑1 (MQ=255)
            ttcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGc                                                                                     <  1:6255930/101‑1 (MQ=255)
              catCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcag                                                                                   <  1:3354923/101‑1 (MQ=255)
                tCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaa                                                                                 <  1:3182422/101‑1 (MQ=255)
                 nnnattaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaac                                                                                <  2:7890960/98‑1 (MQ=255)
                 cTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaac                                                                                >  1:3555690/1‑101 (MQ=255)
                    attaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaacaga                                                                             <  1:6556478/101‑1 (MQ=255)
                     ttaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaacagaa                                                                            >  1:5052912/1‑101 (MQ=255)
                      taAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaa                                                                           >  1:715572/1‑101 (MQ=255)
                      taAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaa                                                                           <  2:2116927/101‑1 (MQ=255)
                       nnnccAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaaa                                                                          <  2:10304962/98‑1 (MQ=255)
                       aaGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaaa                                                                          <  1:13845988/101‑1 (MQ=255)
                         gCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaannn                                                                        >  2:11889204/1‑98 (MQ=255)
                          nnngCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACa                                                                       <  2:14340931/98‑1 (MQ=255)
                          ccAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACa                                                                       <  1:7519255/101‑1 (MQ=255)
                           cAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaa                                                                      >  1:10231636/1‑101 (MQ=255)
                            aGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAAcnnn                                                                     >  2:9465705/1‑98 (MQ=255)
                            aGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaa                                                                     <  1:11017442/101‑1 (MQ=255)
                            aGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaa                                                                     <  2:6345998/101‑1 (MQ=255)
                            aGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaa                                                                     >  1:2360534/1‑101 (MQ=255)
                               aaTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAAt                                                                  >  1:5196696/1‑101 (MQ=255)
                                aTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAAAGTATGGCAGAACAGAAAAACaaaannn                                                                 >  2:12925527/1‑98 (MQ=255)
                                 tCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGc                                                                >  2:740084/1‑101 (MQ=255)
                                   nnngCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACACAAATGccc                                                              <  2:10826337/98‑3 (MQ=255)
                                   gACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCAc                                                              <  2:4276572/101‑1 (MQ=255)
                                     nnntctcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCt                                                            <  2:12670052/98‑1 (MQ=255)
                                      gatctcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCtt                                                           >  2:3195175/3‑101 (MQ=255)
                                      gCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCtt                                                           >  1:6900192/1‑101 (MQ=255)
                                        tctcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTg                                                         <  1:10027498/101‑1 (MQ=255)
                                         ctcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGt                                                        <  1:1780985/101‑1 (MQ=255)
                                          tcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGtt                                                       >  1:12649382/1‑101 (MQ=255)
                                           ccGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttt                                                      <  1:8029474/101‑1 (MQ=255)
                                            cGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGtttt                                                     <  1:3814824/101‑1 (MQ=255)
                                            cGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGtttt                                                     <  1:9936011/101‑1 (MQ=255)
                                             ggAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGGTGGTCACCTCAGAGAGTATGGCAGAAAAGAAAGACAAAAAGGCAGCATTGttttt                                                    >  2:657037/1‑101 (MQ=255)
                                             ggAGTCAGATGACATAAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttttt                                                    >  1:13410210/1‑101 (MQ=255)
                                              gAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTc                                                   <  1:4589747/101‑1 (MQ=255)
                                               aGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCa                                                  <  2:543799/101‑1 (MQ=255)
                                                 tCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTcnnn                                                >  2:9238428/1‑98 (MQ=255)
                                                 tCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAAc                                                >  2:2661197/1‑101 (MQ=255)
                                                 tCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTCTGTTTTTCAAc                                                <  2:5936103/101‑1 (MQ=255)
                                                  taGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCannn                                               >  2:11706021/2‑98 (MQ=255)
                                                   aGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCt                                              >  1:13046631/1‑101 (MQ=255)
                                                   aGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCt                                              >  1:11222225/1‑101 (MQ=255)
                                                   aGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCt                                              >  1:2011932/1‑101 (MQ=255)
                                                     aTGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGa                                            <  1:11809228/101‑1 (MQ=255)
                                                      tGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGaa                                           <  1:9670866/101‑1 (MQ=255)
                                                           ttAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTc                                      >  2:5030441/1‑101 (MQ=255)
                                                           nnnaCGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTc                                      <  2:13508051/98‑1 (MQ=255)
                                                            tAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCa                                     >  2:22580/1‑101 (MQ=255)
                                                             nnngagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAg                                    <  2:7560866/98‑1 (MQ=255)
                                                             aaCGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAg                                    <  2:4127254/101‑1 (MQ=255)
                                                             aaCGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAg                                    <  1:5399657/101‑1 (MQ=255)
                                                             aaCGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAg                                    <  1:9407046/101‑1 (MQ=255)
                                                                gagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCga                                 <  1:4803235/101‑1 (MQ=255)
                                                                gagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCga                                 >  1:8843267/1‑101 (MQ=255)
                                                                 agaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCgag                                >  1:5562804/1‑101 (MQ=255)
                                                                  gaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCgaga                               >  1:8269277/1‑101 (MQ=255)
                                                                     tGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTa                            <  2:5219937/101‑1 (MQ=255)
                                                                     tGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTa                            <  2:715572/101‑1 (MQ=255)
                                                                     nnngCAAAGAGCATATCATTGATGTTCAACTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGagagaa                            <  2:10568613/98‑3 (MQ=255)
                                                                      gAGCAAAGAGCATATCATTGATGTTCACCTCAGAGCATATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaa                           <  1:11985377/101‑1 (MQ=255)
                                                                      gAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaa                           >  1:9336793/1‑101 (MQ=255)
                                                                       nnnaaaGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaa                          <  2:12366393/98‑1 (MQ=255)
                                                                       aGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGtnnn                          >  2:12248232/1‑98 (MQ=255)
                                                                       aGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaa                          <  1:14341492/101‑1 (MQ=255)
                                                                        gCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTannn                         >  2:8490827/1‑98 (MQ=255)
                                                                        gCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaa                         <  1:4126330/101‑1 (MQ=255)
                                                                          aaaGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaaaa                       >  2:3356264/1‑101 (MQ=255)
                                                                           aaGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAAt                      >  2:4709843/1‑101 (MQ=255)
                                                                           aaGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAAt                      <  2:6843647/101‑1 (MQ=255)
                                                                            agagCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAAtt                     <  1:11616396/101‑1 (MQ=255)
                                                                               gCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCGCCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATg                  <  1:10700607/101‑1 (MQ=255)
                                                                                nnnatCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATgg                 <  2:11277062/98‑1 (MQ=255)
                                                                                   atCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGa              <  2:3826765/101‑1 (MQ=255)
                                                                                     cATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAcc            <  2:2949924/101‑1 (MQ=255)
                                                                                     cATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAcc            >  1:13602363/1‑101 (MQ=255)
                                                                                      aTTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCa           <  1:901776/101‑1 (MQ=255)
                                                                                        tGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCAtt         <  1:6730964/101‑1 (MQ=255)
                                                                                          aTGTTTACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCAtttt       <  1:12334739/101‑1 (MQ=255)
                                                                                           tGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTa      >  1:1560219/1‑101 (MQ=255)
                                                                                           tGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGGAAAAAATTATGGGGACCATTTTa      >  1:5781898/1‑101 (MQ=255)
                                                                                             ttCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTAcc    <  2:908966/101‑1 (MQ=255)
                                                                                              tCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATTCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCg   >  1:7835974/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
CCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGT  >  NC_000913/1650257‑1650452

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: