Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,108,133 A→G G208G (GGA→GGG cdh → CDP‑diacylglycerol phosphotidylhydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,108,1330AG99.7% 1051.8 / NA 336G208G (GGA→GGGcdhCDP‑diacylglycerol phosphotidylhydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (128/207);  minor base T (0/1);  total (128/208)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 350 total aligned reads displayed.

GGCGCGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGACGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTC  >  NC_000913/4108043‑4108222
                                                                                          |                                                                                         
ggCGCGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCg                                                                                         >  2:5568177/1‑93 (MQ=255)
 gcgcGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGgcgc                                                                                        <  2:119118/93‑1 (MQ=255)
 gcgcGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGgcgc                                                                                        <  1:6741483/93‑1 (MQ=255)
    cgTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTAc                                                                                     >  2:441828/1‑93 (MQ=255)
    cgTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTAc                                                                                     <  1:3846243/93‑1 (MQ=255)
     gTCGGGTATCGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACg                                                                                    <  2:6188941/93‑1 (MQ=255)
     gTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACg                                                                                    >  1:6177783/1‑93 (MQ=255)
     gTCGGGTAACGGAAAGCGAACTGGTACAACGCAGACCATTTCTGATGCTAGCAGAAGAAGTACCTAAGGCGAGGGAACATATGGGGCGCTACg                                                                                    >  1:2811430/1‑93 (MQ=255)
      tCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACgg                                                                                   >  2:2161339/1‑93 (MQ=255)
        gggTAACGGAAAGCGAACTGGTACAACGCAGCCGATTTATGATGCTGGCAGAAGAAGTGCCTGAGGCGCGGGAACATATGGGGCGCTACGGGc                                                                                 <  1:6337718/93‑1 (MQ=255)
        gggTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGc                                                                                 <  2:1625933/93‑1 (MQ=255)
        gggTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGc                                                                                 <  1:3652425/93‑1 (MQ=255)
        gggTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGACGCGCGGGAACATATGGGGCGCTACGGGc                                                                                 >  1:7808128/1‑93 (MQ=255)
         ggTAGCGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCt                                                                                <  2:3134005/93‑1 (MQ=255)
         ggTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCt                                                                                >  2:7071118/1‑93 (MQ=255)
          gTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTg                                                                               >  2:1594645/1‑93 (MQ=255)
           tAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTgg                                                                              >  1:2898335/1‑93 (MQ=255)
           tAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTgg                                                                              >  2:752526/1‑93 (MQ=255)
           gaaCGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTgg                                                                              <  1:1032382/92‑1 (MQ=255)
              cGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGa                                                                           <  2:4632521/93‑1 (MQ=255)
               ggAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGAt                                                                          <  2:4035927/93‑1 (MQ=255)
               ggAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGAt                                                                          <  2:3423693/93‑1 (MQ=255)
               ggAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGAt                                                                          >  2:4994073/1‑93 (MQ=255)
               ggAAAGCGAACGGGTACAACGCAGCCCAGTTATGATGGTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGAt                                                                          <  2:6147463/93‑1 (MQ=255)
                gAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATg                                                                         >  2:7222525/1‑93 (MQ=255)
                gAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATg                                                                         <  1:7138930/93‑1 (MQ=255)
                 aaaGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATgg                                                                        <  1:5221192/93‑1 (MQ=255)
                   aGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTa                                                                      <  2:596397/93‑1 (MQ=255)
                    gCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTAc                                                                     <  2:572855/93‑1 (MQ=255)
                    gCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTAc                                                                     <  1:4846614/93‑1 (MQ=255)
                      gAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACgg                                                                   <  2:4327906/93‑1 (MQ=255)
                      gAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACgg                                                                   >  2:7675108/1‑93 (MQ=255)
                        aCTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCa                                                                 <  1:7284954/93‑1 (MQ=255)
                          tGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGTCGCTACGGGCTGGTGATGGTACGGCaga                                                               >  1:319074/1‑93 (MQ=255)
                           ggTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGTTACGGGCTGGCGATGGTACGGCagag                                                              <  1:3966139/93‑1 (MQ=255)
                           ggTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCagag                                                              <  2:5868330/93‑1 (MQ=255)
                            gTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGt                                                             <  1:5706287/93‑1 (MQ=255)
                              aCAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGa                                                           <  1:5150892/93‑1 (MQ=255)
                              aCAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGGACGGCAGAGTGa                                                           >  2:544841/1‑93 (MQ=255)
                               cAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGAt                                                          <  2:2291188/93‑1 (MQ=255)
                               cAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGAt                                                          >  1:5459025/1‑93 (MQ=255)
                                   gCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACt                                                      >  1:6854364/1‑93 (MQ=255)
                                   gCAGCCAATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACt                                                      >  2:7786836/1‑93 (MQ=255)
                                    cAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTc                                                     <  1:5080794/93‑1 (MQ=255)
                                      gCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAt                                                   <  2:1935965/93‑1 (MQ=255)
                                      gCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAt                                                   <  1:2538066/93‑1 (MQ=255)
                                        ccATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAttt                                                 >  1:6410986/1‑93 (MQ=255)
                                         cATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTg                                                <  1:654833/93‑1 (MQ=255)
                                           tttATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTa                                              >  2:2947516/1‑93 (MQ=255)
                                           tttATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTa                                              >  1:4172164/1‑93 (MQ=255)
                                           tttATGATCCTGCAAGAAGAAGTGCCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTa                                              <  2:7633669/93‑1 (MQ=255)
                                            ttATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTAt                                             <  1:5948692/93‑1 (MQ=255)
                                             tATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTAtt                                            <  2:2392441/93‑1 (MQ=255)
                                               tgatgCTGGCAGAAGAAGTACCTGAGGCGCGGGAGCATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGc                                          >  1:6934753/1‑93 (MQ=255)
                                               tgatgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGc                                          <  1:2028159/93‑1 (MQ=255)
                                               tgatgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGc                                          >  2:3680994/1‑93 (MQ=255)
                                               tgatgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGc                                          <  1:7562773/93‑1 (MQ=255)
                                                  tgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTgg                                       >  2:3812347/1‑93 (MQ=255)
                                                     tGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGa                                    <  2:7272697/93‑1 (MQ=255)
                                                          gaagaaGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACaa                               <  1:6299344/93‑1 (MQ=255)
                                                          gaagaaGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACaa                               >  1:4954116/1‑93 (MQ=255)
                                                          gaagaaGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACaa                               >  1:1499809/1‑93 (MQ=255)
                                                           aagaagTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAAc                              >  2:2843222/1‑93 (MQ=255)
                                                           aagaagTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAAc                              >  1:7066150/1‑93 (MQ=255)
                                                           aagaagTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAAc                              <  2:6420741/93‑1 (MQ=255)
                                                            agaagTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACg                             >  2:7896356/1‑93 (MQ=255)
                                                             gaagTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGa                            <  1:3280287/93‑1 (MQ=255)
                                                               agTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGaaa                          <  1:7426162/93‑1 (MQ=255)
                                                               agTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGaaa                          <  2:4590503/93‑1 (MQ=255)
                                                                gTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAAt                         <  2:5616659/93‑1 (MQ=255)
                                                                  aCCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCt                       >  2:2481371/1‑93 (MQ=255)
                                                                   ccTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTa                      <  2:1076939/93‑1 (MQ=255)
                                                                   ccTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTa                      >  2:3611479/1‑93 (MQ=255)
                                                                   ccTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTa                      <  1:1906683/93‑1 (MQ=255)
                                                                   ccTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTa                      <  1:343240/93‑1 (MQ=255)
                                                                    cTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTAc                     <  1:1348384/93‑1 (MQ=255)
                                                                     tGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCTACACAACGAAATCTACt                    <  1:1794542/93‑1 (MQ=255)
                                                                     tGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACt                    >  1:1391668/1‑93 (MQ=255)
                                                                     tGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACt                    >  1:1744436/1‑93 (MQ=255)
                                                                     tGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACt                    <  2:7138254/93‑1 (MQ=255)
                                                                     tGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACt                    <  2:7322497/93‑1 (MQ=255)
                                                                     tGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGCATTGCTCGTGACACAACGAAATCTACt                    >  2:7195751/1‑93 (MQ=255)
                                                                       aGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGa                  <  1:4483056/93‑1 (MQ=255)
                                                                       aGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGa                  <  1:2701070/93‑1 (MQ=255)
                                                                        ggCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGGCACAACGAAATCTACTGAc                 <  1:1283484/93‑1 (MQ=255)
                                                                        ggCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGAc                 <  2:1357060/93‑1 (MQ=255)
                                                                            cgGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCtt             <  2:7582416/93‑1 (MQ=255)
                                                                             gggAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTa            <  2:3371110/93‑1 (MQ=255)
                                                                              ggAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTaa           <  1:400014/93‑1 (MQ=255)
                                                                               gAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAAc          >  1:2336032/1‑93 (MQ=255)
                                                                                aaCATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAAcc         <  1:2425439/93‑1 (MQ=255)
                                                                                 aCATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCg        <  2:516369/93‑1 (MQ=255)
                                                                                 aCATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCg        >  1:7693796/1‑93 (MQ=255)
                                                                                  cATATGGGGCGCTCCGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTTTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGt       <  2:6736790/93‑1 (MQ=255)
                                                                                  cATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGt       <  1:4661801/93‑1 (MQ=255)
                                                                                  aatatGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGt       <  1:1584924/92‑1 (MQ=255)
                                                                                    tatGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGc     <  1:3859106/93‑1 (MQ=255)
                                                                                    tatGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGc     <  1:2219795/93‑1 (MQ=255)
                                                                                       ggggCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTc  <  2:1838530/93‑1 (MQ=255)
                                                                                       ggggCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTc  >  1:3081387/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGCGCGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGACGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTC  >  NC_000913/4108043‑4108222

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