Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,041,671 (T)6→7 intergenic (+46/‑108) dsbA → / → b3861 protein disulfide isomerase I/putative GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,041,6651.T98.6% 831.5 / 3.1 207intergenic (+40/‑114)dsbA/b3861protein disulfide isomerase I/putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base T (84/120);  total (86/121)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 258 total aligned reads displayed.

ATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGCTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGATATTCAAAGTTTTG  >  NC_000913/4041570‑4041755
                                                                                                |                                                                                            
aTGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGc                                                                                               <  1:13105316/96‑1 (MQ=255)
 tGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCt                                                                                             <  1:1637618/96‑1 (MQ=255)
  gTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtt                                                                                           <  1:3533703/96‑1 (MQ=255)
  gTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtt                                                                                           <  1:6313978/96‑1 (MQ=255)
  gTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtt                                                                                           <  2:5622952/96‑1 (MQ=255)
   tttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttt                                                                                          <  1:10920050/96‑1 (MQ=255)
   tttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttt                                                                                          <  2:7304132/96‑1 (MQ=255)
   tttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttt                                                                                          <  2:2718973/96‑1 (MQ=255)
   tatttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttt                                                                                          <  1:12875466/94‑1 (MQ=255)
     tttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttttt                                                                                       >  2:14016470/1‑96 (MQ=255)
      ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtttttt                                                                                      <  2:7155756/96‑1 (MQ=255)
      ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtttttt                                                                                      <  2:3951984/96‑1 (MQ=255)
       tGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttttttt                                                                                     <  2:357956/96‑2 (MQ=255)
       tGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttttttt                                                                                     <  1:485253/96‑2 (MQ=255)
        gTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttttttta                                                                                    >  2:10056459/1‑94 (MQ=255)
         ttCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtttttttat                                                                                   <  1:1058287/96‑4 (MQ=255)
          tCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtttttttatt                                                                                  >  2:7605924/1‑92 (MQ=255)
            agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTtattat                                                                                 >  2:8042103/1‑96 (MQ=255)
            agcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTtattat                                                                                 >  2:8352201/1‑96 (MQ=255)
             gcagTATGCTTATACAGTGAAATATCTGTCCTAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTattatt                                                                                <  2:7557124/96‑1 (MQ=255)
               agTATGCTGGTACAGTGAAATAGCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTta                                                                              >  1:1963329/1‑96 (MQ=255)
                 tATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaat                                                                            <  1:8359550/96‑1 (MQ=255)
                 tATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATGACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaat                                                                            <  2:2832912/96‑1 (MQ=255)
                  aTGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaata                                                                           <  2:8757800/96‑1 (MQ=255)
                  aTGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaata                                                                           >  2:3235747/1‑96 (MQ=255)
                    gCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTaataaa                                                                         <  1:5112156/96‑1 (MQ=255)
                    gCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTaataaa                                                                         <  2:6239603/96‑1 (MQ=255)
                    gCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTaataaa                                                                         >  1:3046463/1‑96 (MQ=255)
                      tGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTaataaata                                                                       <  1:1446309/96‑1 (MQ=255)
                      tGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTaataaata                                                                       <  1:5292535/96‑1 (MQ=255)
                       gATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAataaatat                                                                      <  1:6609182/96‑1 (MQ=255)
                        agaCAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAataaatata                                                                     >  1:7540414/3‑96 (MQ=255)
                        aTACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAataaatata                                                                     <  1:12457461/96‑1 (MQ=255)
                        aTACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAataaatata                                                                     >  1:12111215/1‑96 (MQ=255)
                         tACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAataaatataa                                                                    <  2:12523644/96‑1 (MQ=255)
                         tACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAataaatataa                                                                    <  2:6690719/96‑1 (MQ=255)
                           cAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAataaatataaat                                                                  <  1:5841204/96‑1 (MQ=255)
                           cAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAataaatataaat                                                                  <  2:11294831/96‑1 (MQ=255)
                           cAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAAGATTCACGGGCTGTTTTTATTATTTAataaatataaat                                                                  >  1:4014125/1‑96 (MQ=255)
                            aGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAAtaaatataaata                                                                 >  1:5406927/1‑96 (MQ=255)
                            aGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAAtaaatataaata                                                                 >  1:4337361/1‑96 (MQ=255)
                            aGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAAtaaatataaata                                                                 <  1:9737265/96‑1 (MQ=255)
                            aGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAAtaaatataaata                                                                 <  1:13723381/96‑1 (MQ=255)
                             gTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATAc                                                                <  2:1181423/96‑1 (MQ=255)
                             gTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATAc                                                                >  1:7386988/1‑96 (MQ=255)
                              tGAAAGATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACa                                                               <  1:2444093/96‑1 (MQ=255)
                                aaaTATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACAtt                                                             >  2:9470502/1‑96 (MQ=255)
                                  atatCTGTCCCAAAAAAAATAAATCAATGTAAATTACATAAAGCCCGTGAATATTCCCGGGCTTTTTTTATTATTTAATAAATATAAATACATTCt                                                           <  2:5473761/96‑1 (MQ=255)
                                   tatCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTg                                                          <  2:4903941/96‑1 (MQ=255)
                                    atCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGa                                                         >  2:11753874/1‑96 (MQ=255)
                                     tttGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGAt                                                        >  1:10302908/3‑96 (MQ=255)
                                     tCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCCCGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGAt                                                        <  2:6312720/96‑1 (MQ=255)
                                      cTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATa                                                       <  1:14016470/96‑1 (MQ=255)
                                      cTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATa                                                       <  1:8099027/96‑1 (MQ=255)
                                      cTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATa                                                       <  1:10911857/96‑1 (MQ=255)
                                        gTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAAt                                                     >  1:3755559/1‑96 (MQ=255)
                                         tCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATg                                                    >  1:1705779/1‑96 (MQ=255)
                                         tCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATg                                                    <  1:1228612/96‑1 (MQ=255)
                                         tCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATg                                                    <  2:10617945/96‑1 (MQ=255)
                                           cGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCTCGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCa                                                  <  2:12251916/96‑1 (MQ=255)
                                           cGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCa                                                  >  2:9176127/1‑96 (MQ=255)
                                              gaAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATcc                                               <  1:9248734/96‑1 (MQ=255)
                                               aaaaaaTAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCt                                              <  1:9991677/96‑1 (MQ=255)
                                               aaaaaaTAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCt                                              <  2:1885899/96‑1 (MQ=255)
                                                aaaaaTAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTg                                             <  2:11139173/96‑1 (MQ=255)
                                                aaaaaTAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTg                                             >  1:3354103/1‑96 (MQ=255)
                                                aaaaaTAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTg                                             <  2:9096141/96‑1 (MQ=255)
                                                  aaaTAATTCAATGTAAATTACATAAATCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGcc                                           <  2:13274111/96‑1 (MQ=255)
                                                  aaaTAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGcc                                           <  2:12192136/96‑1 (MQ=255)
                                                    ataatTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGc                                         >  2:4618107/1‑96 (MQ=255)
                                                      aatTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTg                                       <  2:14223099/96‑1 (MQ=255)
                                                      aatTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTg                                       >  1:11905076/1‑96 (MQ=255)
                                                      aatTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTg                                       <  1:295926/96‑1 (MQ=255)
                                                        ttCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGa                                     <  2:4357733/96‑1 (MQ=255)
                                                          cAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGAcc                                   <  1:1594457/96‑2 (MQ=255)
                                                          cAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACa                                   <  1:14212426/96‑1 (MQ=255)
                                                           aaTGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACat                                  <  2:13725719/96‑1 (MQ=255)
                                                            aTGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACatt                                 <  2:14312395/96‑1 (MQ=255)
                                                                 aaTTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTc                            <  1:1221252/96‑1 (MQ=255)
                                                                 aaTTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTc                            >  2:653988/1‑96 (MQ=255)
                                                                 aaTTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTc                            <  2:2415837/96‑1 (MQ=255)
                                                                  aTTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCa                           >  2:8510233/1‑96 (MQ=255)
                                                                    tACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGc                         <  2:1284566/96‑1 (MQ=255)
                                                                       aTAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATa                      <  1:6807601/96‑1 (MQ=255)
                                                                        tAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCAtta                     <  1:10811234/96‑3 (MQ=255)
                                                                        tAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATaa                     >  2:11563176/1‑96 (MQ=255)
                                                                          aaGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGt                   <  2:6141730/96‑1 (MQ=255)
                                                                             cccGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGa                >  2:1083111/1‑96 (MQ=255)
                                                                             cccGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGa                <  2:2983796/96‑1 (MQ=255)
                                                                              ccGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGat               >  1:9414898/1‑96 (MQ=255)
                                                                              ccGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGat               >  2:2309139/1‑96 (MQ=255)
                                                                               cGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGata              <  1:2232905/96‑1 (MQ=255)
                                                                                gTGAATATTCACGGGCTTTGTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGatat             <  1:11602107/96‑1 (MQ=255)
                                                                                 tGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGatagt            >  1:10409035/1‑94 (MQ=255)
                                                                                 tGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGATAtt            >  1:9652785/1‑96 (MQ=255)
                                                                                    atATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGATATTCaa         >  2:3888615/2‑96 (MQ=255)
                                                                                          aCGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGATATTCAAAGtttt   >  1:6976017/1‑96 (MQ=255)
                                                                                          aCGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGATATTCAAAGtttt   <  2:6911204/96‑1 (MQ=255)
                                                                                           cGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGATATTCAAAGTTTTg  >  1:4551353/1‑96 (MQ=255)
                                                                                           cGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGATATTCAAAGTTTTg  <  1:13818957/96‑1 (MQ=255)
                                                                                                |                                                                                            
ATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGCTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGATATTCAAAGTTTTG  >  NC_000913/4041570‑4041755

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