Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 529,036 A→G 29.9% pseudogene (397/491 nt) ybbD → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913529,0360AG29.9% 153.7 / 98.1 117pseudogene (397/491 nt)ybbDpseudogene
Reads supporting (aligned to +/- strand):  ref base A (39/43);  new base G (20/15);  total (59/58)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

Only 100 of 119 total aligned reads displayed.

GTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCT  >  NC_000913/528939‑529135
                                                                                                 |                                                                                                    
gTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATg                                                                                                    >  2:111732/1‑100 (MQ=255)
 tGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTgg                                                                                                   <  2:915435/100‑1 (MQ=255)
 tGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATgg                                                                                                   <  2:503995/100‑1 (MQ=255)
   aGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGgtggtg                                                                                                 <  1:1328599/100‑1 (MQ=255)
   aGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGgtg                                                                                                 <  2:596997/100‑1 (MQ=255)
   aGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGgtg                                                                                                 <  2:1810830/100‑1 (MQ=255)
   aGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGgtg                                                                                                 <  2:2981830/100‑1 (MQ=255)
    gggCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGgtgt                                                                                                >  2:1106359/1‑100 (MQ=255)
     ggCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTa                                                                                               >  2:3111058/1‑100 (MQ=255)
      gCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTaa                                                                                              >  2:983885/1‑100 (MQ=255)
      gCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTaa                                                                                              <  2:46177/100‑1 (MQ=255)
      gCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTa                                                                                               <  1:557915/99‑1 (MQ=255)
        cAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAAc                                                                                            <  2:542971/100‑1 (MQ=255)
         aGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAAcc                                                                                           <  1:779011/100‑1 (MQ=255)
            gAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑ttt                                                                                       >  1:1897931/1‑100 (MQ=255)
               cAAGTGGCGAAAATCGTATCAAAGAATGATTCAGATACAAAAAATGTGTGGTGTAAATATGGTAAGCTACCAGGGCAAGGGGATGGTGTAAACC‑tttttt                                                                                    <  2:619127/100‑1 (MQ=255)
                 aGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑tttttttt                                                                                  >  2:2877202/1‑100 (MQ=255)
                 aGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑tttttttt                                                                                  >  1:2461068/1‑100 (MQ=255)
                 aGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑tttttttt                                                                                  >  1:768924/1‑100 (MQ=255)
                 aGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑tttttttt                                                                                  >  2:964187/1‑100 (MQ=255)
                 aGTGGCGAAAATCGGATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑tttttttt                                                                                  >  1:3239210/1‑100 (MQ=255)
                  gTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTttg                                                                                 <  1:3139707/100‑1 (MQ=255)
                  gTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTttg                                                                                 <  1:1411782/100‑1 (MQ=255)
                     gCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTttgttg                                                                              >  1:1539414/1‑100 (MQ=255)
                     gCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTttgttg                                                                              <  2:952589/100‑1 (MQ=255)
                      cGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTgg                                                                             >  1:3424252/1‑100 (MQ=255)
                      cGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTgg                                                                             >  2:3095119/1‑100 (MQ=255)
                      cGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTgg                                                                             >  2:2360126/1‑100 (MQ=255)
                      cGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTgg                                                                             >  1:470052/1‑100 (MQ=255)
                      aaaaaaTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTttgtt                                                                               <  1:998945/96‑1 (MQ=255)
                              gTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACCTTTTTTTTTGTTGGTGAAAt                                                                       >  1:3098588/1‑99 (MQ=255)
                               tATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTaa                                                                    <  1:3207095/100‑1 (MQ=255)
                                aTCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAAt                                                                   <  2:2825616/100‑1 (MQ=255)
                                  cAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGt                                                                 >  1:220985/1‑100 (MQ=255)
                                   aaaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGtt                                                                >  2:1890526/1‑100 (MQ=255)
                                   aaaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGtt                                                                >  1:857819/1‑100 (MQ=255)
                                   aaaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGtt                                                                <  1:2526669/100‑1 (MQ=255)
                                   aaaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGtt                                                                <  2:1644579/100‑1 (MQ=255)
                                       aaTGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGc                                                            <  2:2225066/100‑1 (MQ=255)
                                       aaTGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGc                                                            <  1:1244478/100‑1 (MQ=255)
                                            tCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCattat                                                       >  2:1426360/1‑100 (MQ=255)
                                               aGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGTTTTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTAtttt                                                    <  1:1193013/100‑1 (MQ=255)
                                               aGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTAtttt                                                    >  2:3057405/1‑100 (MQ=255)
                                               aGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTAtttt                                                    <  1:1102299/100‑1 (MQ=255)
                                               aGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTAtttt                                                    <  2:306460/100‑1 (MQ=255)
                                                gATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTta                                                   <  2:2181677/100‑1 (MQ=255)
                                                gATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTta                                                   <  1:914062/100‑1 (MQ=255)
                                                 aTAAAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTta                                                   <  1:2503470/99‑1 (MQ=255)
                                                  tACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTtata                                                 <  2:1238483/100‑1 (MQ=255)
                                                   aCAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATaa                                                <  1:111732/100‑1 (MQ=255)
                                                    cAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAAc                                               <  1:3257789/100‑1 (MQ=255)
                                                    cAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAAc                                               >  1:1695212/1‑100 (MQ=255)
                                                     aaaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACa                                              >  1:1591653/1‑100 (MQ=255)
                                                     aaaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTAGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACa                                              >  1:1921580/1‑100 (MQ=255)
                                                     aaaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACa                                              >  2:3462023/1‑100 (MQ=255)
                                                      aaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACaa                                             <  2:1897947/100‑1 (MQ=255)
                                                      aaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTTTTTGTTTAATTAATGTTACGCATTATTTtata                                                 >  1:2245671/1‑96 (MQ=255)
                                                      aaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTTGTGAAATTAATGTTACGCATTATTTTATAACa                                              >  2:1210793/1‑99 (MQ=255)
                                                      aaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACaa                                             >  2:1251078/1‑100 (MQ=255)
                                                      aaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACaa                                             >  1:1804142/1‑100 (MQ=255)
                                                      aaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACaa                                             >  2:855253/1‑100 (MQ=255)
                                                           aGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATAtt                                        <  1:1226739/100‑1 (MQ=255)
                                                             tgtgGTGTAAATATGGTAAGATGCCAGGGCACGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTgg                                      <  1:2328719/100‑1 (MQ=255)
                                                             tgtgGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTgg                                      <  1:2975462/100‑1 (MQ=255)
                                                              gtggtgTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGa                                     <  1:3011318/100‑1 (MQ=255)
                                                                ggtgTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGc                                   >  2:367323/1‑100 (MQ=255)
                                                                 gtgtAAATATGGTAAGATACCAGGGCAAGGGGATGGTGT‑AACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTg                                 >  1:1948395/1‑100 (MQ=255)
                                                                  tgtAAATATGGTAAGATACCAGGGCAAGGGGATAGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTg                                 <  2:241324/100‑1 (MQ=255)
                                                                    tAAATATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGa                               >  1:2797913/1‑100 (MQ=255)
                                                                     aaaTATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGAt                              >  1:2148020/1‑100 (MQ=255)
                                                                      aaTATGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGAtt                             >  1:1030589/1‑100 (MQ=255)
                                                                      aaTATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGAtt                             >  2:1124155/1‑100 (MQ=255)
                                                                         atGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGcc                          <  1:418355/100‑1 (MQ=255)
                                                                          tGGTAAGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCt                         >  2:283781/1‑100 (MQ=255)
                                                                          tGGTAAGATACCAGGGCAAGGGGGTGATGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCt                         >  2:2111446/1‑100 (MQ=255)
                                                                          tGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCt                         <  2:1212286/100‑1 (MQ=255)
                                                                          cgatctgATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCt                         <  1:1251104/94‑1 (MQ=255)
                                                                              aaGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATg                     <  1:1780755/100‑1 (MQ=255)
                                                                              aaGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATg                     <  1:298850/100‑1 (MQ=255)
                                                                               aGATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTCTTTTATAACAAATATTGGAGCTGGATTGCCTGATGc                    <  2:2443080/100‑1 (MQ=255)
                                                                                gATACCAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCt                   >  1:1793313/1‑100 (MQ=255)
                                                                                gATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCtagatcg                   >  2:1251115/1‑94 (MQ=255)
                                                                                  tACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtg                 >  1:2631222/1‑100 (MQ=255)
                                                                                    ccAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTTGTGAAATTAATGTTACGCATTATTTTATAACAACTATTGGAGCTGGATTGCCTGATGCTtgtg               >  1:2854296/1‑100 (MQ=255)
                                                                                    ccAGGGCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtgtg               <  2:1318989/100‑1 (MQ=255)
                                                                                    ccAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtgtg               >  2:1486353/1‑100 (MQ=255)
                                                                                    ccAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtgtg               >  2:677735/1‑100 (MQ=255)
                                                                                     cAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGc              <  2:470056/100‑1 (MQ=255)
                                                                                     cAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGc              <  1:1106350/100‑1 (MQ=255)
                                                                                       gggCAAGGGGGTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag            <  2:3369790/100‑1 (MQ=255)
                                                                                         gCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCagag          <  2:2585530/100‑1 (MQ=255)
                                                                                         gCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGATGAAATTAATGTTACGCACTATTTTATAACAAATATTGGTTCTCTATTGCCTGTTGcgt                  >  1:2680755/1‑90 (MQ=255)
                                                                                          cAAGGGGATGGTGTAAACC‑TTTTTTTTTTTTTTTTTTTTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGt         >  1:327682/1‑100 (MQ=255)
                                                                                          cAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGt         <  1:2360107/100‑1 (MQ=255)
                                                                                             ggggATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAAt      <  1:1330166/100‑1 (MQ=255)
                                                                                              gggATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAAtt     >  1:2623832/1‑100 (MQ=255)
                                                                                               ggATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTg    >  1:1072453/1‑100 (MQ=255)
                                                                                               ggATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTg    <  2:2696031/100‑1 (MQ=255)
                                                                                               ggATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTg    >  2:3021681/1‑100 (MQ=255)
                                                                                                 aTGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCt  <  1:576496/100‑1 (MQ=255)
                                                                                                 |                                                                                                    
GTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACC‑TTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCT  >  NC_000913/528939‑529135

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