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% 234.8 / NA 72A113A (GCG→GCA) pssAphosphatidylserine synthase, CDP‑diacylglycerol‑serine O‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/40);  total (32/40)

GAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGATGT  >  NC_000913/2722930‑2723199
                                                                                                                                       |                                                                                                                                      
gAAAGCGTTGTCTGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAg                                                                                                                                     <  2:1082718/139‑1 (MQ=255)
  aCGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGa                                                                                                                                   <  2:768369/139‑1 (MQ=255)
      gTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGACATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAAt                                                                                                                               >  2:1184207/1‑139 (MQ=255)
         gTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCg                                                                                                                            >  1:946245/1‑139 (MQ=255)
           aTGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCggg                                                                                                                          <  2:587986/139‑1 (MQ=255)
              aGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGt                                                                                                                       >  2:1360087/1‑139 (MQ=255)
                    aaaGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGt                                                                                                                 >  2:394453/1‑139 (MQ=255)
                    aaaGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGt                                                                                                                 >  1:590167/1‑139 (MQ=255)
                    aaaGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGt                                                                                                                 >  1:326707/1‑139 (MQ=255)
                    aaaGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACCGCCGCATGGCACAGGAAAATCCGGGCGTAGATGt                                                                                                                 >  1:768055/1‑139 (MQ=255)
                     aaGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGtt                                                                                                                <  1:691119/139‑1 (MQ=255)
                               ccGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATg                                                                                                      <  2:1631638/139‑1 (MQ=255)
                                     cTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACGGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTc                                                                                                <  2:186794/139‑1 (MQ=255)
                                        gATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTCCTGCCGCATGGCACAGGAAAATCCGGACGTAGATGTTCCGGTTTATGGCGTTCcaa                                                                                             >  1:748033/1‑139 (MQ=255)
                                        gATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaa                                                                                             >  2:1668899/1‑139 (MQ=255)
                                        gATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaa                                                                                             >  2:510406/1‑139 (MQ=255)
                                         aTGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACCGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaat                                                                                            <  2:768055/139‑1 (MQ=255)
                                             gCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaatcaat                                                                                        >  2:1178832/1‑139 (MQ=255)
                                              cGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATa                                                                                       <  1:1360087/139‑1 (MQ=255)
                                                ggTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGtt                                                                                                                <  2:231680/112‑1 (MQ=255)
                                                ggTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGtt                                                                                                                >  1:231680/1‑112 (MQ=255)
                                                  tGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAg                                                                                                                     >  2:907354/1‑105 (MQ=255)
                                                  tGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAg                                                                                                                     <  1:907354/105‑1 (MQ=255)
                                                   gCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGt                                                                                  <  1:1668899/139‑1 (MQ=255)
                                                    cTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCggg                                                                                                                          <  2:84477/98‑1 (MQ=255)
                                                    cTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCggg                                                                                                                          >  1:84477/1‑98 (MQ=255)
                                                       gTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAg                                                                              >  1:559982/1‑139 (MQ=255)
                                                        tCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGc                                                                             <  2:1284987/139‑1 (MQ=255)
                                                           aCTGACATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTTGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCt                                                                          <  1:1184207/139‑1 (MQ=255)
                                                             tGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCGGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGAGGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTg                                                                        <  1:394453/139‑1 (MQ=255)
                                                                            cGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATc                                                                                                                              >  2:825868/1‑70 (MQ=255)
                                                                            cGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATc                                                                                                                              <  1:825868/70‑1 (MQ=255)
                                                                            cGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTa                                                         <  1:510406/139‑1 (MQ=255)
                                                                             gTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTaa                                                        <  1:1056000/139‑1 (MQ=255)
                                                                             gTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTa                                                         >  1:1630720/1‑138 (MQ=255)
                                                                             gTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTa                                                         <  2:1630720/138‑1 (MQ=255)
                                                                                 aCGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGc                                                    <  2:748033/139‑1 (MQ=255)
                                                                                  cGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAGTCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCt                                                   <  1:1665582/139‑1 (MQ=255)
                                                                                  cGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCt                                                   <  2:326707/139‑1 (MQ=255)
                                                                                     aTTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTa                                                <  2:1627303/139‑1 (MQ=255)
                                                                                      ttGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGa                                                                                <  1:854460/106‑1 (MQ=255)
                                                                                      ttGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGa                                                                                >  2:854460/1‑106 (MQ=255)
                                                                                          cgcTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTatcatc                                           >  1:1279498/1‑139 (MQ=255)
                                                                                               cGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATcgacga                                      >  2:602307/1‑139 (MQ=255)
                                                                                               cGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATcgacga                                      >  1:1230904/1‑139 (MQ=255)
                                                                                               cGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGCGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATcgacga                                      <  1:1178832/139‑1 (MQ=255)
                                                                                                    tCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCg                                 <  2:1291034/139‑1 (MQ=255)
                                                                                                    tCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCg                                 <  1:220618/139‑1 (MQ=255)
                                                                                                    tCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCg                                 <  2:590167/139‑1 (MQ=255)
                                                                                                    tCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCg                                 >  2:220618/1‑139 (MQ=255)
                                                                                                         cacTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTaaa                                                       >  1:1646902/1‑112 (MQ=255)
                                                                                                         cacTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTaaa                                                       <  2:1646902/112‑1 (MQ=255)
                                                                                                         cacTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACtt                            >  2:717291/1‑139 (MQ=255)
                                                                                                          acTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACttt                           <  1:1210318/139‑1 (MQ=255)
                                                                                                              aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGtt                                                                                                                >  1:1564792/1‑50 (MQ=255)
                                                                                                              aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGtt                                                                                                                <  2:1564792/50‑1 (MQ=255)
                                                                                                              aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGtt                                                                                                 >  2:328343/1‑65 (MQ=255)
                                                                                                              aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGtt                                                                                                 <  1:328343/65‑1 (MQ=255)
                                                                                                              aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGt                                <  1:845433/130‑1 (MQ=255)
                                                                                                              aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGt                                >  2:845433/1‑130 (MQ=255)
                                                                                                              aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAg                       <  1:521392/139‑1 (MQ=255)
                                                                                                                   gACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTg                  <  1:424405/139‑1 (MQ=255)
                                                                                                                    aCTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGt                                                                                                           >  1:1137034/1‑49 (MQ=255)
                                                                                                                    aCTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGt                                                                                                           <  2:1137034/49‑1 (MQ=255)
                                                                                                                     cTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGcc                <  2:946245/139‑1 (MQ=255)
                                                                                                                              ccGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGaa        >  2:1642881/1‑138 (MQ=255)
                                                                                                                              ccGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGaa        <  1:1642881/138‑1 (MQ=255)
                                                                                                                                gCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGc                                                               <  2:806961/81‑1 (MQ=255)
                                                                                                                                gCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGc                                                               >  1:806961/1‑81 (MQ=255)
                                                                                                                                 cATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGAt    >  1:1633349/1‑139 (MQ=255)
                                                                                                                                   tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGATGt  <  2:1633349/139‑1 (MQ=255)
                                                                                                                                   tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGATGt  >  1:296721/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GAACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGATGT  >  NC_000913/2722930‑2723199

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