Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,279,197 A→C R145S (AGA→AGC agaA → putative N‑acetylgalactosamine‑6‑phosphate deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,279,1970AC99.6% 794.8 / NA 248R145S (AGA→AGCagaAputative N‑acetylgalactosamine‑6‑phosphate deacetylase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (125/122);  total (125/123)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

GGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGAGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGG  >  NC_000913/3279107‑3279286
                                                                                          |                                                                                         
 gAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTc                                                                                        >  2:2409121/1‑93 (MQ=255)
 gAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTc                                                                                        >  1:1008863/1‑93 (MQ=255)
   aGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGcgtcgt                                                                                      <  1:5029403/93‑1 (MQ=255)
    gctatccatATGGCGTCTCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGtt                                                                                     <  2:2080458/90‑1 (MQ=255)
    gCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGtt                                                                                     >  1:3644093/1‑93 (MQ=255)
    gCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGtt                                                                                     <  2:5489780/93‑1 (MQ=255)
     ccatccatATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTg                                                                                    >  2:4447621/1‑93 (MQ=255)
      catccatATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTgc                                                                                   <  1:4669516/93‑1 (MQ=255)
       atccatATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTgcg                                                                                  >  1:1088660/1‑93 (MQ=255)
       atccatATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTgcg                                                                                  >  1:5813041/1‑93 (MQ=255)
        tccatATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTgcgc                                                                                 >  2:223218/1‑93 (MQ=255)
        tccatATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTgcgc                                                                                 >  1:1587214/1‑93 (MQ=255)
        tccatATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTgcgc                                                                                 <  1:5606367/93‑1 (MQ=255)
         ccatATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCt                                                                                <  1:5215464/93‑1 (MQ=255)
           atatGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTgg                                                                              <  1:2426927/93‑1 (MQ=255)
           atatGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGCTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTgg                                                                              >  1:3384163/1‑93 (MQ=255)
            tatGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGa                                                                             >  1:2979551/1‑93 (MQ=255)
            tatGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGa                                                                             >  2:529979/1‑93 (MQ=255)
             atGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGAt                                                                            >  1:2154663/1‑93 (MQ=255)
             atGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGAt                                                                            <  2:1268571/93‑1 (MQ=255)
             atGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGAt                                                                            <  2:2214999/93‑1 (MQ=255)
              tGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATa                                                                           >  1:5461410/1‑93 (MQ=255)
              tGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATa                                                                           <  1:6006804/93‑1 (MQ=255)
               ggCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAg                                                                          >  2:517550/1‑93 (MQ=255)
                 cGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCg                                                                        <  1:3095644/93‑1 (MQ=255)
                  gTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCgg                                                                       <  1:2485990/93‑1 (MQ=255)
                  gTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCgg                                                                       <  2:5381468/93‑1 (MQ=255)
                   tCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCggg                                                                      >  2:2956724/1‑93 (MQ=255)
                   tCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCggg                                                                      >  1:1037280/1‑93 (MQ=255)
                   tCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCggg                                                                      >  2:5712847/1‑93 (MQ=255)
                     gCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCt                                                                    >  1:836161/1‑93 (MQ=255)
                      cTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTa                                                                   >  2:1368156/1‑93 (MQ=255)
                           tCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACAtgt                                                              >  1:261112/1‑93 (MQ=255)
                             cGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGc                                                            >  1:6233610/1‑93 (MQ=255)
                             cGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGc                                                            >  1:6125264/1‑93 (MQ=255)
                                cgcgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGcaac                                                         <  2:1246708/93‑1 (MQ=255)
                                 gcgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGcaaca                                                        <  1:3138327/93‑1 (MQ=255)
                                 gcgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGcaaca                                                        >  1:2947186/1‑93 (MQ=255)
                                    aaTGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGcaacaact                                                     >  1:536867/1‑91 (MQ=255)
                                     aTGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATc                                                    >  2:1798441/1‑93 (MQ=255)
                                      tGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCt                                                   >  1:4577412/1‑93 (MQ=255)
                                      tGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCt                                                   <  2:71534/93‑1 (MQ=255)
                                       gCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTg                                                  <  1:3755998/93‑1 (MQ=255)
                                         tGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGa                                                >  1:4179354/1‑93 (MQ=255)
                                          gggTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGAt                                               <  2:494332/93‑1 (MQ=255)
                                           ggTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGAtt                                              >  2:4200954/1‑93 (MQ=255)
                                           ggTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGAtt                                              <  2:1469957/93‑1 (MQ=255)
                                             tgtTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCa                                            <  2:168018/93‑1 (MQ=255)
                                              gtTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAg                                           <  2:5204209/93‑1 (MQ=255)
                                               ttGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAgg                                          >  2:618988/1‑93 (MQ=255)
                                               ttGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAgg                                          >  2:2128165/1‑93 (MQ=255)
                                                 gATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGt                                        >  1:3489052/1‑93 (MQ=255)
                                                  aTGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTc                                       >  2:3161601/1‑93 (MQ=255)
                                                    ggTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCaa                                     >  1:732994/1‑93 (MQ=255)
                                                    ggTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCaa                                     >  2:614510/1‑93 (MQ=255)
                                                    ggTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCaa                                     >  1:5955594/1‑93 (MQ=255)
                                                      tgtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAAtt                                   <  2:351179/93‑1 (MQ=255)
                                                      tgtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAAtt                                   <  2:3298992/93‑1 (MQ=255)
                                                      tgtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAAtt                                   <  2:655717/93‑1 (MQ=255)
                                                       gtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTa                                  <  1:6068860/93‑1 (MQ=255)
                                                       gtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTa                                  >  1:1381532/1‑93 (MQ=255)
                                                       gtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTa                                  <  1:1719120/93‑1 (MQ=255)
                                                       gtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTa                                  <  1:2887437/93‑1 (MQ=255)
                                                       gtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTa                                  >  2:2614008/1‑93 (MQ=255)
                                                         tCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGc                                >  2:1430760/1‑93 (MQ=255)
                                                          cTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCt                               >  1:208660/1‑93 (MQ=255)
                                                           tGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCtt                              >  2:4111341/1‑93 (MQ=255)
                                                           tGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCtt                              <  2:1003573/93‑1 (MQ=255)
                                                           tGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCtt                              >  2:1649665/1‑93 (MQ=255)
                                                            gggATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTc                             <  1:2250879/93‑1 (MQ=255)
                                                             ggATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCg                            >  1:3989619/1‑93 (MQ=255)
                                                                tCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGttt                         <  1:4454560/93‑1 (MQ=255)
                                                                   cTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGa                      <  2:289921/93‑1 (MQ=255)
                                                                   cTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGa                      <  1:3582595/93‑1 (MQ=255)
                                                                   cTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGa                      >  2:3895886/1‑93 (MQ=255)
                                                                   cTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGa                      >  2:3621919/1‑93 (MQ=255)
                                                                   cTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGa                      <  2:2567760/93‑1 (MQ=255)
                                                                    ttAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGAt                     >  2:406280/1‑93 (MQ=255)
                                                                    ttAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGAt                     <  2:5412216/93‑1 (MQ=255)
                                                                     tAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATa                    >  2:5950573/1‑93 (MQ=255)
                                                                     tAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATa                    >  1:1395163/1‑93 (MQ=255)
                                                                      aaaCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAg                   >  1:6203042/1‑93 (MQ=255)
                                                                      aaaCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAg                   <  2:4732032/93‑1 (MQ=255)
                                                                       aaCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGt                  >  1:4887922/1‑93 (MQ=255)
                                                                       aaCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGt                  <  1:4842790/93‑1 (MQ=255)
                                                                        aCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGtt                 <  2:934802/93‑1 (MQ=255)
                                                                           gggCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGc              <  2:5907076/93‑1 (MQ=255)
                                                                            ggCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCt             <  2:55943/93‑1 (MQ=255)
                                                                               aaaCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCt          >  2:5070006/1‑93 (MQ=255)
                                                                               aaaCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCt          >  2:1404718/1‑93 (MQ=255)
                                                                                aaCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCtt         >  1:5722927/1‑93 (MQ=255)
                                                                                  cgcgCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTa       <  2:5335231/93‑1 (MQ=255)
                                                                                   gcgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTAt      <  1:2723598/93‑1 (MQ=255)
                                                                                     gcCAGCGTCGTTGCTCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTg    <  2:2773159/93‑1 (MQ=255)
                                                                                     gcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTg    >  2:3010988/1‑93 (MQ=255)
                                                                                     gcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTg    <  1:1167650/93‑1 (MQ=255)
                                                                                      ccAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTgg   >  1:422651/1‑93 (MQ=255)
                                                                                      ccAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTgg   >  1:1644860/1‑93 (MQ=255)
                                                                                      ccAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTgg   <  1:1304807/93‑1 (MQ=255)
                                                                                       cAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTggg  >  2:4568020/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGAGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGG  >  NC_000913/3279107‑3279286

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