Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,410,649 A→G 19.8% I124V (ATC→GTC)  dusB → tRNA‑dihydrouridine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,410,6490AG19.8% 341.6 / 87.1 179I124V (ATC→GTC) dusBtRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base A (70/72);  new base G (18/17);  total (88/91)
Fisher's exact test for biased strand distribution p-value = 8.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

Only 100 of 179 total aligned reads displayed.

CCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAA  >  NC_000913/3410551‑3410747
                                                                                                  |                                                                                                  
cccAGATTATTGATATCAATATGGGTTGCCCGGCTAATAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGAt                                                                                                   <  2:2803056/100‑1 (MQ=255)
cccAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGAt                                                                                                   <  1:1046253/100‑1 (MQ=255)
 ccAGATTATTGATATCAATATGGGTTGCCCGGCTACAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATc                                                                                                  >  2:3077834/1‑100 (MQ=255)
 ccAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATc                                                                                                  >  2:3131393/1‑100 (MQ=255)
  cAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATcc                                                                                                 >  2:2661801/1‑100 (MQ=255)
  cAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATcc                                                                                                 <  1:136355/100‑1 (MQ=255)
      ttattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTAc                                                                                             <  1:2243901/100‑1 (MQ=255)
       tattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTAcc                                                                                            <  1:1636669/100‑1 (MQ=255)
       tattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTAcc                                                                                            <  1:3103700/100‑1 (MQ=255)
        attGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCg                                                                                           <  2:2874499/100‑1 (MQ=255)
          tGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGAAGTACCCGGATGTCGTTAAATCGATCCTTACCGAg                                                                                         >  2:411848/1‑100 (MQ=255)
           gATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAgg                                                                                        >  2:1990540/1‑100 (MQ=255)
           gATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAgg                                                                                        >  2:699067/1‑100 (MQ=255)
            atatCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGgt                                                                                       >  1:376892/1‑100 (MQ=255)
            atatCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGgt                                                                                       >  1:319246/1‑100 (MQ=255)
             tatCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGgtc                                                                                      >  1:3845843/1‑100 (MQ=255)
             tatCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGgtc                                                                                      >  2:1176618/1‑100 (MQ=255)
                cAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGgtcgtc                                                                                   >  1:3344220/1‑100 (MQ=255)
                cAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGgtcgtc                                                                                   <  2:3203766/100‑1 (MQ=255)
                 aaTATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCa                                                                                  >  1:1013539/1‑100 (MQ=255)
                    atGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATg                                                                               <  1:2519123/100‑1 (MQ=255)
                    atGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATg                                                                               <  2:3331261/100‑1 (MQ=255)
                     tgatttGCGCGGCTAAAAAAGTGAATCGCAAGCTAGCAGGTTCAGCCCTTGTGCAGTACCCGGATGTCGTGAGATCGATCCTTACCGAAGTTGCCaagg                                                                               <  2:1433811/95‑3 (MQ=255)
                     tGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGc                                                                              >  2:732050/1‑100 (MQ=255)
                      gggTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCa                                                                             >  2:3626164/1‑100 (MQ=255)
                        gTTGCCCGGCTAAGAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGt                                                                           >  1:2800909/1‑100 (MQ=255)
                        gTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGt                                                                           <  2:785400/100‑1 (MQ=255)
                        gTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGt                                                                           <  1:3648464/100‑1 (MQ=255)
                        gTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTCCAGCCCCCTTACAGTACCCGGATGTCGTTAAATCGCTCCGTACCGAGGTCGTCAATGCAGt                                                                           <  2:3905782/100‑1 (MQ=255)
                             ccGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACg                                                                      >  2:1364260/1‑100 (MQ=255)
                                gCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTc                                                                   >  1:194517/1‑100 (MQ=255)
                                 cTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTc                                                                   >  1:2533140/1‑99 (MQ=255)
                                    aaaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGt                                                               >  1:3321415/1‑100 (MQ=255)
                                     aaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGtt                                                              <  1:2955809/100‑1 (MQ=255)
                                     aaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGtt                                                              >  2:3966406/1‑100 (MQ=255)
                                      aaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTa                                                             >  2:569774/1‑100 (MQ=255)
                                       aaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTAc                                                            >  1:1443770/1‑100 (MQ=255)
                                       aaGTGAAGCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTAc                                                            <  2:3238833/100‑1 (MQ=255)
                                         gTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTAccc                                                          >  2:44448/1‑100 (MQ=255)
                                         gTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTAccc                                                          >  1:444660/1‑100 (MQ=255)
                                           gAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTg                                                        >  2:4014932/1‑100 (MQ=255)
                                             aTCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGaa                                                      >  1:2112401/1‑100 (MQ=255)
                                             aTCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGaa                                                      >  1:4007366/1‑100 (MQ=255)
                                             aTCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATAGATCCTTACCGAGGACGTCAATGCAGTGGACGTTCCTGTTACCCTGaa                                                      >  1:1378826/1‑100 (MQ=255)
                                             aTCGCAAGCTCGCAGGTTCAGCCCGCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGaa                                                      <  2:2329492/100‑1 (MQ=255)
                                              tCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAg                                                     <  2:1501215/100‑1 (MQ=255)
                                               cGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCTGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCt                                                         <  2:1905628/95‑1 (MQ=255)
                                               cGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGa                                                    >  1:49763/1‑100 (MQ=255)
                                                gCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGAt                                                   >  1:444471/1‑100 (MQ=255)
                                                 cAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGAtt                                                  <  2:2076389/100‑1 (MQ=255)
                                                 cAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGAtt                                                  >  2:864532/1‑100 (MQ=255)
                                                  aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTc                                                 <  1:1232831/100‑1 (MQ=255)
                                                  aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTc                                                 <  2:50534/100‑1 (MQ=255)
                                                     cTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCa                                              >  2:1848300/1‑100 (MQ=255)
                                                      tCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCAc                                             <  1:2714529/100‑1 (MQ=255)
                                                      tCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCAc                                             >  2:898610/1‑100 (MQ=255)
                                                       cGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCAcc                                            <  1:3680010/100‑1 (MQ=255)
                                                         cAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCgg                                          >  2:218911/1‑100 (MQ=255)
                                                         cAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCgg                                          >  2:1620853/1‑100 (MQ=255)
                                                          aGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGc                                         <  1:2331465/100‑1 (MQ=255)
                                                           ggTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCt                                        >  1:738332/1‑100 (MQ=255)
                                                            gTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTg                                       <  1:699071/100‑1 (MQ=255)
                                                            gTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTg                                       >  1:3928699/1‑100 (MQ=255)
                                                             ttCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTgg                                      <  1:3491243/100‑1 (MQ=255)
                                                              tCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTggg                                     <  1:1343125/100‑1 (MQ=255)
                                                              tCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTggg                                     >  2:2766978/1‑100 (MQ=255)
                                                               cAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGc                                    <  2:376889/100‑1 (MQ=255)
                                                               cAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGc                                    <  1:1990550/100‑1 (MQ=255)
                                                               cAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGc                                    <  2:2287085/100‑1 (MQ=255)
                                                                 gCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCAc                                  >  1:1849235/1‑100 (MQ=255)
                                                                  cccTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCAcc                                 >  2:1921627/1‑100 (MQ=255)
                                                                    ctctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCgg                               <  1:868668/100‑1 (MQ=255)
                                                                     tctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGa                              >  2:843361/1‑100 (MQ=255)
                                                                     tattGCAGTACCCGGAAGTCCTTAAATCGATCCTTACCGATTACGCCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTCGGCAc                                  >  2:4019302/3‑96 (MQ=255)
                                                                       ttGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAac                            <  1:2127464/100‑1 (MQ=255)
                                                                       ttGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAac                            <  1:512189/100‑1 (MQ=255)
                                                                        tGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAaca                           >  2:95085/1‑100 (MQ=255)
                                                                        tGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAaca                           >  2:87710/1‑100 (MQ=255)
                                                                        tGCAGTACCCGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAaca                           >  1:587632/1‑100 (MQ=255)
                                                                         gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAacac                          >  1:3998391/1‑100 (MQ=255)
                                                                         gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAacac                          >  2:1134880/1‑100 (MQ=255)
                                                                         gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAacac                          >  1:2164706/1‑100 (MQ=255)
                                                                           aGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCATCGGCTGGGCACCGGAACACCg                        <  2:1730472/100‑1 (MQ=255)
                                                                               cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAAc                    >  2:2590426/1‑100 (MQ=255)
                                                                               cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAAc                    >  2:558301/1‑100 (MQ=255)
                                                                               cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAAc                    >  1:2036582/1‑100 (MQ=255)
                                                                                 cGGATGTCGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTg                  <  1:3473551/100‑1 (MQ=255)
                                                                                 cGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACtn                  <  2:3209038/100‑2 (MQ=255)
                                                                                 cGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTg                  <  1:3737452/100‑1 (MQ=255)
                                                                                  ggATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGc                 <  2:1407671/100‑1 (MQ=255)
                                                                                   gATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCg                <  2:2342234/100‑1 (MQ=255)
                                                                                      gTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAag             <  1:360783/100‑1 (MQ=255)
                                                                                       tCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAaga            >  1:1722243/1‑100 (MQ=255)
                                                                                        cGTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAagag           >  1:2783133/1‑100 (MQ=255)
                                                                                        cGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAagag           >  1:4143172/1‑100 (MQ=255)
                                                                                         gTTAAATCGGTCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAgaga          >  1:1348742/1‑100 (MQ=255)
                                                                                            aaaTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTg       <  1:2592400/100‑1 (MQ=255)
                                                                                            aaaTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTg       <  1:712406/100‑1 (MQ=255)
                                                                                                cGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCa   >  2:1017154/1‑100 (MQ=255)
                                                                                                 gATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCaa  >  2:3590901/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
CCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAA  >  NC_000913/3410551‑3410747

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