Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,704,845 T→G K155N (AAA→AAC)  lepB ← leader peptidase (signal peptidase I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,704,8450TG90.6% 447.4 / 23.5 149K155N (AAA→AAC) lepBleader peptidase (signal peptidase I)
Reads supporting (aligned to +/- strand):  ref base T (4/10);  new base G (67/68);  total (71/78)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 152 total aligned reads displayed.

CGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACC  >  NC_000913/2704748‑2704940
                                                                                                 |                                                                                               
cGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACgtt                                                                                               <  2:851687/100‑1 (MQ=255)
cGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACgtt                                                                                               <  2:851659/100‑1 (MQ=255)
cGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTGTTGAGACCGGATCGTAAGTGACgtt                                                                                               >  1:2872748/1‑100 (MQ=255)
 gTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTa                                                                                              <  1:576863/100‑1 (MQ=255)
   ttGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATc                                                                                            <  1:2180701/100‑1 (MQ=255)
     gAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGc                                                                                          >  1:4073887/1‑100 (MQ=255)
     gAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGc                                                                                          <  2:1007877/100‑1 (MQ=255)
     gAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGACCGTAAGTGACGTTATCGc                                                                                          >  2:1793870/1‑100 (MQ=255)
       gTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGccc                                                                                        >  1:2796312/1‑100 (MQ=255)
       gTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGccc                                                                                        >  2:2943039/1‑100 (MQ=255)
         aGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCgg                                                                                      >  2:3835468/1‑100 (MQ=255)
         aGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCgg                                                                                      <  1:1506448/100‑1 (MQ=255)
          ggTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGt                                                                                     >  1:3174739/1‑100 (MQ=255)
          ggTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGt                                                                                     >  1:1319995/1‑100 (MQ=255)
          ggTGACCGGCAGCGCGTTTTCACACGCCCGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGt                                                                                     >  2:2330338/1‑100 (MQ=255)
           gTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTa                                                                                    <  1:8043/100‑1 (MQ=255)
              aCCGGCAGCGGGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAACCGCCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAAc                                                                                 <  2:2468335/100‑1 (MQ=255)
               ccGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAAcc                                                                                >  1:3378609/1‑100 (MQ=255)
               ccGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAAcc                                                                                >  2:3240503/1‑100 (MQ=255)
               ccGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAAcc                                                                                >  2:1394310/1‑100 (MQ=255)
                cGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAAccc                                                                               >  1:3006495/1‑100 (MQ=255)
                 ggCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACGCa                                                                              >  2:3613614/1‑100 (MQ=255)
                     gcgcGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACcgc                                                                          <  2:3550295/100‑1 (MQ=255)
                      cgcgTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACcgcg                                                                         <  1:1513112/100‑1 (MQ=255)
                       gcgTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACcgcgc                                                                        >  1:2549828/1‑100 (MQ=255)
                        cgTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACcgcgcg                                                                       <  2:2544436/100‑1 (MQ=255)
                        cgTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACcgcgcg                                                                       >  1:1942637/1‑100 (MQ=255)
                         ggtttCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTGAGCTCTTTTGAGGCCGGACCGTAAGTGCCGTTATCGCCCGGTAAACACACCgcgcgc                                                                      >  2:2436132/3‑100 (MQ=255)
                          ttttCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCt                                                                     >  2:3296385/1‑100 (MQ=255)
                          ttttCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCt                                                                     >  2:3373275/1‑100 (MQ=255)
                          ttttCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCt                                                                     <  1:1486524/100‑1 (MQ=255)
                            ttCACACGCCTGGCCGGAACTGCAGCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTg                                                                   <  2:551873/100‑1 (MQ=255)
                              cacaCGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGAt                                                                 <  2:1102232/100‑1 (MQ=255)
                               acacGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATg                                                                >  2:314410/1‑100 (MQ=255)
                               acacGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATg                                                                <  2:451818/100‑1 (MQ=255)
                               acacGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATg                                                                >  1:2288027/1‑100 (MQ=255)
                                 acGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTa                                                              <  1:543502/100‑1 (MQ=255)
                                     cTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCa                                                          >  1:3351451/1‑100 (MQ=255)
                                     cTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCa                                                          <  1:2409393/100‑1 (MQ=255)
                                      tGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCaa                                                         >  1:1190055/1‑100 (MQ=255)
                                       ggCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAg                                                        >  2:1615521/1‑100 (MQ=255)
                                         ccGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCt                                                      >  1:311847/1‑100 (MQ=255)
                                          cGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCtt                                                     >  1:2306302/1‑100 (MQ=255)
                                            gAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTg                                                   <  1:1599338/100‑1 (MQ=255)
                                            gAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTg                                                   <  2:1631195/100‑1 (MQ=255)
                                            gAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTg                                                   >  1:2023561/1‑100 (MQ=255)
                                            gAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTg                                                   <  1:1355939/100‑1 (MQ=255)
                                              aCTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCTCGCTTGATGTAATCAAGCTTTGGa                                                 <  2:1510563/100‑1 (MQ=255)
                                              aCTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGa                                                 <  2:630701/100‑1 (MQ=255)
                                                tGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATc                                               >  2:3561548/1‑100 (MQ=255)
                                                   aTCCCGGTTGAATCGTCAGCTCTTTTGATACCGGATCGGAAGTGACGTTATCGCCCCGTAAACCCACCGGGCGCTTGATGTAATCGAGCTTTGTAACTTc                                            >  2:2877170/1‑100 (MQ=255)
                                                   aTCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCCTGATGTAATCAAGCCTTGGCTCTTc                                            >  2:148812/1‑100 (MQ=255)
                                                    tCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTcc                                           <  1:507440/100‑1 (MQ=255)
                                                    tCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTcc                                           <  2:322880/100‑1 (MQ=255)
                                                    tCCCGGTTGAAACGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTcc                                           <  1:3639015/100‑1 (MQ=255)
                                                      ccGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCgg                                         <  1:2800409/100‑1 (MQ=255)
                                                      ccGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCgg                                         <  1:694915/100‑1 (MQ=255)
                                                      ccGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCgg                                         >  1:2059073/1‑100 (MQ=255)
                                                         gTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGata                                      >  1:1084511/1‑100 (MQ=255)
                                                         gTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGata                                      <  2:2228293/100‑1 (MQ=255)
                                                          ttGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGatat                                     <  2:2339403/100‑1 (MQ=255)
                                                            gAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATAttt                                   >  1:1200109/1‑100 (MQ=255)
                                                             aaTCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTa                                  <  1:2808476/100‑1 (MQ=255)
                                                                cGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAg                               <  1:2652437/100‑1 (MQ=255)
                                                                cGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAg                               <  2:1870014/100‑1 (MQ=255)
                                                                cGTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCCCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAg                               <  2:2919581/100‑1 (MQ=255)
                                                                 gTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGa                              >  2:2949877/1‑100 (MQ=255)
                                                                 gTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACTCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGa                              >  1:2849255/1‑100 (MQ=255)
                                                                 gTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTACTCAAGCTTTGGATCTTCCGGATATTTAAAGa                              >  1:583266/1‑100 (MQ=255)
                                                                 gTCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGa                              >  1:3790594/1‑100 (MQ=255)
                                                                  tCAGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGAc                             <  2:3056245/100‑1 (MQ=255)
                                                                    aGCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCa                           >  2:2903914/1‑100 (MQ=255)
                                                                     gCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCAc                          >  2:958292/1‑100 (MQ=255)
                                                                     gCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCAc                          >  1:2620606/1‑100 (MQ=255)
                                                                     gCTCTTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCAc                          <  2:2177958/100‑1 (MQ=255)
                                                                      ctctTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACg                         >  2:2698643/1‑100 (MQ=255)
                                                                      ctctTTTGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACg                         >  2:2850572/1‑100 (MQ=255)
                                                                         ttttGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCAcagat                      >  1:1017387/1‑96 (MQ=255)
                                                                         ttttGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGata                      <  1:1481018/100‑1 (MQ=255)
                                                                         ttttGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGata                      >  2:3700497/1‑100 (MQ=255)
                                                                           ttGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGCATATTTAAAGACCACGATATc                    <  2:2609265/100‑1 (MQ=255)
                                                                            tGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATcg                   <  1:656315/100‑1 (MQ=255)
                                                                            tGAGACCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATcg                   >  1:4053372/1‑100 (MQ=255)
                                                                              agaCCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATcgcc                 >  1:180123/1‑100 (MQ=255)
                                                                              agaCCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATcgcc                 <  1:2686904/100‑1 (MQ=255)
                                                                               gaCCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATcgccg                <  1:390280/100‑1 (MQ=255)
                                                                               gaCCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATcgccg                >  2:1917351/1‑100 (MQ=255)
                                                                               gaCCGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATcgccg                >  1:2411678/1‑100 (MQ=255)
                                                                                 ccGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCcgcg              <  2:1414963/100‑1 (MQ=255)
                                                                                 ccGGATCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCcgcg              >  1:425989/1‑100 (MQ=255)
                                                                                     aTCGTAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTc          >  1:2562079/1‑100 (MQ=255)
                                                                                      gcGTAAGTGACGTTAGCGCCCGGTAAACGCGCCGCGCGCTTGATGTAGTCACGCTGTGGATCTTCCGGATATTTAAAGACCACGATGTCGCCGCGGTTCg         <  2:1149218/99‑1 (MQ=255)
                                                                                       cGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCgg        <  2:224116/100‑1 (MQ=255)
                                                                                         tAAGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGAt      <  2:648825/100‑1 (MQ=255)
                                                                                          aaGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTGGATGTAATCAAGCTTTGGATCTTCCGGATATTTACAGACCACGATATCGCCGCGTTTCGGATg     >  1:174518/1‑100 (MQ=255)
                                                                                           aGTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGa    >  1:1579042/1‑100 (MQ=255)
                                                                                            gTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGAc   <  1:834287/100‑1 (MQ=255)
                                                                                            gTGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGAc   >  1:903577/1‑100 (MQ=255)
                                                                                             tGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGAcc  <  2:1009957/100‑1 (MQ=255)
                                                                                             tGACGTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCACGTTTCGGATGAcc  <  2:3873738/100‑1 (MQ=255)
                                                                                                 |                                                                                               
CGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACC  >  NC_000913/2704748‑2704940

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