Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,148,639 A→G Y340C (TAC→TGC)  plsX → putative phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,148,6390AG96.5% 292.0 / ‑0.6 87Y340C (TAC→TGC) plsXputative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (39/45);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATA  >  NC_000913/1148544‑1148734
                                                                                               |                                                                                               
cAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCATCGAATTGCCGCTCGCCTGGAATCTGTATGCCCa                                                                                             <  2:1743024/100‑1 (MQ=255)
  aTCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGc                                                                                           <  1:1565109/100‑1 (MQ=255)
   tCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCt                                                                                          <  1:2011673/100‑1 (MQ=255)
   tCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAAGCTGTATGCCCAGCt                                                                                          <  1:967091/100‑1 (MQ=255)
      gCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGt                                                                                       >  1:608168/1‑100 (MQ=255)
         aGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGtttt                                                                                    <  1:1921222/100‑1 (MQ=255)
         aGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGtttn                                                                                    <  1:1010889/100‑2 (MQ=255)
          gCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTg                                                                                   <  2:608166/100‑1 (MQ=255)
            ttttGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAg                                                                                 <  1:344781/100‑1 (MQ=255)
            ttttGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAg                                                                                 <  2:732275/100‑1 (MQ=255)
             tttGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAgc                                                                                <  1:1404955/100‑1 (MQ=255)
             tttGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAgc                                                                                <  1:1936716/100‑1 (MQ=255)
             tttGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAgc                                                                                >  1:1559554/1‑100 (MQ=255)
             tttGCGGTCGCGATTGAACAGGCAGTGCAGG‑TGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAgct                                                                               >  2:212824/1‑100 (MQ=255)
               tGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAgctg                                                                              >  2:2256776/1‑100 (MQ=255)
                gCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAgctgc                                                                             >  1:1624395/1‑100 (MQ=255)
                  ggTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGctgctg                                                                           >  1:1336390/1‑100 (MQ=255)
                  ggTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGctgctg                                                                           >  2:119817/1‑100 (MQ=255)
                   gTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTgg                                                                          >  1:2318860/1‑100 (MQ=255)
                       cgATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACgg                                                                      >  2:1213190/1‑100 (MQ=255)
                        gATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGt                                                                     <  1:1496089/100‑1 (MQ=255)
                         aTTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTg                                                                    >  2:1988654/1‑100 (MQ=255)
                          ttGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTgg                                                                   <  2:243310/100‑1 (MQ=255)
                          ttGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTgg                                                                   <  2:1408595/100‑1 (MQ=255)
                          ttGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTgg                                                                   <  2:1807900/100‑1 (MQ=255)
                           tGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGc                                                                  <  1:90484/100‑1 (MQ=255)
                             aaCAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGGCGGTGGcga                                                                >  2:627666/1‑98 (MQ=255)
                             aaCAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCaa                                                                >  1:1878678/1‑100 (MQ=255)
                             aaCAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCaa                                                                >  1:2425501/1‑100 (MQ=255)
                             aaCAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCaa                                                                >  1:2888568/1‑100 (MQ=255)
                               cAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCaaaa                                                              >  2:2065954/1‑100 (MQ=255)
                               cAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCaaaa                                                              >  1:1933749/1‑100 (MQ=255)
                                  gCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCg                                                           <  1:1693243/100‑1 (MQ=255)
                                   cAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGcag                                                          <  1:1944231/100‑3 (MQ=255)
                                    aGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGa                                                         >  1:2066833/1‑100 (MQ=255)
                                    aGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGa                                                         <  1:703626/100‑1 (MQ=255)
                                    aGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGa                                                         <  1:2606857/100‑1 (MQ=255)
                                    aGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGa                                                         >  1:568768/1‑100 (MQ=255)
                                     gTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGaa                                                        >  2:31984/1‑100 (MQ=255)
                                       gCAGGCGGTGCAGCGACAAGTTCCTCGGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAAct                                                      >  1:1764860/1‑100 (MQ=255)
                                          ggCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTg                                                   <  1:512002/100‑1 (MQ=255)
                                          ggCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTg                                                   >  1:2068075/1‑100 (MQ=255)
                                            cGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCg                                                 <  2:106471/100‑1 (MQ=255)
                                            cGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCg                                                 <  1:242201/100‑1 (MQ=255)
                                             ggTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCgg                                                <  1:202584/100‑1 (MQ=255)
                                                 cAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGc                                            <  2:1211541/100‑1 (MQ=255)
                                                 cAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGc                                            <  1:2623088/100‑1 (MQ=255)
                                                   gCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAg                                          >  1:1641649/1‑100 (MQ=255)
                                                   gCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAg                                          <  2:1064296/100‑1 (MQ=255)
                                                   gCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAg                                          >  1:954884/1‑100 (MQ=255)
                                                     gACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGa                                        <  1:31984/100‑1 (MQ=255)
                                                      aCAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGAc                                       <  2:2786573/100‑1 (MQ=255)
                                                       cAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACg                                      >  2:360540/1‑100 (MQ=255)
                                                        aaGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGc                                     <  2:1228971/100‑1 (MQ=255)
                                                            tCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGcc                                 <  2:1399419/100‑1 (MQ=255)
                                                             ccTCAGCGAATTGCCGCTCGCCTGGAATCCGTATGCCCAGCTGGCTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGcc                                 >  1:2293761/1‑99 (MQ=255)
                                                              cTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAg                               <  1:2249880/100‑1 (MQ=255)
                                                              cTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAg                               >  2:2087080/1‑100 (MQ=255)
                                                              cTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAg                               <  1:2011348/100‑1 (MQ=255)
                                                               tCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGc                              <  1:912737/100‑1 (MQ=255)
                                                               tCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGc                              >  1:2211028/1‑100 (MQ=255)
                                                                cAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCg                             <  2:2092520/100‑1 (MQ=255)
                                                                cAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCg                             <  1:1709505/100‑1 (MQ=255)
                                                                cAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCg                             >  1:552475/1‑100 (MQ=255)
                                                                 aGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGa                            <  1:1831264/100‑1 (MQ=255)
                                                                  gCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGaa                           >  2:906042/1‑100 (MQ=255)
                                                                      aTTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCg                       >  2:1986985/1‑100 (MQ=255)
                                                                       ttGCCGCTCGGCTGGAATCTGTATGCCCAGCTGGTTGTGAGCTGCGGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGc                      <  2:1234698/100‑1 (MQ=255)
                                                                       ttGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGc                      >  1:2487027/1‑100 (MQ=255)
                                                                       ttGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGc                      <  2:2273148/100‑1 (MQ=255)
                                                                         gCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAg                    >  2:1264794/1‑100 (MQ=255)
                                                                          ccGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGt                   <  1:2256789/100‑1 (MQ=255)
                                                                          ccGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGt                   >  1:858337/1‑100 (MQ=255)
                                                                           cGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGtt                  >  1:1439389/1‑100 (MQ=255)
                                                                           cGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGACAGCGAACTCGCAGtt                  >  2:1696740/1‑100 (MQ=255)
                                                                            gCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCa                     <  2:2013345/96‑1 (MQ=255)
                                                                              tCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTg                <  2:453132/99‑1 (MQ=255)
                                                                                 ccTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAg            <  1:205106/100‑1 (MQ=255)
                                                                                 ccTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAg            <  2:1649474/100‑1 (MQ=255)
                                                                                    ggAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGa         <  1:2553481/100‑1 (MQ=255)
                                                                                      aaTCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACg       <  1:2891538/100‑1 (MQ=255)
                                                                                       aTCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACgg      >  1:1966624/1‑100 (MQ=255)
                                                                                          tGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCGGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGtat   >  1:2020638/1‑100 (MQ=255)
                                                                                          tGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGtat   >  2:1049184/1‑100 (MQ=255)
                                                                                           gTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGtata  >  2:1923875/1‑100 (MQ=255)
                                                                                           gTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGtata  <  1:2510107/100‑1 (MQ=255)
                                                                                           gTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGtata  >  1:2657565/1‑100 (MQ=255)
                                                                                               |                                                                                               
CAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATA  >  NC_000913/1148544‑1148734

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