Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 364,024 C→T 26.6% R761H (CGC→CAC)  lacZ ← beta‑D‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913364,0240CT26.6% 278.2 / 107.8 169R761H (CGC→CAC) lacZbeta‑D‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (59/65);  new base T (18/27);  total (77/92)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 170 total aligned reads displayed.

GTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACAC  >  NC_000913/363927‑364123
                                                                                                 |                                                                                                   
gTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCgg                                                                                                   <  1:3487834/100‑1 (MQ=255)
gTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCg                                                                                                    <  2:413257/99‑1 (MQ=255)
   gTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTa                                                                                                >  1:969983/1‑100 (MQ=255)
   gTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTa                                                                                                <  1:2048279/100‑1 (MQ=255)
        ccggCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTg                                                                                           <  1:2650847/97‑1 (MQ=255)
        ccAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTg                                                                                           <  1:1297373/100‑1 (MQ=255)
         cAGCGGTGCACGGGTGAACTGATCGCGCAGGGGAGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGc                                                                                          >  2:2711949/1‑100 (MQ=255)
         cAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGc                                                                                          >  1:3754806/1‑100 (MQ=255)
          aGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCTGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGcc                                                                                         <  1:2907846/100‑1 (MQ=255)
          aGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGcc                                                                                         <  1:2762063/100‑1 (MQ=255)
          aGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGcc                                                                                         <  1:460642/100‑1 (MQ=255)
               tGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATATGTGAAAGCAAGCCTGACTGGCGGTTAAATTGCCAACGc                                                                                    >  2:839784/1‑100 (MQ=255)
                   cGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCttat                                                                                <  1:3880094/100‑1 (MQ=255)
                     ggTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCttatta                                                                              >  1:1456233/1‑100 (MQ=255)
                      gTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTAc                                                                             >  2:1896902/1‑100 (MQ=255)
                      gTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTAc                                                                             <  1:1165859/100‑1 (MQ=255)
                        gAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTAccc                                                                           <  2:1793982/100‑1 (MQ=255)
                        gAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTAccc                                                                           <  1:4371/100‑1 (MQ=255)
                        gAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTAccc                                                                           <  2:600139/100‑1 (MQ=255)
                        gAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTAccc                                                                           <  1:3366750/100‑1 (MQ=255)
                        gAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTAccc                                                                           >  2:831447/1‑100 (MQ=255)
                         aaCTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCa                                                                          <  1:1563713/100‑1 (MQ=255)
                          aCTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAg                                                                         >  2:1700403/1‑100 (MQ=255)
                          aCTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAg                                                                         >  2:1264260/1‑100 (MQ=255)
                           cTGATCGCGGGGCGGCGTCAGCAGTTGTCTTTAATCGCCAATCCACATCTGTGACAACAAGCCTGACTGCCCGTGAAATTGCCAACGCCTATTAccc                                                                           >  2:909163/1‑97 (MQ=255)
                           cTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGc                                                                        >  1:1977164/1‑100 (MQ=255)
                           cTGATCGCGCAGCGGCGTCAGCAGGTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGc                                                                        <  2:1152193/100‑1 (MQ=255)
                            tGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCt                                                                       >  1:4038520/1‑100 (MQ=255)
                            tGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGATAAATTGCCAACGCTTATTACCCAGCt                                                                       >  1:2287409/1‑100 (MQ=255)
                             gATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTc                                                                      >  1:886447/1‑100 (MQ=255)
                             gATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTc                                                                      >  2:4131268/1‑100 (MQ=255)
                                cgcgCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGa                                                                    <  2:1879030/99‑1 (MQ=255)
                                 gcgcAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATg                                                                  <  1:2466104/100‑1 (MQ=255)
                                 gcgcAGCGGCGTCAGCAGTTGTTTGTTAGCGACAGTCCACATCTGTGAAAGAAAGTCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATg                                                                  <  2:1115205/100‑1 (MQ=255)
                                     aGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCaan                                                              <  1:2599177/100‑2 (MQ=255)
                                     aGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCaaa                                                              <  2:1983503/100‑1 (MQ=255)
                                      gcggcgTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCaaaa                                                             >  1:1055997/1‑100 (MQ=255)
                                        ggcgTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAAt                                                           >  2:1976438/1‑100 (MQ=255)
                                        ggcgTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAAt                                                           >  2:387676/1‑100 (MQ=255)
                                          cgTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATcc                                                         <  1:2096016/100‑1 (MQ=255)
                                            tgagcagTGGTTTTTTAACGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCAt                                                       <  2:354293/98‑1 (MQ=255)
                                            tCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCAt                                                       <  1:2278329/100‑1 (MQ=255)
                                             cagcagTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCAtt                                                      <  1:3347992/100‑1 (MQ=255)
                                             cagcaaTTGTTTTTTATCACCTATCCACATCTGTGAAAGAAAGCTTGACTGGCGGCTATATTGCCAACGCTTATTACCCAGCTTGGTGCAAAATTCTAtt                                                      >  2:1334481/1‑100 (MQ=255)
                                              agcagTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCAttt                                                     >  2:664054/1‑100 (MQ=255)
                                              agaagTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCTTGACTGGTGGTTAAATTGCCATCGCTTATTACCCAGCTCGATGCAAAAATCCAt                                                       <  2:3109355/95‑1 (MQ=255)
                                               gcagTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTc                                                    <  1:2108968/100‑1 (MQ=255)
                                                cagTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCg                                                   >  1:2477514/1‑100 (MQ=255)
                                                    tgttTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCACAAATCCATTTCGc                                                  >  2:1453840/1‑97 (MQ=255)
                                                    tgttTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCtgg                                               >  1:1330321/1‑100 (MQ=255)
                                                      ttttttATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGCTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCtggt                                              >  2:3700405/1‑99 (MQ=255)
                                                       tttttATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCtggtgg                                            >  2:2201431/1‑100 (MQ=255)
                                                       tttttATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATACATTTCGCtggt                                              >  1:277493/1‑98 (MQ=255)
                                                          ttATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCa                                         <  1:3553122/100‑1 (MQ=255)
                                                           tATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAg                                        >  1:1847136/1‑100 (MQ=255)
                                                              cGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATg                                     <  1:1929949/100‑1 (MQ=255)
                                                               gCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGc                                    >  2:3361866/1‑100 (MQ=255)
                                                               gCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGc                                    <  1:3183840/100‑1 (MQ=255)
                                                                ccAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCg                                   >  1:2428295/1‑100 (MQ=255)
                                                                  aaTCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCggg                                 <  2:2152902/100‑1 (MQ=255)
                                                                  aaTCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCggg                                 <  1:2735320/100‑1 (MQ=255)
                                                                  aaTCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTAAAAAAAAAAAAAAAAAAATCCATTTCGCTGGTGGTCAGATGc                                    <  2:3571855/97‑1 (MQ=255)
                                                                   aTCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGa                                >  1:683311/1‑100 (MQ=255)
                                                                   aTCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTCCCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGa                                <  2:2298251/100‑1 (MQ=255)
                                                                     ccACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCCAAAATCCATTTCGCTGGTGGTCAGATGCGGGATg                              <  1:2888440/100‑1 (MQ=255)
                                                                      cacaTCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATgg                             <  2:3586364/100‑1 (MQ=255)
                                                                      cacaTCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATgg                             >  1:759002/1‑100 (MQ=255)
                                                                      cacaTCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATgg                             <  2:1095379/100‑1 (MQ=255)
                                                                       acaTCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGc                            <  1:3994852/100‑1 (MQ=255)
                                                                        caTCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCg                           <  1:1417570/100‑1 (MQ=255)
                                                                          tCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTg                         <  2:3988089/100‑1 (MQ=255)
                                                                           cTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGATGGTGGTCAGATGCGGGATGGCGTgg                        >  2:1891995/1‑100 (MQ=255)
                                                                           cTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTgg                        >  2:4016583/1‑100 (MQ=255)
                                                                           cTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTgg                        <  1:4016609/100‑1 (MQ=255)
                                                                            tgtgAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTggg                       >  1:2963394/1‑100 (MQ=255)
                                                                            tgtgAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTAGTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTggg                       <  2:54158/100‑1 (MQ=255)
                                                                            tgtgAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCTAGTTCGCTGGTGGTCAGATGCGGGATGGCGTggg                       >  2:1310576/1‑100 (MQ=255)
                                                                            tgtgAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTggg                       >  1:3354488/1‑100 (MQ=255)
                                                                             ntgAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGa                      >  2:1812407/2‑100 (MQ=255)
                                                                             gtgAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACTCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGa                      <  1:926858/100‑1 (MQ=255)
                                                                             gtgAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCGTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGa                      <  2:938893/100‑1 (MQ=255)
                                                                             gtgAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCGTATTACCCAGCTCGAGGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTgg                        <  2:1677441/98‑1 (MQ=255)
                                                                               gaaagaaaGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGAcg                    <  1:364765/100‑1 (MQ=255)
                                                                               gaaagaaaGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGAcg                    >  1:2721157/1‑100 (MQ=255)
                                                                                aaagaaagCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGAcgc                   >  1:3198306/1‑100 (MQ=255)
                                                                                  agaaagCCTGACTGGCGGTTAAATTGCCAACACTTATTACCCAGCGCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGAGTGGGACgcgg                 >  2:2608102/1‑100 (MQ=255)
                                                                                   gaaagCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACgcggc                >  1:1897359/1‑100 (MQ=255)
                                                                                    aaagCCTGACTGGTGGTAAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACgcggcg               <  2:2388625/100‑1 (MQ=255)
                                                                                     aagCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGcggcgg              <  2:1210214/100‑1 (MQ=255)
                                                                                     aagCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGcggcgg              <  1:1659689/100‑1 (MQ=255)
                                                                                       gCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCgggg            >  1:715685/1‑100 (MQ=255)
                                                                                         cTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAg          >  1:3898162/1‑100 (MQ=255)
                                                                                         cTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAg          <  1:1480920/100‑1 (MQ=255)
                                                                                         cTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAg          <  2:516443/100‑1 (MQ=255)
                                                                                          tGACTGGCGGTTAAATGGCCGACGATTATTACACAGTTCTATGCAAGAATCCATTTCTGTGGTGGTCAGCTGCGGTATGGCCTGGGACacggcg               <  2:309398/94‑1 (MQ=255)
                                                                                           gACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCg        <  1:1063649/100‑1 (MQ=255)
                                                                                             cTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTc      >  1:2669181/1‑100 (MQ=255)
                                                                                              tGGCGGTTAAAGTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTca     <  2:3882587/100‑1 (MQ=255)
                                                                                                gCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTcaca   >  1:2572056/1‑100 (MQ=255)
                                                                                                 cGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCacac  <  1:3559550/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
GTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACAC  >  NC_000913/363927‑364123

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: ‑