Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,173,076 C→A 100% pseudogene (287/442 nt) gatC ← pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,0760CA94.7% 352.4 / 6.2 116pseudogene (287/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base A (45/64);  total (47/69)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 116 total aligned reads displayed.

TGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCA  >  NC_000913/2172980‑2173172
                                                                                                |                                                                                                
tGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGcc                                                                                               >  1:1081398/1‑100 (MQ=255)
tGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGcc                                                                                               <  1:2278468/100‑1 (MQ=255)
tGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGcc                                                                                               >  2:454981/1‑100 (MQ=255)
 gAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGccc                                                                                              <  1:1780351/100‑1 (MQ=255)
   aaaTACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTg                                                                                            <  2:1381358/100‑1 (MQ=255)
     aTACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGAt                                                                                          <  1:858735/100‑1 (MQ=255)
      tACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATc                                                                                         <  1:118856/100‑1 (MQ=255)
      tACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATc                                                                                         <  1:2781939/100‑1 (MQ=255)
      tACCGGTCAGATAAATTGCACAGATAATAATGAAACCGGGAATATTTTTCGGGGAGAAAAACTAAATCAATAACCAGGTAATAAGAAAACAGcccgaatc                                                                                         >  2:471852/1‑95 (MQ=255)
      tACCGGTCAGAGAAATTGCACCGATAATAATGAAACCGGGAATATGTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATc                                                                                         <  2:2269704/100‑1 (MQ=255)
       aCCGGTCAGATAAATTGCATCGATAATAATGTAACCGGGAATATTTTGTGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATcc                                                                                        <  1:918398/100‑1 (MQ=255)
       aCCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGa                                                                                           <  2:40246/97‑1 (MQ=255)
       aCCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATc                                                                                         <  2:835634/99‑1 (MQ=255)
        ccGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCa                                                                                       <  1:374187/100‑1 (MQ=255)
         cGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCAt                                                                                      <  1:965981/100‑1 (MQ=255)
          ggTCAGATAACTTGTACCGATAATAATGAAACCGGGAATAGTATGCGGTGAGAAAACCTGAATCAGTAACCAGATAATGGGAGAATAGCTCTTATCCAtt                                                                                     >  2:283308/1‑100 (MQ=255)
          ggTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCAtt                                                                                     <  2:1510010/100‑1 (MQ=255)
           gTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTg                                                                                    <  2:2829600/100‑1 (MQ=255)
           gTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTg                                                                                    <  2:411672/100‑1 (MQ=255)
            tCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGa                                                                                   <  1:2034375/100‑1 (MQ=255)
             cAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGaa                                                                                  <  1:2656204/100‑1 (MQ=255)
             cAGAGAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGCGAGAACAGCCCTGATCCAtt                                                                                     <  2:1851169/97‑1 (MQ=255)
              aGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAg                                                                                 <  2:1544505/100‑1 (MQ=255)
              aGATAAATTGAACCGATAATAATGAAACCGGGAAGATTTGGAGGGGAGAAAGCATGAATCAGTAACGAGGTAATGGGAGAACAGTCCTGATCGATTGAAg                                                                                 <  2:1159471/100‑1 (MQ=255)
               gctAAATTGCACCGATAATAATAAAACCGGGAATATTTACCGCGGGGAAAACCTGAATCAGTAACCAGGTAATTGGAGAACAGCCCTGATACAt                                                                                      >  2:1469038/3‑94 (MQ=255)
               gATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATAAGTAAAAAGGTAATGGGAGAACAGAACTGATAAATTGaa                                                                                  >  1:873407/1‑98 (MQ=255)
                   aaTTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGccacc                                                                            <  2:2723507/100‑1 (MQ=255)
                   aaTTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGccacc                                                                            <  1:2761799/100‑1 (MQ=255)
                   aaTTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGccacc                                                                            >  1:179537/1‑100 (MQ=255)
                    aTTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGccacc                                                                            <  1:1762857/99‑1 (MQ=255)
                      tGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGATAATAGAAGAACAGCCCTGATCCATTGAAGCCACCATa                                                                         >  1:102708/1‑100 (MQ=255)
                          ccGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATAcccc                                                                     >  1:1856800/1‑100 (MQ=255)
                              taataaTGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGc                                                                 >  1:2217841/1‑100 (MQ=255)
                              taataaTGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATACATTGAAGCCACCATACCCCCg                                                                   >  2:1469610/1‑98 (MQ=255)
                               aataatGTAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCATGATGCATTGAAGCCACCATAcccc                                                                     >  2:374384/1‑95 (MQ=255)
                                 taatGAAACCGGGAATATCTTGCAGGGAGAAAACATGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATGGAAGCCACCATACCCCCgcc                                                                 >  2:533094/1‑95 (MQ=255)
                                  aatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCTGTAATCAGTTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCtttt                                                             <  1:819453/100‑1 (MQ=255)
                                  aatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCtttt                                                             >  1:95832/1‑100 (MQ=255)
                                   atGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTa                                                            <  1:2712638/100‑1 (MQ=255)
                                    tgttactgtgAATTTTTTGCTTGGAGTTATCCTGAATTTGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTaa                                                           <  1:1755667/91‑1 (MQ=255)
                                    tGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACACCCGGCTTTTaa                                                           >  1:1977114/1‑100 (MQ=255)
                                     gAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAAc                                                          <  1:993395/100‑1 (MQ=255)
                                     gAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAAc                                                          >  1:634764/1‑100 (MQ=255)
                                     gAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATAGAAGCCACCATACCCCCGGCTTTTAAc                                                          >  1:659490/1‑100 (MQ=255)
                                      aaaCCGGGAATATTTTGCGGGGAGAATACCTGAATCATTAATCATGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACg                                                         <  1:502499/100‑1 (MQ=255)
                                      aaaCCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCACCCTGAACCATTGAAGCCACCTTACCCCCGGCTTTTa                                                            >  2:1032851/1‑97 (MQ=255)
                                           gggAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCa                                                    <  2:1401998/100‑1 (MQ=255)
                                           gggAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAACCACCATACCCCCGGCTTTTAACGCCCCa                                                    >  1:825218/1‑100 (MQ=255)
                                            ggAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAg                                                   <  2:273928/100‑1 (MQ=255)
                                              aaTATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCa                                                 <  2:1156380/100‑1 (MQ=255)
                                               atatTTTGCGGGGAGAAAACCTGATTCAGTAACCAGGTAATTGGAGAACAGCTCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCAt                                                <  2:662161/100‑1 (MQ=255)
                                                tatTTTGCGGTGAGAAATCCTGAATCAGTATCCACGTAATTGGAGAATAGCTCTGATCTATTGAAGCCATTATACTCTCGGTTTTTAACGCCCCAGCAt                                                <  2:2650394/99‑1 (MQ=255)
                                                tatTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCAtt                                               <  2:402509/100‑1 (MQ=255)
                                                tatTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAACCACCATACCCCCGGCTTTTAACGCCCCAGCAtt                                               <  1:2151629/100‑1 (MQ=255)
                                                    ttGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgcc                                           >  2:584812/1‑100 (MQ=255)
                                                     tGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgccg                                          >  1:1276762/1‑100 (MQ=255)
                                                      gCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgccgc                                         >  1:225142/1‑100 (MQ=255)
                                                       cGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgccgcc                                        <  1:1450345/100‑1 (MQ=255)
                                                       cGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTAAAACCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgccgcc                                        >  1:2482919/1‑100 (MQ=255)
                                                         gggAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAg                                      <  1:2937349/100‑1 (MQ=255)
                                                          ggAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGt                                     <  1:2787030/100‑1 (MQ=255)
                                                            agaAAACCTTAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTg                                   <  2:583249/100‑1 (MQ=255)
                                                             gaAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTgg                                  >  2:359150/1‑100 (MQ=255)
                                                               acaCCTGAATCAGTAACCAGGAAATGGGAGAACAGCCATGATCCATTGAAACCACCATACCACCGGCTTTTAAAGCCCCAGCATTAGCCGCCAGTTcgg                                 >  2:1019502/3‑99 (MQ=255)
                                                               aaaCCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGgt                                >  1:1938793/1‑100 (MQ=255)
                                                               aaaCCTGAATCAGTAACCAGGTAATGGGAAAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGgt                                <  1:1165582/100‑1 (MQ=255)
                                                                aaCCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGgtg                               <  1:3634/100‑1 (MQ=255)
                                                                   cTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGc                            <  2:705714/100‑1 (MQ=255)
                                                                    tGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCa                           >  2:2922264/1‑100 (MQ=255)
                                                                     gAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCaa                          <  2:16997/100‑1 (MQ=255)
                                                                     gAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCACCAGTGTGGTGTGCaa                          >  1:1422928/1‑100 (MQ=255)
                                                                      aaTCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCaaa                         >  1:1784911/1‑100 (MQ=255)
                                                                       aTCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAAc                        <  2:827599/100‑1 (MQ=255)
                                                                        gatctAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAAcc                       <  2:1466498/96‑1 (MQ=255)
                                                                          aGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCaa                     >  2:2866813/1‑100 (MQ=255)
                                                                          aGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACACCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCaa                     >  2:2422349/1‑100 (MQ=255)
                                                                           gTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTCGGGTGTGCAAACCAAt                    >  2:2892628/1‑100 (MQ=255)
                                                                             aaCCAGGTAATGGGAGAACAGCCCTGTTCTATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAg                  <  1:1016529/100‑1 (MQ=255)
                                                                             aaCCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCagatc                  >  1:1466513/1‑96 (MQ=255)
                                                                             aaCCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAg                  <  2:1794571/100‑1 (MQ=255)
                                                                             aaCCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAg                  <  1:45159/100‑1 (MQ=255)
                                                                              aCCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGt                 >  1:2588867/1‑100 (MQ=255)
                                                                              aCCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGt                 >  2:44019/1‑100 (MQ=255)
                                                                                cAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGttt               <  1:2511948/100‑1 (MQ=255)
                                                                                 aGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGGTGGGTGTGCAAACCAATAGTTTg              >  2:1703073/1‑100 (MQ=255)
                                                                                  ggTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGc             >  1:1496493/1‑100 (MQ=255)
                                                                                  ggTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGc             >  1:2182423/1‑100 (MQ=255)
                                                                                   gTAATGGTAGAACAGCCTTGATCCATTCAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGTCAGTTGGGTGTGCAAACCAATAGTTTGCg            <  2:2946796/100‑1 (MQ=255)
                                                                                   gTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCg            <  2:1969394/100‑1 (MQ=255)
                                                                                    tAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTg              <  2:95831/97‑1 (MQ=255)
                                                                                      aTGGGAGAACAGCCCTGATCCATTGAAGCCACCATTCCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACTATTGGTTTGcacac         <  1:2533645/100‑5 (MQ=255)
                                                                                      aTGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACTATTGGTTTGcacac         <  2:2863562/100‑5 (MQ=255)
                                                                                      aTGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTcg         <  2:2230150/100‑1 (MQ=255)
                                                                                       tGGGAGAACAGTTTTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTcn         <  1:1250423/99‑2 (MQ=255)
                                                                                       tGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTcgc        <  2:1938770/100‑1 (MQ=255)
                                                                                       tGGGAGAACAGCCCTGATCAATTGTAGCAACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACTATTGGTTTGcacac         <  2:2863469/99‑5 (MQ=255)
                                                                                         ggAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGa      >  1:2864150/1‑100 (MQ=255)
                                                                                            gaACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATcc   >  2:1681108/1‑100 (MQ=255)
                                                                                            gaACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATcc   <  2:699399/100‑1 (MQ=255)
                                                                                             aaCAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCca  >  2:1384572/1‑100 (MQ=255)
                                                                                                |                                                                                                
TGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCA  >  NC_000913/2172980‑2173172

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