Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,224,242 A→T 15.5% intergenic (+335/+37) ycgI → / ← minE pseudogene/cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,224,2420AT15.5% 326.3 / 50.3 129intergenic (+335/+37)ycgI/minEpseudogene/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base A (46/63);  new base T (9/11);  total (55/74)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-02

Only 100 of 129 total aligned reads displayed.

GGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCA  >  NC_000913/1224146‑1224339
                                                                                                |                                                                                                 
gggggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAActc                                                                                                >  2:1255427/1‑100 (MQ=255)
   ggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGc                                                                                             <  1:1617812/100‑1 (MQ=255)
   ggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGc                                                                                             <  2:1105438/100‑1 (MQ=255)
      ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCtt                                                                                          <  2:664450/100‑1 (MQ=255)
      ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCtt                                                                                          >  2:113347/1‑100 (MQ=255)
       gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTg                                                                                         <  2:819094/100‑1 (MQ=255)
       gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTg                                                                                         <  1:1489963/100‑1 (MQ=255)
       gTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTg                                                                                         <  2:952892/100‑1 (MQ=255)
        tGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGa                                                                                        <  1:2842538/100‑1 (MQ=255)
        tGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGa                                                                                        >  1:2035001/1‑100 (MQ=255)
        tGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGa                                                                                        <  2:736640/100‑1 (MQ=255)
         ggAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGaa                                                                                       >  1:1499600/1‑100 (MQ=255)
         ggAAATAGCAATGAGGAGTATCAGCAAGAATACGCGCCGCTTTACCACAACGTGGACGAGAGCGATGAAAACCTCAAGGCAGAGTTATCTCTGCCt                                                                                           <  1:1371782/96‑1 (MQ=255)
          gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGagt                                                                                      <  1:1184369/100‑3 (MQ=255)
          gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGagt                                                                                      <  2:2587215/100‑3 (MQ=255)
          gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGagt                                                                                      <  2:2587279/100‑3 (MQ=255)
          gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAg                                                                                      <  2:491117/100‑1 (MQ=255)
            aaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGagttt                                                                                    <  1:658868/100‑5 (MQ=38)
             aTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGagtttt                                                                                   <  2:1609760/100‑6 (MQ=255)
             aTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATa                                                                                   >  2:349965/1‑100 (MQ=255)
             aTAGAAATGACGAGTCCCAGCAAGAATAATCGCCGTATTACCACAACGGGAATAAAAAAAATAAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATa                                                                                   <  1:879763/100‑1 (MQ=255)
                gCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCTTGagtttttca                                                                                <  2:2284558/100‑9 (MQ=37)
                gCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                                >  2:685716/1‑100 (MQ=255)
                  aaTGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgc                                                                              >  1:137093/1‑100 (MQ=255)
                   aTGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcg                                                                             >  2:2529026/1‑100 (MQ=255)
                   aTGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgc                                                                              >  1:1300225/1‑99 (MQ=255)
                    tGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcgc                                                                            <  1:349916/100‑1 (MQ=255)
                     gaggagTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCt                                                                           <  2:2781825/100‑1 (MQ=255)
                         agTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGTGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                                                       >  2:2620694/1‑100 (MQ=255)
                          gTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTTAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTAc                                                                      <  2:1279538/100‑1 (MQ=255)
                             tCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGc                                                                   <  2:1218277/100‑1 (MQ=255)
                             tCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGc                                                                   >  1:1044013/1‑100 (MQ=255)
                               aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                                                 <  2:2275727/100‑1 (MQ=255)
                               aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCgg                                                                 >  2:367153/1‑100 (MQ=255)
                               aGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCg                                                                  >  1:402883/1‑99 (MQ=255)
                                gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCggg                                                                >  2:1173608/1‑100 (MQ=255)
                                gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCggg                                                                <  1:113346/100‑1 (MQ=255)
                                   aGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCtt                                                             >  1:2526367/1‑100 (MQ=255)
                                        aCTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTc                                                        <  2:2831525/100‑1 (MQ=255)
                                        aCTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTc                                                        <  2:2084260/100‑1 (MQ=255)
                                        aCTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTc                                                        >  2:2245919/1‑100 (MQ=255)
                                         cTCGCCGCTTTACCACAACGTTGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                       <  1:2102140/100‑1 (MQ=255)
                                         cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                       <  2:165661/100‑1 (MQ=255)
                                         cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                       <  2:2483366/100‑1 (MQ=255)
                                         cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                       >  2:2772479/1‑100 (MQ=255)
                                          tCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAg                                                      <  2:1986235/100‑1 (MQ=255)
                                           cgccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                                                     <  2:1823934/100‑1 (MQ=255)
                                           cgccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                                                     >  2:1740834/1‑100 (MQ=255)
                                           cgccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                                                     >  2:1950212/1‑100 (MQ=255)
                                           cgccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                                                     <  2:1634785/100‑1 (MQ=255)
                                           cgccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                                                     <  1:1116089/100‑1 (MQ=255)
                                           cgccgcTTTACCACAACGTGGATGAGAGGCATGAAAAACTCAACGCAGAGATAACTCTGCCTTGCAGATAAATACACTTTTACAACGGGCTTATTTCAGc                                                     >  1:837875/1‑100 (MQ=255)
                                             ccgcTTTACCACAACGTCGATGAAAGGGATGAAAAACTCAAGGCAGAGACAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCa                                                       <  1:139122/96‑1 (MQ=255)
                                               gcTTTACTAAAACGTGGATGAGATGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGc                                                     <  1:1939985/96‑1 (MQ=255)
                                               gcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtctt                                                 >  2:926/1‑100 (MQ=255)
                                                cTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc                                                <  2:78946/100‑1 (MQ=255)
                                                cTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc                                                <  1:2888353/100‑1 (MQ=255)
                                                     ccACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCtt                                           >  1:2477944/1‑100 (MQ=255)
                                                     ccACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCtt                                           >  1:2814436/1‑100 (MQ=255)
                                                     ccACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCtt                                           >  2:2144900/1‑100 (MQ=255)
                                                       acaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTcc                                         <  1:1281191/100‑1 (MQ=255)
                                                         aaCGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCgg                                       <  1:26463/100‑1 (MQ=255)
                                                             tGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAg                                   <  2:2169561/100‑1 (MQ=255)
                                                              ggATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAgg                                  <  1:2144887/100‑1 (MQ=255)
                                                              ggATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAgg                                  <  2:2249267/100‑1 (MQ=255)
                                                                  gagaGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCAc                              <  2:717952/100‑1 (MQ=255)
                                                                    gagGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGt                            <  2:1985489/100‑1 (MQ=255)
                                                                     agGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGtt                           <  2:2912366/100‑1 (MQ=255)
                                                                      gggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTc                          <  1:633076/100‑1 (MQ=255)
                                                                      gggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTc                          >  1:2095503/1‑100 (MQ=255)
                                                                      gggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTc                          >  2:2354800/1‑100 (MQ=255)
                                                                       ggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCa                         <  2:256299/100‑1 (MQ=255)
                                                                        gATGAAAAACTCAAGGAAGAGATAACTCTGCCTTGAAGATAAAGGCGCTTTTACAGCGGGTTTATTTCAGCTGTTCTGTTTCCGGTAAGGTCACGTTCa                         <  1:1632566/99‑1 (MQ=255)
                                                                         aTGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAgn                       <  2:1430128/100‑2 (MQ=255)
                                                                           gAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTc                     <  1:2184739/100‑1 (MQ=255)
                                                                            aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCa                    <  2:2231073/100‑1 (MQ=255)
                                                                             aaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCaa                   <  2:1294562/100‑1 (MQ=255)
                                                                              aaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAg                  <  2:2642490/100‑1 (MQ=255)
                                                                               tcttCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGa                 >  1:86140/4‑100 (MQ=255)
                                                                                cttCAAGGCAGCGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGaa                >  1:1530421/3‑100 (MQ=255)
                                                                                cttCAAGGCAGAGTTATCTCTGCCTTGAGGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGaa                >  1:1460983/3‑100 (MQ=38)
                                                                                 ttCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAAt               >  1:682521/2‑100 (MQ=255)
                                                                                 ttCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAAt               >  2:1822469/2‑100 (MQ=255)
                                                                                  tCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATa              >  2:2327991/1‑100 (MQ=255)
                                                                                  tCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCGTTTACAGCGGGCTTATTTCAGCTATTCTGGTTCCGGTAAGGTCACGTTCAGCTCAAGAAt               >  1:2178473/1‑99 (MQ=255)
                                                                                   cAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAg             >  2:1187930/1‑100 (MQ=255)
                                                                                   cAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAg             >  2:57535/1‑100 (MQ=255)
                                                                                    aaGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa            >  1:1272185/1‑100 (MQ=255)
                                                                                    aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa            <  2:1271249/100‑1 (MQ=255)
                                                                                    aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGa            >  1:1813594/1‑100 (MQ=255)
                                                                                     aGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaa           <  2:50588/100‑1 (MQ=255)
                                                                                      ggCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGaaa          <  2:75556/100‑1 (MQ=255)
                                                                                        cAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAata        >  2:930855/1‑100 (MQ=255)
                                                                                           agaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCg     <  1:1863319/100‑1 (MQ=255)
                                                                                            gaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGTGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGc    <  2:2730730/100‑1 (MQ=255)
                                                                                            gaTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGc    <  1:878822/100‑1 (MQ=255)
                                                                                             aTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGcc   <  1:1950203/100‑1 (MQ=255)
                                                                                             aTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGcc   <  2:161397/100‑1 (MQ=255)
                                                                                             aTAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGc    <  2:1955847/99‑1 (MQ=255)
                                                                                              tAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCa  >  1:2111475/1‑100 (MQ=255)
                                                                                                |                                                                                                 
GGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCA  >  NC_000913/1224146‑1224339

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: