Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 321,039 G→A D70N (GAC→AAC)  ykgE → putative dehydrogenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913321,0390GA99.3% 521.6 / ‑5.1 146D70N (GAC→AAC) ykgEputative dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (66/79);  total (66/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 155 total aligned reads displayed.

GGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACGACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCATC  >  NC_000913/320949‑321127
                                                                                          |                                                                                        
ggATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATaacaac                                                                                        >  1:1417370/1‑93 (MQ=255)
 gATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAacaacg                                                                                       >  2:1746505/1‑93 (MQ=255)
 gATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAacaacg                                                                                       >  2:3281734/1‑93 (MQ=255)
  aTGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAAcaacga                                                                                      >  1:949711/1‑93 (MQ=255)
   tgctgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGAt                                                                                     <  2:4066250/93‑1 (MQ=255)
     ctgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATcc                                                                                   >  2:4643217/1‑93 (MQ=255)
      tgcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATccc                                                                                  >  1:3293545/1‑93 (MQ=255)
       gcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCa                                                                                 <  1:4270928/93‑1 (MQ=255)
       gcGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCa                                                                                 >  2:3408026/1‑93 (MQ=255)
        cGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCat                                                                                <  1:1321690/93‑1 (MQ=255)
        cGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCat                                                                                >  1:4251344/1‑93 (MQ=255)
        cGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCat                                                                                <  2:1586527/93‑1 (MQ=255)
         ggTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCatt                                                                               >  1:156977/1‑93 (MQ=255)
            cAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCattatt                                                                            <  2:3935013/93‑1 (MQ=255)
             aGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTAttt                                                                           <  2:508393/93‑1 (MQ=255)
              gCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTc                                                                          >  1:3220728/1‑93 (MQ=255)
              gCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTc                                                                          >  1:3136364/1‑93 (MQ=255)
              gCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTc                                                                          <  1:2854114/93‑1 (MQ=255)
               ccTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCa                                                                         <  1:3904933/93‑1 (MQ=255)
                cTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCAc                                                                        <  1:2467822/93‑1 (MQ=255)
                  gCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCg                                                                      <  2:1789754/93‑1 (MQ=255)
                    gATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggc                                                                    >  2:1116691/1‑93 (MQ=255)
                    gATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggc                                                                    >  1:2594778/1‑93 (MQ=255)
                     aTCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggct                                                                   <  1:2129294/93‑1 (MQ=255)
                      tCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctg                                                                  >  1:4009463/1‑93 (MQ=255)
                      tCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctg                                                                  <  2:1055373/93‑1 (MQ=255)
                        aaTAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctggc                                                                <  2:162802/93‑1 (MQ=255)
                        aaTAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctggc                                                                >  2:2231101/1‑93 (MQ=255)
                        aaTAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCggctggc                                                                <  1:218346/93‑1 (MQ=255)
                        aaTAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCagctgga                                                                >  2:2042899/1‑92 (MQ=255)
                          tAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGctc                                                              <  1:1626794/93‑1 (MQ=255)
                            gCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATGCCATTATTTCACCGGCTGGCTCtt                                                            >  2:1904494/1‑93 (MQ=255)
                            gCGGTGATACCAAAGAAGCGATTCCGGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGCCTGGCTCtt                                                            <  2:3197086/93‑1 (MQ=255)
                              ggTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGc                                                          <  2:24208/93‑1 (MQ=255)
                               gTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCa                                                         <  2:2555644/93‑1 (MQ=255)
                               gTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCa                                                         <  1:1920797/93‑1 (MQ=255)
                               gTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCa                                                         <  1:1859874/93‑1 (MQ=255)
                                ttATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCAc                                                        <  1:4613194/93‑1 (MQ=255)
                                    atCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTAt                                                    <  2:243964/93‑1 (MQ=255)
                                    atCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTAt                                                    <  1:2159616/93‑1 (MQ=255)
                                    atCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTAt                                                    <  1:201891/93‑1 (MQ=255)
                                    atCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTAt                                                    >  2:4487140/1‑93 (MQ=255)
                                     tCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGAGCTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTATTGCACCTATg                                                   <  2:1522458/93‑1 (MQ=255)
                                       aaaGAAGCTATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACGACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGcc                                                 <  1:1359356/93‑1 (MQ=255)
                                         agaagCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGt                                               >  1:1782825/1‑93 (MQ=255)
                                             gCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTaaaa                                           >  1:2227155/1‑93 (MQ=255)
                                             gCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTaaaa                                           >  1:724336/1‑93 (MQ=255)
                                              cGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTaaaaa                                          <  1:1051905/93‑1 (MQ=255)
                                              cGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTaaaaa                                          <  1:3381338/93‑1 (MQ=255)
                                               gATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAg                                         <  2:3258948/93‑1 (MQ=255)
                                               gATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAg                                         <  1:3388052/93‑1 (MQ=255)
                                                 ttCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGtt                                       >  2:1591470/1‑93 (MQ=255)
                                                 ttCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGtt                                       >  1:3519138/1‑93 (MQ=255)
                                                  tCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCGCCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTa                                      <  1:1985594/93‑1 (MQ=255)
                                                   ccAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAc                                     <  1:629870/93‑1 (MQ=255)
                                                   ccAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAc                                     <  2:1110403/93‑1 (MQ=255)
                                                    cAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCGCCTATGCCGTAAAAAGTTAcc                                    <  1:3076835/93‑1 (MQ=255)
                                                    cAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAcc                                    <  2:1287537/93‑1 (MQ=255)
                                                     aGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAccc                                   <  1:3788999/93‑1 (MQ=255)
                                                     aGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAccc                                   <  1:1342993/93‑1 (MQ=255)
                                                      gggATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCg                                  >  2:877723/1‑93 (MQ=255)
                                                      gggATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCg                                  >  1:2848224/1‑93 (MQ=255)
                                                      gggATGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCg                                  >  1:2327356/1‑93 (MQ=255)
                                                          tGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGt                              >  1:325482/1‑93 (MQ=255)
                                                          tGAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGt                              <  2:3323381/93‑1 (MQ=255)
                                                           gAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTa                             >  2:4303337/1‑93 (MQ=255)
                                                           gAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTa                             >  1:2992645/1‑93 (MQ=255)
                                                           gAAAAATCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTCCACCTATGCCGTAAAAAGTTACCCGACGTa                             <  1:3290202/93‑1 (MQ=255)
                                                            aaaaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTAt                            <  2:2266006/93‑1 (MQ=255)
                                                            aaaaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGATACCCGACGTAt                            >  1:2531477/1‑93 (MQ=255)
                                                             aaaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATc                           <  1:1005507/93‑1 (MQ=255)
                                                             aaaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATc                           <  1:634914/93‑1 (MQ=255)
                                                             aaaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATc                           <  2:299361/93‑1 (MQ=255)
                                                             aaaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATc                           <  2:2181781/93‑1 (MQ=255)
                                                             aaaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATc                           <  1:3028139/93‑1 (MQ=255)
                                                              aaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCt                          >  2:904174/1‑93 (MQ=255)
                                                              aaaTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTATTGCACCTATGCCGTAAAAAGTTACCCGACGTATCt                          <  2:1205801/93‑1 (MQ=255)
                                                                aTCTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTgg                        >  2:3355143/1‑93 (MQ=255)
                                                                  cTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCg                      <  2:2684700/93‑1 (MQ=255)
                                                                  cTGATCGCCGCACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCg                      <  2:1613825/93‑1 (MQ=255)
                                                                       cgccgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGa                 >  1:2569806/1‑93 (MQ=255)
                                                                        gccgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGaa                <  2:3449557/93‑1 (MQ=255)
                                                                        gccgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGaa                <  1:4338333/93‑1 (MQ=255)
                                                                        gccgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGaa                <  2:4294008/93‑1 (MQ=255)
                                                                         ccgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGGAc               >  2:846870/1‑93 (MQ=255)
                                                                         ccgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAAc               <  2:4492564/93‑1 (MQ=255)
                                                                         ccgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTAACCGACGTATCTGGCGGATGAAc               >  2:3886706/1‑93 (MQ=255)
                                                                          cgcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAAcc              <  1:948156/93‑1 (MQ=255)
                                                                           gcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCt             <  2:2739851/93‑1 (MQ=255)
                                                                           gcACTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCt             >  2:2852397/1‑93 (MQ=255)
                                                                              cTGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGaa          <  1:2917336/93‑1 (MQ=255)
                                                                               tGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAAt         <  1:973539/93‑1 (MQ=255)
                                                                               tGGAGGATAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAAt         <  1:1359099/93‑1 (MQ=255)
                                                                                ggaggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGCCGTATCTGGCGGATGAACCTGAATg        >  1:4528270/1‑93 (MQ=255)
                                                                                ggaggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATg        <  2:2517905/93‑1 (MQ=255)
                                                                                 gaggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATgg       <  1:4107286/93‑1 (MQ=255)
                                                                                 gaggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATgg       >  2:2859719/1‑93 (MQ=255)
                                                                                  aggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATggg      <  1:567872/93‑1 (MQ=255)
                                                                                   ggaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGc     >  2:2983031/1‑93 (MQ=255)
                                                                                    gaTAACAACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCa    >  1:1487927/1‑93 (MQ=255)
                                                                                          |                                                                                        
GGATGCTGCGGTCAGCCTGCGATCAATAGCGGTTATATCAAAGAAGCGATTCCAGGGATGAAAAATCTGATCGCCGCACTGGAGGATAACGACGATCCCATTATTTCACCGGCTGGCTCTTGCACCTATGCCGTAAAAAGTTACCCGACGTATCTGGCGGATGAACCTGAATGGGCATC  >  NC_000913/320949‑321127

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