Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 259,346 C→A R252L (CGC→CTC)  phoE ← outer membrane phosphoporin protein E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913259,3460CA97.9% 316.1 / ‑2.4 93R252L (CGC→CTC) phoEouter membrane phosphoporin protein E
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (44/47);  total (45/48)
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.98e-01

CAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTT  >  NC_000913/259249‑259442
                                                                                                 |                                                                                                
cAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAgt                                                                                                <  2:444499/100‑1 (MQ=255)
    gATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGttt                                                                                            <  1:306570/100‑1 (MQ=255)
        gACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGa                                                                                        >  2:2164149/1‑100 (MQ=255)
         aCGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGaa                                                                                       <  2:1479002/100‑1 (MQ=255)
         aCGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGaa                                                                                       >  1:56537/1‑100 (MQ=255)
                ccAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAg                                                                                >  2:1151424/1‑100 (MQ=255)
                ccAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAg                                                                                >  1:2662235/1‑100 (MQ=255)
                ccAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAg                                                                                <  1:2385741/100‑1 (MQ=255)
                 cAAAGTCAAACTGGTATTTAGCGACCGCTTTAAATTTCTGTGTCTTATTGGCAAAGTCGCCAGTTATTGGCGTCCTTTTGAGTGTTTCAGAATAGAAAGt                                                                               <  1:1222116/100‑1 (MQ=255)
                 cAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGt                                                                               <  2:1391768/100‑1 (MQ=255)
                 cAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGt                                                                               <  2:780307/100‑1 (MQ=255)
                 cAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGt                                                                               <  1:919507/100‑1 (MQ=255)
                 cAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGt                                                                               <  1:797251/100‑1 (MQ=255)
                  aaaGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCACTAGTTACTGCCGTCATTTTGAGTGTTTCAGAATAGAAAGt                                                                               <  2:1903664/99‑1 (MQ=255)
                     gTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGcc                                                                           <  1:2499090/100‑1 (MQ=255)
                       cAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAg                                                                         <  2:2766752/100‑1 (MQ=255)
                       cAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAg                                                                         >  2:1131176/1‑100 (MQ=255)
                          aCTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATa                                                                      <  2:1249075/100‑1 (MQ=255)
                          aCTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATa                                                                      <  1:1671668/100‑1 (MQ=255)
                             ggTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAat                                                                   >  2:2684777/1‑100 (MQ=255)
                                aTTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTTATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAAtatt                                                                >  2:1987551/1‑100 (MQ=255)
                                aTTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAAtatt                                                                >  2:1210537/1‑100 (MQ=255)
                                 ttGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAAtatta                                                               >  1:2212116/1‑100 (MQ=255)
                                   gAGCGACCGCGTCAAAGTTCTGTGTCTTATTGGCAAAGGCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGGTGCCAGATAAATattatt                                                             >  2:2595088/1‑100 (MQ=255)
                                    aGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTg                                                            <  1:2692836/100‑1 (MQ=255)
                                    aGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTg                                                            >  2:875187/1‑100 (MQ=255)
                                      tccgatcTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGc                                                          <  2:43694/94‑1 (MQ=255)
                                      cGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGc                                                          >  2:1583198/1‑100 (MQ=255)
                                      cGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGc                                                          <  2:962925/100‑1 (MQ=255)
                                       gACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCa                                                         <  1:784489/100‑1 (MQ=255)
                                        aCCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCAt                                                        >  2:2307926/1‑100 (MQ=255)
                                        aCCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCAt                                                        >  2:2446897/1‑100 (MQ=255)
                                        aCCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCAt                                                        >  1:1516341/1‑100 (MQ=255)
                                         ccGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATc                                                       >  2:878340/1‑100 (MQ=255)
                                          cGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCg                                                      >  2:2380826/1‑100 (MQ=255)
                                           gCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGt                                                     >  2:17616/1‑100 (MQ=255)
                                            cTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTa                                                    >  1:2794589/1‑100 (MQ=255)
                                            cTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTa                                                    <  2:1312710/100‑1 (MQ=255)
                                              tCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCagatcgga                                                  >  1:43694/1‑93 (MQ=255)
                                              tCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTAtt                                                  >  2:2341564/1‑100 (MQ=255)
                                              tCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTAtt                                                  >  1:366935/1‑100 (MQ=255)
                                              tCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTAtt                                                  <  1:2823514/100‑1 (MQ=255)
                                              tCAAAGTTCTGTGTCTTATGTGCAAATCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTAt                                                   <  1:960737/99‑1 (MQ=255)
                                                aaaGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTc                                                <  2:2811101/100‑1 (MQ=255)
                                                 aaGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTc                                                <  2:2596607/99‑1 (MQ=255)
                                                  aGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAg                                              <  2:96973/100‑1 (MQ=255)
                                                     tCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAGGTTGCCAGATAAATATTATTGGCATCGTATTTCAGAcc                                           <  2:1779596/100‑1 (MQ=255)
                                                     tCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGAcc                                           >  2:2503651/1‑100 (MQ=255)
                                                     tCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGAcc                                           <  2:791393/100‑1 (MQ=255)
                                                     tCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGAcc                                           >  1:469741/1‑100 (MQ=255)
                                                          gtCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTg                                      <  1:2078858/100‑1 (MQ=255)
                                                          gtCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTg                                      >  1:468187/1‑100 (MQ=255)
                                                           tCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGc                                     <  1:1210541/100‑1 (MQ=255)
                                                           tCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGc                                     <  2:2598371/100‑1 (MQ=255)
                                                           tCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGc                                     >  1:1233583/1‑100 (MQ=255)
                                                            cTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGcc                                    <  2:2748164/100‑1 (MQ=255)
                                                            cTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGcc                                    <  1:1845432/100‑1 (MQ=255)
                                                             ttATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGccc                                   <  2:1910426/100‑1 (MQ=255)
                                                             ttATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGcc                                    <  2:2840758/99‑1 (MQ=255)
                                                              tATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCa                                  >  2:2244475/1‑100 (MQ=255)
                                                              tATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCa                                  <  1:2500898/100‑1 (MQ=255)
                                                                ttGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATg                                <  2:541287/100‑1 (MQ=255)
                                                                 tGGCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGc                               >  1:2231633/1‑100 (MQ=255)
                                                                  ggCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCt                              >  1:1754393/1‑100 (MQ=255)
                                                                  ggCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCt                              <  2:1661374/100‑1 (MQ=255)
                                                                  ggCAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCt                              <  1:889591/100‑1 (MQ=255)
                                                                    cAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTc                            >  2:722586/1‑100 (MQ=255)
                                                                    cAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTc                            >  2:720217/1‑100 (MQ=255)
                                                                    cAAAGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTc                            <  1:2446913/100‑1 (MQ=255)
                                                                       aGCCGCCAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGc                         >  1:2723779/1‑100 (MQ=255)
                                                                        gccgccAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCa                        <  1:2378452/100‑1 (MQ=255)
                                                                        gccgccAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCa                        <  1:2468031/100‑1 (MQ=255)
                                                                         ccgccAGTTATTGGCGTTATTTTGAGTGTTTCAGAATAGTAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGGTGCCCATGCTTCTGCa                        <  2:817060/99‑1 (MQ=255)
                                                                          cgccAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACg                      >  1:2046390/1‑100 (MQ=255)
                                                                          cgccAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACg                      <  2:1598741/100‑1 (MQ=255)
                                                                           gccAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGc                     >  1:632931/1‑100 (MQ=255)
                                                                           gccAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGc                     >  1:1283239/1‑100 (MQ=255)
                                                                             cAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCtt                   >  1:758906/1‑100 (MQ=255)
                                                                             cAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCtt                   >  2:535165/1‑100 (MQ=255)
                                                                             cAGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTACACACtt                   >  2:2313563/1‑100 (MQ=255)
                                                                              aGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTg                  <  1:596184/100‑1 (MQ=255)
                                                                              aGTTATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTg                  <  2:2794568/100‑1 (MQ=255)
                                                                                ttATTGGCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGcc                <  2:299642/100‑1 (MQ=255)
                                                                                      gCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGcc          <  1:2088172/100‑1 (MQ=255)
                                                                                      gCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGcc          >  2:532122/1‑100 (MQ=255)
                                                                                      gCGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGcc          >  1:2097685/1‑100 (MQ=255)
                                                                                       cGTCATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCa         >  2:713342/1‑100 (MQ=255)
                                                                                          cATTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACgg      <  2:738895/100‑1 (MQ=255)
                                                                                           aTTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGc     >  2:2573071/1‑100 (MQ=255)
                                                                                           aTTTTGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGc     >  2:354299/1‑100 (MQ=255)
                                                                                            ttttGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCt    >  2:1933110/1‑100 (MQ=255)
                                                                                             tttGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCtt   <  1:1583207/100‑1 (MQ=255)
                                                                                              ttGAGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCttt  >  1:1736394/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CAGCGATGGACGCAGACCAAAGTCAAACTGGTATTGAGCGACCGCTTCAAAGTTCTGTGTCTTATTGGCAAAGCCGCCAGTTATTGGCGTCATTTTGCGTGTTTCAGAATAGAAAGTTGCCAGATAAATATTATTGGCATCGTATTTCAGACCTGTTGCCCATGCTTCTGCACGCTTGCCTGTGCCACGGCTTT  >  NC_000913/259249‑259442

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