Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,494,883:1 (A)7→8 100% coding (261/564 nt) idnK → D‑gluconate kinase, thermosensitive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,494,8761.A96.8% 330.3 / 3.2 93coding (254/564 nt)idnKD‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base A (49/41);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 112 total aligned reads displayed.

ATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  NC_000913/4494778‑4494970
                                                                                                   |                                                                                                 
aTCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a                                                                                                 <  1:1156810/100‑1 (MQ=255)
aTCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a                                                                                                 >  2:789098/1‑100 (MQ=255)
aTCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑a                                                                                                 <  1:1898912/100‑1 (MQ=255)
  cTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaa                                                                                               >  1:1967541/1‑100 (MQ=255)
    gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa                                                                                             <  1:621118/100‑1 (MQ=255)
    gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa                                                                                             >  1:1270466/1‑100 (MQ=255)
    gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaa                                                                                             >  1:2042600/1‑100 (MQ=255)
     atgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaa                                                                                            <  1:705989/100‑1 (MQ=255)
      tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaa                                                                                           >  1:1421347/1‑100 (MQ=255)
      tgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaa                                                                                           <  1:1860438/100‑1 (MQ=255)
        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                         <  1:541925/100‑3 (MQ=255)
        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaac                                                                                         <  1:1340505/100‑3 (MQ=255)
        aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATAAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaa                                                                                          >  2:1064884/1‑98 (MQ=255)
         aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaaca                                                                                        >  1:993417/1‑97 (MQ=255)
         aGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATAATT‑‑aaaaaaaac                                                                                         <  2:836899/99‑3 (MQ=255)
          gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaacag                                                                                       >  1:1536527/1‑96 (MQ=255)
          gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑aaaaaaaacag                                                                                       <  1:455070/100‑5 (MQ=255)
                  tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGt                                                                                >  2:905236/1‑98 (MQ=255)
                  tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGa                                                                              <  1:853410/100‑1 (MQ=255)
                  tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGa                                                                              >  2:288561/1‑100 (MQ=255)
                  tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGa                                                                              <  1:417676/100‑1 (MQ=255)
                       ggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAATAATGAAACAGGATTTATTGTCTGTTAATCATTAAAAAAAACA‑‑GTATCGTGATAttt                                                                         <  2:238401/100‑1 (MQ=255)
                       ggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATAttt                                                                         >  2:1085183/1‑100 (MQ=255)
                        gctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATAtttt                                                                        >  1:968118/1‑100 (MQ=255)
                          tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTAc                                                                      <  1:1185510/100‑1 (MQ=255)
                          tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTAc                                                                      <  1:121889/100‑1 (MQ=255)
                          tggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTAc                                                                      >  1:2099774/1‑100 (MQ=255)
                           ggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACg                                                                     >  2:440294/1‑100 (MQ=255)
                            gAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTAc                                                                      >  2:1444175/1‑98 (MQ=255)
                             aaCGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCa                                                                   <  1:1980192/100‑1 (MQ=255)
                               cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAg                                                                 <  1:846580/100‑1 (MQ=255)
                                gCTTAAATGATGCTTCATACAGGCTTTATAAAAAGAATGAAACATGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAg                                                                 <  2:1160363/99‑1 (MQ=255)
                                 cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTATGCAAgg                                                                >  2:946054/1‑99 (MQ=255)
                                 cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAggg                                                               >  2:355561/1‑100 (MQ=255)
                                 cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAggg                                                               <  2:282934/100‑1 (MQ=255)
                                     aaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGc                                                           >  1:1398274/1‑100 (MQ=255)
                                       tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGccc                                                         >  2:2270854/1‑100 (MQ=255)
                                        gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGcccc                                                        <  1:2058888/100‑1 (MQ=255)
                                        gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGcccc                                                        <  1:1200291/100‑1 (MQ=255)
                                        gatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGcccc                                                        >  2:130348/1‑100 (MQ=255)
                                         atgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGccccc                                                       >  1:1865566/1‑100 (MQ=255)
                                          tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                      <  2:993412/100‑1 (MQ=255)
                                           gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                      >  1:553398/1‑99 (MQ=255)
                                           gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                     >  2:1197726/1‑100 (MQ=255)
                                           gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGGAGccccc                                                       >  1:658678/1‑98 (MQ=255)
                                            cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                    <  2:374958/100‑1 (MQ=255)
                                            cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                    <  1:51366/100‑1 (MQ=255)
                                             ttCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                    <  1:228/100‑1 (MQ=255)
                                              tCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                                 <  2:881738/100‑1 (MQ=255)
                                               cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGtt                                                  <  1:2220847/99‑1 (MQ=255)
                                                 tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                               >  1:2054138/1‑100 (MQ=255)
                                                 tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                               <  2:1640722/100‑1 (MQ=255)
                                                   cAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAttt                                             >  2:1597796/1‑100 (MQ=255)
                                                   cAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAttt                                             >  2:1596998/1‑100 (MQ=255)
                                                    aGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTc                                            <  2:1775556/100‑1 (MQ=255)
                                                    aGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAAAAAGGAATCTTATTTTTTTACCCAAGGTAAGCCCCCATGTTCAttt                                             >  2:1627425/1‑100 (MQ=255)
                                                        tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctct                                        <  2:2087907/100‑1 (MQ=255)
                                                         ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                       >  1:1602172/1‑100 (MQ=255)
                                                         ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                       >  1:810543/1‑100 (MQ=255)
                                                          tataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                      <  1:1761887/100‑1 (MQ=255)
                                                          tataAAAAGAATGAAACAAGATTTATTGTCTGTTCAGCATTAAAAAAAACA‑‑GTATCGTAATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCct                                          <  2:1908747/96‑1 (MQ=255)
                                                           ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAATGGTCGCCCCCATGTTCATTTCCTCTg                                       >  2:1106518/1‑98 (MQ=255)
                                                           ataAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGt                                     >  2:273969/1‑100 (MQ=255)
                                                             aaaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                    <  2:161606/99‑1 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGTGTAGCCCCCATGTTCATCTCCTCTGGTTAg                                  >  2:1707284/1‑100 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                  >  2:829170/1‑100 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                  >  1:2264144/1‑100 (MQ=255)
                                                              aaaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCCGGTTAg                                  >  1:2519/1‑100 (MQ=255)
                                                               aaaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                 >  2:402536/1‑100 (MQ=255)
                                                               aaaGAATGAAACAGGATTTACTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCTTCTGGTTata                                 >  2:353312/1‑98 (MQ=255)
                                                                aaGAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGAt                                <  2:1769106/100‑1 (MQ=255)
                                                                  gAATGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCACAATGTTCATTTACTCTGGTTAGATgg                              >  2:1168284/1‑100 (MQ=255)
                                                                   aaTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                             <  1:2185436/100‑1 (MQ=255)
                                                                   aaTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGc                             <  1:1576005/100‑1 (MQ=255)
                                                                    aTGAAACAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCg                            <  1:2251105/100‑1 (MQ=255)
                                                                      gAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                         <  1:384825/100‑1 (MQ=255)
                                                                       aaaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                         <  2:968113/100‑1 (MQ=255)
                                                                       aaaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                         >  1:187885/1‑100 (MQ=255)
                                                                        aaCAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCCTTTCCTCTGGTTAGATGGCg                            >  2:1757351/1‑96 (MQ=255)
                                                                         aCAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTAt                       <  2:1084281/100‑1 (MQ=255)
                                                                         aCAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTAt                       <  1:1986628/100‑1 (MQ=255)
                                                                          cAGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                      <  2:2054122/100‑1 (MQ=255)
                                                                           aGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                     <  2:157618/100‑1 (MQ=255)
                                                                           aGGATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGa                     <  2:702607/100‑1 (MQ=255)
                                                                             gATTTATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaa                    >  2:1488763/1‑99 (MQ=255)
                                                                              aTTTATTATCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCTAGGGTAGCCCCCATATTCATTTCCTCTGGTTAGATGGCGCCTATgta                    >  2:639843/1‑96 (MQ=255)
                                                                                ttATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTa                >  2:1760169/1‑100 (MQ=255)
                                                                                 tATTGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAt               <  1:1760182/100‑1 (MQ=255)
                                                                                    tGTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtc             >  2:59855/1‑99 (MQ=255)
                                                                                     gTCTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc           >  2:67575/1‑100 (MQ=255)
                                                                                       cTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgc         >  1:21326/1‑100 (MQ=255)
                                                                                       cTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgc         >  2:937948/1‑100 (MQ=255)
                                                                                       cTGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgc         >  1:818098/1‑100 (MQ=255)
                                                                                        tGTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcg        >  1:782842/1‑100 (MQ=255)
                                                                                         gTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTCTGAAACTATTCTcgcgc       >  2:674039/1‑100 (MQ=255)
                                                                                         gTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc       <  1:2152164/100‑1 (MQ=255)
                                                                                         gTTCATCATTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGACACCATGTTTATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCgcgcgc       >  2:2261913/1‑100 (MQ=255)
                                                                                            catcatTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGaa    >  1:2276173/1‑100 (MQ=255)
                                                                                             atcatTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt   <  1:789101/100‑1 (MQ=255)
                                                                                              tcatTAAAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATg  <  2:1185294/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
ATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCATCATT‑‑AAAAAAACA‑‑GTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATG  >  NC_000913/4494778‑4494970

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