Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,345,920 C→T intergenic (+41/‑173) ptsO → / → yrbL phosphocarrier protein HPr‑like NPr/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,345,9200CT98.8% 894.3 / ‑4.4 256intergenic (+41/‑173)ptsO/yrbLphosphocarrier protein HPr‑like NPr/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (143/110);  total (144/112)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 269 total aligned reads displayed.

CACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGA  >  NC_000913/3345830‑3346009
                                                                                          |                                                                                          
   tGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACaa                                                                                      >  1:3935207/1‑93 (MQ=255)
    ggCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAAt                                                                                     >  1:283064/1‑93 (MQ=255)
    ggCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAAt                                                                                     >  1:3810180/1‑93 (MQ=255)
     gccgccGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATa                                                                                    >  1:456207/1‑93 (MQ=255)
      gcgccgTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATaa                                                                                   <  1:1002377/92‑1 (MQ=255)
      ccgccgTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATaa                                                                                   <  1:3884915/93‑1 (MQ=255)
      ccgccgTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATaa                                                                                   <  2:2758800/93‑1 (MQ=255)
        gccgTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGa                                                                                 >  1:3590614/1‑93 (MQ=255)
        gccgTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGa                                                                                 >  2:1047681/1‑93 (MQ=255)
        gccgTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGa                                                                                 >  2:4400164/1‑93 (MQ=255)
        gccgTTATCCCCCTCTTTAATTCGGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGa                                                                                 <  2:3197008/93‑1 (MQ=255)
           gTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGaaaa                                                                              >  2:3319838/1‑93 (MQ=255)
           gTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGaaaa                                                                              >  1:1237161/1‑93 (MQ=255)
           gTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGaaaa                                                                              >  1:4133959/1‑93 (MQ=255)
           gTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGaaaa                                                                              >  1:2403764/1‑93 (MQ=255)
            ttATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAAt                                                                             >  1:3436677/1‑93 (MQ=255)
            ttATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAAt                                                                             >  2:3856400/1‑93 (MQ=255)
            ttATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAAt                                                                             >  1:595478/1‑93 (MQ=255)
             tATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATc                                                                            >  1:4397565/1‑93 (MQ=255)
              aTCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCa                                                                           >  1:3873990/1‑93 (MQ=255)
                cGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACg                                                                         <  2:817852/93‑1 (MQ=255)
                 gCCCTCTTTAATTCTGGTTTTGATGCAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGt                                                                        >  1:4719049/1‑93 (MQ=255)
                 gCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACga                                                                        >  1:243401/1‑92 (MQ=255)
                 gCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGt                                                                        >  1:4240873/1‑93 (MQ=255)
                 gCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGt                                                                        <  2:2106185/93‑1 (MQ=255)
                  cccTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCTGACAATAAGAAAATCACGt                                                                        >  2:3443795/1‑93 (MQ=255)
                   ccTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTAc                                                                      >  1:1301967/1‑93 (MQ=255)
                    ctctTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACt                                                                     >  1:2560771/1‑93 (MQ=255)
                     tctTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTg                                                                    >  1:909908/1‑93 (MQ=255)
                     tctTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTg                                                                    <  2:1768075/93‑1 (MQ=255)
                     tctTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTg                                                                    >  1:3273878/1‑93 (MQ=255)
                     tctTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTg                                                                    >  1:3644113/1‑93 (MQ=255)
                      ctTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGa                                                                   >  2:2350823/1‑93 (MQ=255)
                      ctTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGa                                                                   >  1:64663/1‑93 (MQ=255)
                       tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGaa                                                                  >  2:4728492/1‑93 (MQ=255)
                       tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCACCTGACAATAAGAAAATCACGTACTGaa                                                                  >  1:989209/1‑93 (MQ=255)
                        ttAATCCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGaaa                                                                 >  2:1654830/1‑93 (MQ=255)
                              cTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTtc                                                           <  2:2250955/93‑1 (MQ=255)
                              cTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTtc                                                           <  2:4033344/93‑1 (MQ=255)
                               tGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTtct                                                          >  1:1423405/1‑93 (MQ=255)
                                   tttGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCtcaat                                                      <  1:1553635/93‑1 (MQ=255)
                                    ttGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCtcaatc                                                     >  2:3087022/1‑93 (MQ=255)
                                       atgaAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAAc                                                  <  1:1169801/93‑1 (MQ=255)
                                       atgaAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAAc                                                  <  2:2221400/93‑1 (MQ=255)
                                        tgaAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACg                                                 >  2:1282303/1‑93 (MQ=255)
                                        tgaAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACg                                                 <  1:271023/93‑1 (MQ=255)
                                         gaAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGt                                                <  2:4171488/93‑1 (MQ=255)
                                          aaGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTc                                               >  2:2940560/1‑93 (MQ=255)
                                           aGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCa                                              <  2:4007826/93‑1 (MQ=255)
                                            gATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCACCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCAt                                             >  1:2937362/1‑93 (MQ=255)
                                             aTTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCAtt                                            <  2:3550080/93‑1 (MQ=255)
                                             aTTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCAtt                                            <  1:2777941/93‑1 (MQ=255)
                                              ttAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCAttt                                           >  2:1377624/1‑93 (MQ=255)
                                              ttAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCAttt                                           >  1:3758232/1‑93 (MQ=255)
                                               tAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTg                                          <  1:3516205/93‑1 (MQ=255)
                                                aaTCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGt                                         >  1:1860758/1‑93 (MQ=255)
                                                 aTCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTa                                        >  1:3469460/1‑93 (MQ=255)
                                                    ttCATCACTTTTACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACAt                                     >  1:2117710/1‑93 (MQ=255)
                                                    ttCATCACTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACAt                                     <  1:4372883/93‑1 (MQ=255)
                                                     tcatcaCTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACAtt                                    >  2:110540/1‑93 (MQ=255)
                                                     tcatcaCTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACAtt                                    <  2:2530341/93‑1 (MQ=255)
                                                      catcaCTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACAttt                                   >  2:4271162/1‑93 (MQ=255)
                                                      catcaCTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACAttt                                   <  2:621704/93‑1 (MQ=255)
                                                       atcaCTTTGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACAtttt                                  >  1:3090732/1‑93 (MQ=255)
                                                            tttGACATACAATTCCCTTCAAACTCCCTCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTgct                             <  1:1145147/93‑2 (MQ=255)
                                                            tttGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTgcg                             >  1:1462456/1‑93 (MQ=255)
                                                            tttGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTgcg                             >  1:741238/1‑93 (MQ=255)
                                                              tGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACTTCATTTGTACATTTTGTGCGCt                           >  1:380373/1‑93 (MQ=255)
                                                              tGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCt                           >  2:2869696/1‑93 (MQ=255)
                                                              tGACATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCt                           <  2:3168983/93‑1 (MQ=255)
                                                                aCATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCttt                         >  1:4026143/1‑93 (MQ=255)
                                                                 cATACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCtttt                        <  1:3343833/93‑1 (MQ=255)
                                                                  aTACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTc                       <  1:2611703/93‑1 (MQ=255)
                                                                  aTACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTc                       >  1:227082/1‑93 (MQ=255)
                                                                  aTACAACTACCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTc                       >  1:1746847/1‑93 (MQ=255)
                                                                   tACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCa                      >  2:2697408/1‑93 (MQ=255)
                                                                   tACAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCa                      >  2:4130815/1‑93 (MQ=255)
                                                                    aCAACTCCCTTCAACCTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCAc                     <  1:109547/93‑1 (MQ=255)
                                                                    aCAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCAc                     >  2:376922/1‑93 (MQ=255)
                                                                    aCAACTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCAc                     <  1:2635692/93‑1 (MQ=255)
                                                                      aaCTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACtt                   <  1:693653/93‑1 (MQ=255)
                                                                       aCTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACttt                  <  2:2271626/93‑1 (MQ=255)
                                                                       aCTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACttt                  <  2:1910144/93‑1 (MQ=255)
                                                                       aCTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACttt                  >  1:1707047/1‑93 (MQ=255)
                                                                       aCTCCCTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACttt                  >  2:4506476/1‑93 (MQ=255)
                                                                        cTCCCTACAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTc                 >  2:4608300/1‑93 (MQ=255)
                                                                           ccTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCaga              >  2:3739875/1‑93 (MQ=255)
                                                                            cTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCagaa             >  2:708660/1‑93 (MQ=255)
                                                                            cTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCagaa             >  2:3922216/1‑93 (MQ=255)
                                                                            cTTCAAACTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCagaa             <  1:4601425/93‑1 (MQ=255)
                                                                                 aaCTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAAcc        <  2:1450726/93‑1 (MQ=255)
                                                                                  aCTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCt       <  2:2056654/93‑1 (MQ=255)
                                                                                  aCTCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCt       <  2:1006726/93‑1 (MQ=255)
                                                                                    tCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCAc                     >  1:1992810/1‑77 (MQ=255)
                                                                                    tCCCCCTGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTa     >  2:536673/1‑93 (MQ=255)
                                                                                     ccccccGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTaa    <  1:4519441/93‑1 (MQ=255)
                                                                                      cccctGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAg   >  1:247553/1‑93 (MQ=255)
                                                                                       ccctGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGa  <  1:3024651/93‑1 (MQ=255)
                                                                                       ccctGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGa  <  2:2731311/93‑1 (MQ=255)
                                                                                       ccctGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGa  <  2:2550571/93‑1 (MQ=255)
                                                                                          |                                                                                          
CACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGA  >  NC_000913/3345830‑3346009

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