Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,723,065 G→A A113A (GCG→GCA)  pssA → phosphatidylserine synthase, CDP‑diacylglycerol‑serine O‑phosphatidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,723,0650GA100.0% 271.5 / NA 82A113A (GCG→GCA) pssAphosphatidylserine synthase, CDP‑diacylglycerol‑serine O‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (42/40);  total (42/40)

ATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGATGT  >  NC_000913/2722941‑2723199
                                                                                                                            |                                                                                                                                      
aTGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCggg                                                                                                                          >  2:1051039/1‑139 (MQ=255)
aTGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTAGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCggg                                                                                                                          >  1:1054095/1‑139 (MQ=255)
         aaaGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAAt                                                                                                                               <  1:1332305/125‑1 (MQ=255)
         aaaGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAAt                                                                                                                               >  2:1332305/1‑125 (MQ=255)
         aaaGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGt                                                                                                                 >  1:520883/1‑139 (MQ=255)
         aaaGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGt                                                                                                                 >  1:798404/1‑139 (MQ=255)
         aaaGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGt                                                                                                                 >  1:432976/1‑139 (MQ=255)
                gCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGtt                                                                                                          >  1:28862/1‑139 (MQ=255)
                  gTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTa                                                                                                        <  1:1680718/139‑1 (MQ=255)
                    ccGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATg                                                                                                      <  2:1365669/139‑1 (MQ=255)
                    ccGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATg                                                                                                      <  1:289557/139‑1 (MQ=255)
                    ccGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATg                                                                                                      <  2:28862/139‑1 (MQ=255)
                          cTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTc                                                                                                <  2:798404/139‑1 (MQ=255)
                             gATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATc                                                                                                                              >  2:1222303/1‑106 (MQ=255)
                             gATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATc                                                                                                                              <  1:1222303/106‑1 (MQ=255)
                             gATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaa                                                                                             >  1:1586044/1‑139 (MQ=255)
                                 tgCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaatcaa                                                                                         <  2:432976/139‑1 (MQ=255)
                                    gggTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATAc                                                                                      <  2:1645775/139‑1 (MQ=255)
                                       tGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCg                                                                                   >  2:134442/1‑139 (MQ=255)
                                       tGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCg                                                                                   >  1:652757/1‑139 (MQ=255)
                                        gCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGt                                                                                  <  2:652757/139‑1 (MQ=255)
                                         cTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGtt                                                                                                 >  2:1002978/1‑123 (MQ=255)
                                         cTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGtt                                                                                                 <  1:1002978/123‑1 (MQ=255)
                                           ggTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGtt                                                                                                                <  2:794979/106‑1 (MQ=255)
                                           ggTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGtt                                                                                                                >  1:794979/1‑106 (MQ=255)
                                            gTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGtt                                                                                                                >  1:627350/1‑105 (MQ=255)
                                            gTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGtt                                                                                                                <  2:627350/105‑1 (MQ=255)
                                                aCTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCt                                                                          <  1:70445/139‑1 (MQ=255)
                                                aCTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCt                                                                          <  2:834139/139‑1 (MQ=255)
                                                    gCATCGTGCACAACGTGGACGCTTTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGTCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGgt                                                                      >  2:942569/1‑139 (MQ=255)
                                                    gCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGgt                                                                      >  1:859986/1‑139 (MQ=255)
                                                          tGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGc                                                                             >  1:584982/1‑126 (MQ=255)
                                                          tGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGc                                                                             <  2:584982/126‑1 (MQ=255)
                                                              caACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACt                                                            >  1:941250/1‑139 (MQ=255)
                                                                 cGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTa                                                         <  2:1586044/139‑1 (MQ=255)
                                                                 cGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTa                                                         <  2:859986/139‑1 (MQ=255)
                                                                 cGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTa                                                         >  1:164463/1‑139 (MQ=255)
                                                                 cGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTa                                                         >  1:504917/1‑139 (MQ=255)
                                                                    ggACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAg                                                      >  2:1646782/1‑139 (MQ=255)
                                                                        gCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCtt                                                  <  2:1054095/139‑1 (MQ=255)
                                                                                   gCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATcgacg                                       <  1:134442/139‑1 (MQ=255)
                                                                                    cGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATcgacga                                      >  2:293032/1‑139 (MQ=255)
                                                                                      gCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATa                                    <  2:520883/139‑1 (MQ=255)
                                                                                         tCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCg                                 <  1:293032/139‑1 (MQ=255)
                                                                                                 tAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTg                                                                >  1:660915/1‑100 (MQ=255)
                                                                                                 tAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTg                                                                <  2:660915/100‑1 (MQ=255)
                                                                                                   aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaatc                                                                                           >  2:1410304/1‑71 (MQ=255)
                                                                                                   aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaatc                                                                                           <  1:1410304/71‑1 (MQ=255)
                                                                                                   aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTg                                                                >  1:890978/1‑98 (MQ=255)
                                                                                                   aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTg                                                                <  2:890978/98‑1 (MQ=255)
                                                                                                    cGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTa                                                <  1:1583705/113‑1 (MQ=255)
                                                                                                    cGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTa                                                >  2:1583705/1‑113 (MQ=255)
                                                                                                    cGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGc                      >  2:1320236/1‑139 (MQ=255)
                                                                                                        gACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTg                  <  1:1400709/139‑1 (MQ=255)
                                                                                                        gACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTg                  <  1:1051039/139‑1 (MQ=255)
                                                                                                         aCTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTa                                                         >  1:542557/1‑99 (MQ=255)
                                                                                                         aCTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTa                                                         <  2:542557/99‑1 (MQ=255)
                                                                                                         aCTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGc                 >  2:1348170/1‑139 (MQ=255)
                                                                                                               aCTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCt           >  1:1346559/1‑139 (MQ=255)
                                                                                                                cTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCt                                                                          >  1:738687/1‑75 (MQ=255)
                                                                                                                cTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCt                                                                          <  2:738687/75‑1 (MQ=255)
                                                                                                                 tGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGc                                  >  1:1689395/1‑114 (MQ=255)
                                                                                                                 tGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGc                                  <  2:1689395/114‑1 (MQ=255)
                                                                                                                  gCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGaa        <  1:1320236/139‑1 (MQ=255)
                                                                                                                   ccGCATGGCACAGGAAAATCCGGGCTTAGATGTTCCGGTTTATGGCGTTc                                                                                                <  2:1175737/50‑1 (MQ=255)
                                                                                                                   ccGCATGGCACAGGAAAATCAGGGCGTAGATGTTCCGGTTTATGGCGTTc                                                                                                >  1:1175737/1‑50 (MQ=255)
                                                                                                                    cGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACg      <  2:1323790/139‑1 (MQ=255)
                                                                                                                     gCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTaaa                                                       <  1:715169/89‑1 (MQ=255)
                                                                                                                     gCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTaaa                                                       >  2:715169/1‑89 (MQ=255)
                                                                                                                     gCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTa                                                <  2:1423/96‑1 (MQ=255)
                                                                                                                     gCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTa                                                >  1:1423/1‑96 (MQ=255)
                                                                                                                     gCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGa     >  1:1280427/1‑139 (MQ=255)
                                                                                                                       aTGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATg                                                                                                      >  1:514866/1‑40 (MQ=255)
                                                                                                                       aTGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATg                                                                                                      <  2:514866/40‑1 (MQ=39)
                                                                                                                       aTGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGTCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCGGaa        <  2:1117026/134‑1 (MQ=255)
                                                                                                                       aTGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTa                               <  1:280884/111‑1 (MQ=255)
                                                                                                                       aTGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTa                               >  2:280884/1‑111 (MQ=255)
                                                                                                                       aTGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGaa        >  1:1117026/1‑134 (MQ=255)
                                                                                                                        tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATg                                                                                                      >  2:74903/1‑39 (MQ=38)
                                                                                                                        tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATg                                                                                                      <  1:74903/39‑1 (MQ=38)
                                                                                                                        tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGATGt  >  1:1311492/1‑139 (MQ=255)
                                                                                                                        tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGATGt  >  2:1343677/1‑139 (MQ=255)
                                                                                                                            |                                                                                                                                      
ATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGATGT  >  NC_000913/2722941‑2723199

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