Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,529,947 C→T A36A (GCC→GCT)  ydcC → putative receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,529,9470CT99.4% 516.4 / ‑5.5 157A36A (GCC→GCT) ydcCputative receptor
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (75/81);  total (76/81)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 157 total aligned reads displayed.

GATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGA  >  NC_000913/1529857‑1530006
                                                                                          |                                                           
gATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGt                                                           >  1:1110400/1‑93 (MQ=255)
gATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGt                                                           >  1:2866397/1‑93 (MQ=255)
 aTGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGtt                                                          <  1:1606101/93‑1 (MQ=255)
  tGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTa                                                         >  1:546309/1‑93 (MQ=255)
  tGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTa                                                         <  1:2298076/93‑1 (MQ=255)
    ggACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTAtt                                                       <  2:646803/93‑1 (MQ=255)
    ggACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTAtt                                                       <  1:2553171/93‑1 (MQ=255)
    ggACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTAtt                                                       >  2:4143594/1‑93 (MQ=255)
    ggACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTAtt                                                       >  1:1838676/1‑93 (MQ=255)
      aCATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTc                                                     <  1:2962906/93‑1 (MQ=255)
      aCATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTc                                                     <  1:391987/93‑1 (MQ=255)
       cATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTAtttat                                                    <  2:1932633/93‑3 (MQ=255)
       cATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCt                                                    >  1:3216172/1‑93 (MQ=255)
        atatTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTg                                                   <  2:1992663/93‑1 (MQ=255)
        atatTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTg                                                   <  2:2015400/93‑1 (MQ=255)
        atatTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTg                                                   >  1:1588947/1‑93 (MQ=38)
         tatTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTgg                                                  >  2:1923463/1‑93 (MQ=255)
         tatTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTgg                                                  <  1:3604722/93‑1 (MQ=255)
          atTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGt                                                 >  1:4205371/1‑93 (MQ=37)
          atTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGt                                                 >  1:1364728/1‑93 (MQ=37)
          atTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGt                                                 <  2:1288011/93‑1 (MQ=255)
           tttCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTg                                                <  2:3449667/93‑1 (MQ=255)
            ttCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGc                                               <  1:3318675/93‑1 (MQ=255)
               tattatCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGa                                            <  2:3092071/93‑1 (MQ=255)
               tattatCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGa                                            >  2:1268024/1‑93 (MQ=37)
               tattatCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGa                                            >  1:1766911/1‑93 (MQ=255)
                attatCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGaa                                           >  1:1078528/1‑93 (MQ=255)
                attatCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGaa                                           <  1:3432826/93‑1 (MQ=255)
                attatCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGaa                                           <  1:3432825/93‑1 (MQ=255)
                attatCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGaa                                           <  1:2226573/93‑1 (MQ=255)
                attatCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGaa                                           <  2:3047657/93‑1 (MQ=255)
                 ttatCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAg                                          >  2:2451465/1‑93 (MQ=37)
                 ttatCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAg                                          <  1:4562800/93‑1 (MQ=255)
                  tatCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAgg                                         >  2:3870480/1‑93 (MQ=255)
                  tatCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAgg                                         <  1:3224392/93‑1 (MQ=255)
                   atCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGc                                        <  2:4305049/93‑1 (MQ=255)
                    tCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCt                                       >  1:4184164/1‑93 (MQ=37)
                    tCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCt                                       <  1:1434880/93‑1 (MQ=37)
                    tCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCt                                       >  2:1330889/1‑93 (MQ=255)
                    tCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCt                                       >  1:2037979/1‑93 (MQ=37)
                     ccccGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTg                                      <  2:3660678/93‑1 (MQ=37)
                      cccGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTgg                                     <  2:1842019/93‑1 (MQ=255)
                      cccGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTgg                                     <  2:1627796/93‑1 (MQ=255)
                      cccGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTgg                                     <  2:2934201/93‑1 (MQ=255)
                      cccGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTgg                                     <  1:4478601/93‑1 (MQ=255)
                        cGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGa                                   <  2:1541513/93‑1 (MQ=255)
                        cGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGa                                   <  2:790867/93‑1 (MQ=255)
                         gATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGaa                                  >  2:4366639/1‑93 (MQ=37)
                         gATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGaa                                  >  2:565192/1‑93 (MQ=255)
                         gATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGaa                                  >  2:526092/1‑93 (MQ=255)
                         gATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGaa                                  <  1:4036194/93‑1 (MQ=255)
                          aTTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAg                                 >  1:4648384/1‑93 (MQ=255)
                           ttACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGa                                <  1:1947169/93‑1 (MQ=37)
                           ttACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGa                                >  2:1342452/1‑93 (MQ=255)
                            tACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGat                               <  1:1648637/93‑1 (MQ=255)
                            tACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGat                               <  2:573245/93‑1 (MQ=255)
                             aCAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGata                              <  1:2577593/93‑1 (MQ=255)
                              cAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGatat                             <  1:4717208/93‑1 (MQ=255)
                               aGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGAtata                            >  2:4364418/1‑93 (MQ=255)
                               aGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGAtata                            >  2:494582/1‑93 (MQ=37)
                                gACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATag                           >  2:4041338/1‑93 (MQ=37)
                                gACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATag                           <  2:1145075/93‑1 (MQ=37)
                                 aCAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATaga                          <  1:3299499/93‑1 (MQ=37)
                                 aCAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATaga                          <  2:2161093/93‑1 (MQ=255)
                                 aCAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATaga                          <  2:3349752/93‑1 (MQ=255)
                                 aCAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATaga                          <  2:1070636/93‑1 (MQ=255)
                                  cAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATagag                         >  2:3419390/1‑93 (MQ=255)
                                  cAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATagag                         >  1:3948075/1‑93 (MQ=255)
                                   aaGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAgg                        >  2:2405550/1‑93 (MQ=255)
                                   aaGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAgg                        >  1:1604284/1‑93 (MQ=255)
                                   aaGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAgg                        <  1:4621869/93‑1 (MQ=255)
                                    aGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGa                       <  1:407631/93‑1 (MQ=37)
                                    aGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGa                       >  2:1587363/1‑93 (MQ=255)
                                    aGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGa                       <  2:2171920/93‑1 (MQ=37)
                                     gCCTGGAAAATGGAACATAAGTTNTCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGAt                      <  1:1938142/93‑1 (MQ=255)
                                     gCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGAt                      <  2:657818/93‑1 (MQ=255)
                                      ccTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGAtt                     >  2:2567562/1‑93 (MQ=255)
                                      ccTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGAtt                     >  2:1778223/1‑93 (MQ=255)
                                      ccTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGAtt                     >  1:73229/1‑93 (MQ=255)
                                       cTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGAttt                    <  2:3020290/93‑1 (MQ=255)
                                       cTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGAttt                    >  2:113115/1‑93 (MQ=255)
                                       cTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGAttt                    >  1:1715054/1‑93 (MQ=255)
                                        tGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGAtttt                   <  2:1716485/93‑1 (MQ=37)
                                        tGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGAtttt                   <  2:4580275/93‑1 (MQ=37)
                                          gAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTgg                 >  2:2624023/1‑93 (MQ=255)
                                          gAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTgg                 <  2:2574801/93‑1 (MQ=255)
                                            aaaTGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTATTTTCTGGTTCAGAGGGTTGGGAGGATATAGAGGATTTTgggg               >  1:4575314/1‑93 (MQ=11)
                                             aaTGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGa              <  1:4204041/93‑1 (MQ=25)
                                               tGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGCTTTTGGGGaaa            >  2:1000921/1‑93 (MQ=16)
                                               tGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGaac            >  1:1526748/1‑92 (MQ=19)
                                               tGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGaaa            <  1:2965319/93‑1 (MQ=18)
                                                 gAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGGAGGCTGGGAAGATATAGAGGATTTTGGGGAcaca          >  1:1058109/1‑93 (MQ=11)
                                                 gAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAaca          >  1:2244147/1‑93 (MQ=21)
                                                 gAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAaca          >  1:2379162/1‑93 (MQ=21)
                                                 gAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAaca          >  2:450568/1‑93 (MQ=21)
                                                       aaGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGTTGCAGAAGGCTGGGAAGACATAGAGGATTTTGGGGAAACACATccc    <  2:449670/93‑1 (MQ=255)
                                                       aaGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATccc    <  2:4116071/93‑1 (MQ=18)
                                                       aaGTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATccc    <  2:3659808/93‑1 (MQ=18)
                                                         gTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGa  <  1:3166114/93‑1 (MQ=14)
                                                         gTTATCGGATATTCTACTGTTGACTATTTGTGCTGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAAAACATCCCGa  >  1:3775324/1‑93 (MQ=9)
                                                                                          |                                                           
GATGGGACATATTTCTATTATCCCCGATTACAGACAAGCCTGGAAAATGGAACATAAGTTATCGGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGCTGGGAAGATATAGAGGATTTTGGGGAAACACATCCCGA  >  NC_000913/1529857‑1530006

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