Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 656,995 C→T L147F (CTC→TTC)  pagP → Lipid IVA palmitoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913656,9950CT99.7% 1356.3 / ‑4.6 317L147F (CTC→TTC) pagPLipid IVA palmitoyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (180/135);  total (180/137)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 319 total aligned reads displayed.

GGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTT  >  NC_000913/656897‑657092
                                                                                                  |                                                                                                 
  gAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTccc                                                                                               <  2:1231648/101‑1 (MQ=255)
   aaaGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCg                                                                                              >  2:1232423/1‑101 (MQ=255)
    aaGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCgg                                                                                             >  2:1232427/1‑101 (MQ=255)
    aaGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCgg                                                                                             >  1:1232425/1‑101 (MQ=255)
    aaGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCgg                                                                                             <  1:1231965/101‑1 (MQ=255)
      gTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGtt                                                                                           >  1:1232432/1‑101 (MQ=255)
        aCCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCt                                                                                         >  1:1232435/1‑101 (MQ=255)
        aCCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCt                                                                                         <  1:1231807/101‑1 (MQ=255)
         ccTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTa                                                                                        >  2:1232439/1‑101 (MQ=255)
               gACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCa                                                                                  >  1:1232443/1‑101 (MQ=255)
                aCCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCAt                                                                                 >  2:1232444/1‑101 (MQ=255)
                 ccGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCAtt                                                                                <  1:1232102/101‑1 (MQ=255)
                 ccGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCAtt                                                                                >  1:1232447/1‑101 (MQ=255)
                  cGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTg                                                                               <  1:1231979/101‑1 (MQ=255)
                  cGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTg                                                                               <  1:1231997/101‑1 (MQ=255)
                   gCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTgg                                                                              >  2:1232453/1‑101 (MQ=255)
                    cTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGc                                                                             >  1:1232454/1‑101 (MQ=255)
                     tGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACTCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGcc                                                                            <  2:1232186/101‑1 (MQ=255)
                     tGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGcc                                                                            >  1:1232461/1‑101 (MQ=255)
                     tGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGcc                                                                            >  2:1232459/1‑101 (MQ=255)
                      ggCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCt                                                                           >  1:1232465/1‑101 (MQ=255)
                        cGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTcc                                                                         <  1:1231337/101‑1 (MQ=255)
                         ggATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCg                                                                        >  2:1232469/1‑101 (MQ=255)
                         ggATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCg                                                                        <  1:1231717/101‑1 (MQ=255)
                         ggATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCg                                                                        <  1:1231757/101‑1 (MQ=255)
                          gATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCg                                                                        <  1:1231869/100‑1 (MQ=255)
                          gATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGt                                                                       >  1:1232471/1‑101 (MQ=255)
                           aTGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTg                                                                      >  1:1232474/1‑101 (MQ=255)
                           aTGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTg                                                                      <  2:1231701/101‑1 (MQ=255)
                            tGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTgg                                                                     >  2:1232479/1‑101 (MQ=255)
                              aaaaTTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGt                                                                   >  1:1232482/1‑101 (MQ=255)
                               aaaTTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGtt                                                                  <  1:1231828/101‑1 (MQ=255)
                                aaTTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTa                                                                 >  2:1232483/1‑101 (MQ=255)
                                 aTTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTAt                                                                <  2:1231882/101‑1 (MQ=255)
                                 aTTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTAt                                                                >  2:1232487/1‑101 (MQ=255)
                                  ttttCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATg                                                               <  2:1231989/101‑1 (MQ=255)
                                   tttCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGt                                                                   >  1:1232491/1‑96 (MQ=255)
                                      cATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGccc                                                           >  2:1232494/1‑101 (MQ=255)
                                         ttAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGcc                                                            >  2:1232498/1‑97 (MQ=255)
                                         ttAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGt                                                        <  2:1232327/101‑1 (MQ=255)
                                          tAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTg                                                       >  2:1232503/1‑101 (MQ=255)
                                          tAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTg                                                       >  1:1232501/1‑101 (MQ=255)
                                          tAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTg                                                       <  1:1231914/101‑1 (MQ=255)
                                            ggTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGAc                                                     >  2:1232508/1‑101 (MQ=255)
                                            ggTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGAc                                                     >  1:1232507/1‑101 (MQ=255)
                                              tCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACtt                                                   <  1:1231855/101‑1 (MQ=255)
                                                tGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACtttt                                                 >  2:1232513/1‑101 (MQ=255)
                                                  ggATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCa                                               >  1:1232516/1‑101 (MQ=255)
                                                   gATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAg                                              >  2:1232520/1‑101 (MQ=255)
                                                   gATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAg                                              <  2:1232034/101‑1 (MQ=255)
                                                   gATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAg                                              >  1:1232522/1‑101 (MQ=255)
                                                    aTTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTc                                                                          <  1:1231779/72‑1 (MQ=255)
                                                    aTTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGa                                             >  2:1232526/1‑101 (MQ=255)
                                                    aTTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGa                                             <  2:1232023/101‑1 (MQ=255)
                                                       cACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGa                                          <  2:1231896/101‑1 (MQ=255)
                                                          cGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGAt                                            >  2:1232530/1‑96 (MQ=255)
                                                           gCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTa                                      <  2:1231861/101‑1 (MQ=255)
                                                           gCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTa                                      >  1:1232531/1‑101 (MQ=255)
                                                            cTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTAc                                     <  2:1231909/101‑1 (MQ=255)
                                                            cTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTAc                                     <  2:1232063/101‑1 (MQ=255)
                                                             tGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACa                                    >  2:1232536/1‑101 (MQ=255)
                                                             tGGCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACa                                    <  2:1232146/101‑1 (MQ=255)
                                                               gCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGa                                             >  2:1232539/1‑90 (MQ=255)
                                                               gCGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACAtt                                  <  1:1232164/101‑1 (MQ=255)
                                                                cGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTc                                 <  2:1231952/101‑1 (MQ=255)
                                                                cGTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTc                                 >  1:1232542/1‑101 (MQ=255)
                                                                 gTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTcc                                >  2:1232545/1‑101 (MQ=255)
                                                                 gTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTcc                                <  2:1231804/101‑1 (MQ=255)
                                                                 gTAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTcc                                <  1:1232122/101‑1 (MQ=255)
                                                                  tAACGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCg                               >  1:1232547/1‑101 (MQ=255)
                                                                     cGGCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGt                            >  2:1232552/1‑101 (MQ=255)
                                                                       gCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTAc                          <  2:1231752/101‑1 (MQ=255)
                                                                       gCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTAc                          <  1:1231910/101‑1 (MQ=255)
                                                                       gCACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTAc                          >  1:1232554/1‑101 (MQ=255)
                                                                        cACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTAcc                         <  2:1231930/101‑1 (MQ=255)
                                                                        cACGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTAcc                         >  2:1232556/1‑101 (MQ=255)
                                                                         aCGCGATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCt                        <  1:1232082/101‑1 (MQ=255)
                                                                           gcgATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTac                      >  2:1232561/1‑101 (MQ=255)
                                                                           gcgATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTac                      >  1:1232562/1‑101 (MQ=255)
                                                                            cgATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTaca                     <  2:1232002/101‑1 (MQ=255)
                                                                             gATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTacaa                    >  2:1232565/1‑101 (MQ=255)
                                                                             gATAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTacaa                    >  1:1232572/1‑101 (MQ=255)
                                                                              aTAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTacaac                   >  2:1232574/1‑101 (MQ=255)
                                                                               tAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTacaaca                  <  2:1232173/101‑1 (MQ=255)
                                                                               tAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTacaaca                  >  1:1232578/1‑101 (MQ=255)
                                                                               tAACTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTacaaca                  >  1:1232581/1‑101 (MQ=255)
                                                                                aaCTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTAcaacaa                 >  2:1232583/1‑101 (MQ=255)
                                                                                  cTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATg               <  1:1232010/101‑1 (MQ=255)
                                                                                  cTGGAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATg               >  2:1232587/1‑101 (MQ=255)
                                                                                     gAATTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCa            <  1:1232278/101‑1 (MQ=255)
                                                                                      aaTTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCaa           <  2:1232100/101‑1 (MQ=255)
                                                                                      aaTTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCaa           >  2:1232590/1‑101 (MQ=255)
                                                                                       aTTACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAAt          <  1:1232046/101‑1 (MQ=255)
                                                                                         tACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTacaa                    <  2:1232077/89‑1 (MQ=255)
                                                                                         tACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAAtgt        <  1:1231931/101‑1 (MQ=255)
                                                                                         tACATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAAtgt        >  1:1232592/1‑101 (MQ=255)
                                                                                          aCATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAAtgtg       >  2:1232594/1‑101 (MQ=255)
                                                                                           cATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATgtgt      <  1:1231884/101‑1 (MQ=255)
                                                                                           cATCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATgtgt      <  2:1232223/101‑1 (MQ=255)
                                                                                             tCCCTTTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAg                                              >  1:1232597/1‑59 (MQ=255)
                                                                                                  |                                                                                                 
GGGAAAGTACCTGGCGACCGCTGGCGGATGAAAATTTTCATTTAGGTCTGGGATTCACCGCTGGCGTAACGGCACGCGATAACTGGAATTACATCCCTCTCCCGGTTCTACTGCCATTGGCCTCCGTGGGTTATGGCCCAGTGACTTTTCAGATGACCTACATTCCGGGTACCTACAACAATGGCAATGTGTACTT  >  NC_000913/656897‑657092

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