Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,672,562 T→C L209P (CTG→CCG)  b1601 → putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,672,5620TC99.5% 619.9 / ‑5.4 197L209P (CTG→CCG) b1601putative transport protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (89/107);  total (89/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 209 total aligned reads displayed.

CGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATCCCGAATA  >  NC_000913/1672472‑1672651
                                                                                          |                                                                                         
 ggTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGt                                                                                        <  2:963371/93‑1 (MQ=255)
 ggTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGt                                                                                        <  2:6334784/93‑1 (MQ=255)
  gTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTc                                                                                       >  1:439588/1‑93 (MQ=255)
   ttGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCg                                                                                      >  2:3302503/1‑93 (MQ=255)
   ttGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCg                                                                                      >  1:5772139/1‑93 (MQ=255)
      aaGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCt                                                                                   <  2:3127595/93‑1 (MQ=255)
      aaGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCt                                                                                   <  1:3514322/93‑1 (MQ=255)
      aaGGGAGGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCt                                                                                   <  2:1409358/93‑1 (MQ=255)
       aGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTg                                                                                  <  2:984246/93‑1 (MQ=255)
       aGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTg                                                                                  <  2:3897483/93‑1 (MQ=255)
       aGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTg                                                                                  <  2:3833373/93‑1 (MQ=255)
       aGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTg                                                                                  >  1:5689777/1‑93 (MQ=255)
        gggATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTgg                                                                                 >  2:4600462/1‑93 (MQ=255)
           aTGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCg                                                                              >  2:2163410/1‑93 (MQ=255)
            tGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGa                                                                             <  2:6307915/93‑1 (MQ=255)
            tGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGa                                                                             <  1:4369673/93‑1 (MQ=255)
            tGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGa                                                                             <  2:6985767/93‑1 (MQ=255)
            tGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGa                                                                             >  2:8626005/1‑93 (MQ=255)
            tGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGAGCGa                                                                             >  2:9239573/1‑93 (MQ=255)
             ggcggcGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGAt                                                                            <  2:4066100/93‑1 (MQ=255)
             ggcggcGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGAt                                                                            >  1:572377/1‑93 (MQ=255)
              gcggggATTCCAAGTGCGATTGGCCGTTTTTCCTCTTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATg                                                                           <  2:8712215/93‑1 (MQ=255)
              gcggcgATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATg                                                                           >  1:4768705/1‑93 (MQ=255)
                  cgATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTcg                                                                       <  2:2389200/93‑1 (MQ=255)
                    aTTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTcgcc                                                                     >  1:9115288/1‑93 (MQ=255)
                     ttCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTcgccc                                                                    >  2:7976307/1‑92 (MQ=255)
                       cAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCa                                                                  >  2:10015327/1‑93 (MQ=255)
                        aaCGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCAc                                                                 >  1:8530914/1‑93 (MQ=255)
                         aCGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACt                                                                <  2:5675004/93‑1 (MQ=255)
                           gTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCg                                                              >  2:4199960/1‑93 (MQ=255)
                                aTTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGtt                                                         >  1:263039/1‑93 (MQ=255)
                                 ttGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttc                                                        >  1:4555209/1‑93 (MQ=255)
                                 ttGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttc                                                        >  2:6900699/1‑93 (MQ=255)
                                  tGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcg                                                       <  2:5469906/93‑1 (MQ=255)
                                   gACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgc                                                      <  1:4983962/93‑1 (MQ=255)
                                   gACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgc                                                      <  2:7592776/93‑1 (MQ=255)
                                    aCAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgct                                                     >  1:7231425/1‑93 (MQ=255)
                                    aCAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgct                                                     >  2:251774/1‑93 (MQ=255)
                                    aCAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgct                                                     <  2:7854225/93‑1 (MQ=255)
                                     cAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgctt                                                    <  2:6476235/93‑1 (MQ=255)
                                      aGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgcttc                                                   >  1:7023544/1‑93 (MQ=255)
                                       gtgtTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgcttcg                                                  >  2:3690207/1‑93 (MQ=255)
                                        tgtTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGttcgcttcgc                                                 <  1:5336027/93‑1 (MQ=255)
                                          tttCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTcgcttcgctt                                               >  1:1796310/1‑93 (MQ=255)
                                            tctcATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCtttt                                             <  1:2977899/93‑1 (MQ=255)
                                            tctcATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCtttt                                             <  2:1179326/93‑1 (MQ=255)
                                            tctcATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCtttt                                             <  2:1522248/93‑1 (MQ=255)
                                            tctcATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCtttt                                             >  2:9744007/1‑93 (MQ=255)
                                             ctcATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTg                                            <  1:1651067/93‑1 (MQ=255)
                                              tcATTATCTTGGGCTGAAAACAACCATCAGAATCATCACCCGCCCGGTCGGCTTGGGGGTGCTCGCCGCACTCGATGGTTGCTTCGCTTTTGt                                           >  2:7090331/1‑93 (MQ=255)
                                                aTTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCt                                         <  1:1525767/93‑1 (MQ=255)
                                                aTTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCt                                         <  2:1510007/93‑1 (MQ=255)
                                                 ttATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTg                                        <  1:4089962/93‑1 (MQ=255)
                                                  tATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTgg                                       <  1:3809750/93‑1 (MQ=255)
                                                  tATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTgg                                       <  1:686845/93‑1 (MQ=255)
                                                    tCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTgggg                                     >  2:7201186/1‑93 (MQ=255)
                                                    tCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTgggg                                     >  1:8570918/1‑93 (MQ=255)
                                                     cTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGa                                    <  2:2617949/93‑1 (MQ=255)
                                                     cTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGa                                    <  1:3950188/93‑1 (MQ=255)
                                                     cTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGa                                    <  1:7516116/93‑1 (MQ=255)
                                                      tGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGAt                                   <  2:2746617/93‑1 (MQ=255)
                                                      tGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGAt                                   >  2:6881857/1‑93 (MQ=255)
                                                         tGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGc                                <  1:4753407/93‑1 (MQ=255)
                                                         tGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGc                                >  1:5253833/1‑93 (MQ=255)
                                                         tGCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGc                                >  1:1616296/1‑93 (MQ=255)
                                                          gCTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCt                               >  1:8692964/1‑93 (MQ=255)
                                                           cTGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTg                              <  2:19446/93‑1 (MQ=255)
                                                            tGAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTgg                             <  1:6067658/93‑1 (MQ=255)
                                                             gAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGc                            >  1:3998677/1‑93 (MQ=255)
                                                             gAAAACAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGc                            >  2:1918360/1‑93 (MQ=255)
                                                               aaaCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCt                          <  1:6436697/93‑1 (MQ=255)
                                                                aaCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCtt                         <  2:5889429/93‑1 (MQ=255)
                                                                aaCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCtt                         >  1:7591796/1‑93 (MQ=255)
                                                                 aCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCttt                        >  1:2430560/1‑93 (MQ=255)
                                                                 aCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCttt                        <  2:1600432/93‑1 (MQ=255)
                                                                 aCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCttt                        <  2:9076248/93‑1 (MQ=255)
                                                                 aCAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCttt                        <  2:6090349/93‑1 (MQ=255)
                                                                  cAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCCGGCCTTTg                       <  1:9335871/93‑1 (MQ=255)
                                                                  cAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTg                       >  2:7426159/1‑93 (MQ=255)
                                                                  cAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTg                       <  1:1591179/93‑1 (MQ=255)
                                                                  cAGCCATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTg                       >  1:2395034/1‑93 (MQ=255)
                                                                      cATCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCt                   >  2:2237992/1‑93 (MQ=255)
                                                                       aTCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCtt                  >  2:8923512/1‑93 (MQ=255)
                                                                        tCAGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTa                 >  1:3435520/1‑93 (MQ=255)
                                                                          cgCATCATCCGCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAAt               <  1:1932683/92‑1 (MQ=255)
                                                                          aGCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAAt               <  2:1236047/93‑1 (MQ=255)
                                                                           gCATCATCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAAtt              >  1:4869922/1‑93 (MQ=255)
                                                                            catcatCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTa             >  2:3612508/1‑93 (MQ=255)
                                                                            catcatCACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTa             >  1:9475726/1‑93 (MQ=255)
                                                                             atcatcACCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTAc            >  1:7624668/1‑93 (MQ=255)
                                                                              tcatcaCCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACa           <  2:3174735/93‑1 (MQ=255)
                                                                                atcgCCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATc         <  1:8320111/93‑1 (MQ=255)
                                                                                 tcaCCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATcc        <  2:3992231/93‑1 (MQ=255)
                                                                                 acaCCGGCCCGGTCGCCTGAGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTAGGGATTGCTGGCCTTTGCGCTTAATTACATca        <  1:8969700/92‑2 (MQ=255)
                                                                                  caCCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATccc       <  2:8937788/93‑1 (MQ=255)
                                                                                  caCCGGCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATccc       <  2:4495675/93‑1 (MQ=255)
                                                                                      ggCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATCCCGAat   >  1:9136687/1‑93 (MQ=255)
                                                                                      ggCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATCCCGAat   >  1:8158903/1‑93 (MQ=255)
                                                                                       gCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATCCCGAata  <  1:1171225/93‑1 (MQ=255)
                                                                                       gCCCGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATCCCGAata  >  2:5034231/1‑93 (MQ=255)
                                                                                          |                                                                                         
CGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTGTTTCTCATTATCTGGTGCTGAAAACAGCCATCAGCATCATCACCGGCCTGGTCGCCTGGGCGATGCTCGCCGCACTCGATGTTCGCTTCGCTTTTGTCTGGGGATTGCTGGCCTTTGCGCTTAATTACATCCCGAATA  >  NC_000913/1672472‑1672651

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