Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,387,877 A→G L21P (CTG→CCG)  psd ← phosphatidylserine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,387,8770AG86.1% 662.6 / 90.6 253L21P (CTG→CCG) psdphosphatidylserine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (14/21);  new base G (105/112);  total (120/133)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

Only 100 of 265 total aligned reads displayed.

CCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCC  >  NC_000913/4387785‑4387969
                                                                                            |                                                                                            
ccTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCa                                                                                              >  1:6104067/1‑93 (MQ=255)
  tctTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCggg                                                                                            <  2:1951327/93‑1 (MQ=255)
   ctTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGc                                                                                           >  2:4511709/1‑93 (MQ=255)
     ttCATGTCGACCTTGTAGTATTTAACGACCACATCGATAACCAGTGTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGa                                                                                         <  1:1696750/93‑1 (MQ=255)
      tCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAg                                                                                        >  2:2647737/1‑93 (MQ=255)
      tCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAg                                                                                        <  2:1782335/93‑1 (MQ=255)
        aTGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTa                                                                                      <  2:827503/93‑1 (MQ=255)
         tGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTaa                                                                                     <  2:4719922/93‑1 (MQ=255)
             gACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCa                                                                                 >  2:2912541/1‑93 (MQ=255)
             gACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCa                                                                                 >  2:2069057/1‑93 (MQ=255)
              aCCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAATAAGCCAt                                                                                <  1:6501597/93‑1 (MQ=255)
              aCCTTGTAGCATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCAt                                                                                <  2:5632855/93‑1 (MQ=255)
               ccTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATa                                                                               >  1:5772437/1‑93 (MQ=255)
                 ttGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGt                                                                             >  2:4269603/1‑93 (MQ=255)
                    tcgtaTTTAACGAACAGATCGATACCCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCACCCCGCCGGGCGAGTAAGCCATAGTTTc                                                                          <  1:1310942/91‑1 (MQ=255)
                    tagtaTTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTc                                                                          <  2:5209847/93‑1 (MQ=255)
                     agtaTTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCg                                                                         <  2:3719046/93‑1 (MQ=255)
                     agtaTTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCg                                                                         >  2:2669732/1‑93 (MQ=255)
                       taTTTAACGAACAGATCGATAACCAGTTTTATCAGCCATCCAGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGc                                                                       <  1:6653822/93‑1 (MQ=255)
                         tttAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAg                                                                     >  1:6548230/1‑93 (MQ=255)
                         tttAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAg                                                                     >  2:117700/1‑93 (MQ=255)
                         tttAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAg                                                                     <  1:594268/93‑1 (MQ=255)
                          ttAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGa                                                                    <  1:2836680/93‑1 (MQ=255)
                          ttAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGa                                                                    <  1:422726/93‑1 (MQ=255)
                          ttAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGa                                                                    <  1:4189823/93‑1 (MQ=255)
                           tAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGaa                                                                   >  1:3372628/1‑93 (MQ=255)
                                aaCAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTAc                                                              <  2:5052659/93‑1 (MQ=255)
                                aaCAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTAc                                                              >  1:2776420/1‑93 (MQ=255)
                                aaCAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTAc                                                              >  1:3270880/1‑93 (MQ=255)
                                aaCAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTAc                                                              <  1:2986392/93‑1 (MQ=255)
                                aaCAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTAc                                                              >  2:272904/1‑93 (MQ=255)
                                aaCAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTAc                                                              <  2:3893444/93‑1 (MQ=255)
                                 aCAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACt                                                             <  2:2845633/93‑1 (MQ=255)
                                 aCAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACt                                                             >  2:4927162/1‑93 (MQ=255)
                                 aCAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACt                                                             >  2:3496250/1‑93 (MQ=255)
                                   aGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGt                                                           >  2:1247034/1‑93 (MQ=255)
                                   aGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGt                                                           >  1:5720418/1‑93 (MQ=255)
                                   aGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGt                                                           <  1:2712451/93‑1 (MQ=255)
                                   aGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGt                                                           >  1:1421852/1‑93 (MQ=255)
                                    gATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTa                                                          <  2:5790030/93‑1 (MQ=255)
                                      tCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGc                                                        <  1:4329628/93‑1 (MQ=255)
                                       cGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTcgcg                                                       >  1:1332173/1‑93 (MQ=255)
                                       cGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCg                                                       >  2:226397/1‑93 (MQ=255)
                                        gATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGa                                                      <  2:842113/93‑1 (MQ=255)
                                         aTAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGaa                                                     >  1:5101204/1‑93 (MQ=255)
                                         aTAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGaa                                                     <  1:4006485/93‑1 (MQ=255)
                                          tAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGaaa                                                    <  1:221236/93‑1 (MQ=255)
                                           gaCCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAg                                                   >  2:6268461/2‑93 (MQ=255)
                                           aaCCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAg                                                   <  1:1752798/93‑1 (MQ=255)
                                              cAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGttt                                                >  1:4849980/1‑93 (MQ=255)
                                              cAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGttt                                                <  1:6242297/93‑1 (MQ=255)
                                              cAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGttt                                                >  2:5456750/1‑93 (MQ=255)
                                               aGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTa                                               >  2:3446153/1‑93 (MQ=255)
                                                gTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTaa                                              >  1:6150296/1‑93 (MQ=255)
                                                gTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTaa                                              >  2:6389723/1‑93 (MQ=255)
                                                 ttttGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTaaa                                             >  2:4606769/1‑93 (MQ=255)
                                                 ttttGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTaaa                                             <  2:2049716/93‑1 (MQ=255)
                                                  tttGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAAt                                            <  2:4391792/93‑1 (MQ=255)
                                                  tttGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAAt                                            >  2:4472103/1‑93 (MQ=255)
                                                      tCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAAt                                        <  2:2601182/93‑1 (MQ=255)
                                                       cAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCCCCAGGCGGGTAAGCCATAGTTTCGGCGGAATGTACTGTAGCGAAAGTTTAAATGAAtt                                       >  1:943138/1‑93 (MQ=255)
                                                          ccATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTaa                                    <  1:865825/93‑1 (MQ=255)
                                                           cATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAAc                                   <  2:1751/93‑1 (MQ=255)
                                                           cATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAAc                                   >  2:4808458/1‑93 (MQ=255)
                                                           cATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAAc                                   <  1:4031172/93‑1 (MQ=255)
                                                            aTCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACa                                  <  1:2076035/93‑1 (MQ=255)
                                                            aTCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACa                                  >  1:2291187/1‑93 (MQ=255)
                                                            aTCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACa                                  >  1:55025/1‑93 (MQ=255)
                                                            aTCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACa                                  >  1:789427/1‑93 (MQ=255)
                                                             tCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACaa                                 <  2:595366/93‑1 (MQ=255)
                                                             tCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACaa                                 <  1:4604885/93‑1 (MQ=255)
                                                             tCCTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACaa                                 <  2:38065/93‑1 (MQ=255)
                                                              ccTGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAg                                <  1:2568187/93‑1 (MQ=255)
                                                                tGCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGt                              <  2:6611059/93‑1 (MQ=255)
                                                                 gCCCGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTa                             <  2:3110739/93‑1 (MQ=255)
                                                                  cccGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAg                            <  1:1035352/93‑1 (MQ=255)
                                                                   ccGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGc                           >  1:4725183/1‑93 (MQ=255)
                                                                    cGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGcc                          <  1:5135505/93‑1 (MQ=255)
                                                                    cGCTTGCTTGCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGcc                          <  1:3552522/93‑1 (MQ=255)
                                                                    cGCTTGCTTGCGCCCCAACCCGCCGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCt                         >  2:3008803/1‑93 (MQ=255)
                                                                      cttgcttgCGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTAGGCAGAATGTACTGTAGCGAAAGATTAAATGAATTTAACAAGGTAGCCTc                        >  2:3554755/1‑93 (MQ=255)
                                                                      cttgcttgCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTc                        >  1:5617194/1‑93 (MQ=255)
                                                                       ttgcttgcGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTcc                       <  2:1038095/93‑1 (MQ=255)
                                                                           ttgcGCCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGGAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGc                   >  1:5405310/1‑93 (MQ=255)
                                                                             gcgcCCCGACCCGCCGGGCGAGTAAGCCATCGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCCCCAGGCCa                 <  1:3457370/93‑1 (MQ=255)
                                                                               gcCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCAtt               >  2:3189767/1‑93 (MQ=255)
                                                                               gcCCCAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCAtt               >  1:3740414/1‑93 (MQ=255)
                                                                               gcCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCAtt               <  1:4973070/93‑1 (MQ=255)
                                                                                ccccAACCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTg              <  2:365865/93‑1 (MQ=255)
                                                                                    aaCCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGtttt          <  1:1025807/93‑1 (MQ=255)
                                                                                     aCCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTg         >  2:1842166/1‑93 (MQ=255)
                                                                                     aCCCGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTg         <  1:3018064/93‑1 (MQ=255)
                                                                                     aCCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTg         >  1:3898141/1‑93 (MQ=255)
                                                                                     aCCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTg         <  2:2529800/93‑1 (MQ=255)
                                                                                       ccgccgGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTc       <  1:4087866/93‑1 (MQ=255)
                                                                                        cGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCg      >  1:6308353/1‑93 (MQ=255)
                                                                                        cGCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCg      <  2:1399267/93‑1 (MQ=255)
                                                                                         gCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGt     >  2:501757/1‑93 (MQ=255)
                                                                                         gCCGGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGt     >  1:5882111/1‑93 (MQ=255)
                                                                                           cAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGACGTTc   >  2:1628650/1‑93 (MQ=255)
                                                                                            |                                                                                            
CCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCC  >  NC_000913/4387785‑4387969

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: