Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,795,165 G→T 8.2% F441L (TTC→TTA) ‡ pvdH ← diaminobutyrate‑2‑oxoglutarate transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,795,1650GT 8.2% 254.7 / 8.6 85F441L (TTC→TTA) ‡pvdHdiaminobutyrate‑2‑oxoglutarate transaminase
Reads supporting (aligned to +/- strand):  ref base G (38/40);  new base T (2/5);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

ACTCATGTTTCGCTCTCCTTGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGA  >  NC_002947/4795018‑4795309
                                                                                                                                                   |                                                                                                                                                
acTCATGTTTCGCTCTCCTTGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGaa                                                                                                                                                >  2:1032924/1‑150 (MQ=255)
acTCATGTTTCGCTCTCCTTGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGaa                                                                                                                                                <  1:1384563/150‑1 (MQ=255)
       tttCGCTCTCCTTGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGACGCGCggg                                                                                                                                         >  1:965098/1‑150 (MQ=255)
       tttCGCTCTCCTTGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTgg                                                                                                                                         >  1:1991103/1‑150 (MQ=255)
        ttCGCTCTCCTTTTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTTCAGCGCTgca                                                                                                                                        <  2:803806/150‑3 (MQ=255)
        ttCGCTCTCCTTGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTTCAGCGCTgca                                                                                                                                        <  1:435621/150‑3 (MQ=255)
        ttCGCTCTCCGTGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTTCAGCGCTgca                                                                                                                                        <  1:608684/150‑3 (MQ=255)
          cGCTCTCCTTGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCa                                                                                                                                      <  2:2159482/150‑1 (MQ=255)
            ctctCCTTGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAAc                                                                                                                                    >  2:1155628/1‑150 (MQ=255)
             tctcCTTGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAAcc                                                                                                                                   <  1:934748/150‑1 (MQ=255)
               tcCTTGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCtc                                                                                                                                 >  2:1723355/1‑150 (MQ=255)
                ctttGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCtcg                                                                                                                                <  2:1165580/148‑1 (MQ=255)
                  ttGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCTGCATCGAAAATGCCCTAGGTCTGCTTCGACTGCAGCGCTGAAGCGCTGGGCAACCtcgtc                                                                                                                              <  1:656994/150‑1 (MQ=255)
                  ttGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCtcgtc                                                                                                                              <  2:41912/150‑1 (MQ=255)
                  ttGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCtcgtc                                                                                                                              >  2:1875833/1‑150 (MQ=255)
                     tCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGAt                                                                                                                           >  1:813321/1‑150 (MQ=255)
                      cGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTgg                                                                                                                                         <  2:316418/135‑1 (MQ=255)
                      cGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTgg                                                                                                                                         >  1:316418/1‑135 (MQ=255)
                           ggCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTc                                                                                                                     <  2:157025/150‑1 (MQ=255)
                           ggCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTc                                                                                                                     <  2:813321/150‑1 (MQ=255)
                            gCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCg                                                                                                                    >  2:995352/1‑150 (MQ=255)
                             cTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCgg                                                                                                                   >  2:1499565/1‑150 (MQ=255)
                              tCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGc                                                                                                                  >  1:786846/1‑150 (MQ=255)
                                cGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCAc                                                                                                                <  2:1541823/150‑1 (MQ=255)
                                   agagACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACtga                                                                                                             <  2:1816541/150‑1 (MQ=255)
                                    gagaCCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACtgat                                                                                                            <  1:1723355/150‑1 (MQ=255)
                                     agaCCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACtgatg                                                                                                           <  2:189175/150‑1 (MQ=255)
                                       aCCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTgatgat                                                                                                         <  1:1032924/150‑1 (MQ=255)
                                            cGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAAcg                                                                                                    <  2:568974/150‑1 (MQ=255)
                                                   cTTGGGTTGGCGGACAGACTTGCAGGCCGACTGCGTCACACAGTGGCATGAAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAg                                                                                             <  2:127517/150‑1 (MQ=255)
                                                    ttGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAgg                                                                                            >  1:1134512/1‑150 (MQ=255)
                                                      gggTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTTCAGCGCTgca                                                                                                                                        >  2:805730/1‑102 (MQ=255)
                                                      gggTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTTCAGCGCTgca                                                                                                                                        <  1:805730/104‑3 (MQ=255)
                                                        gTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAAc                                                                                        >  2:1028111/1‑150 (MQ=255)
                                                        gTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAAc                                                                                        <  2:1211202/150‑1 (MQ=255)
                                                         ttGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACg                                                                                       >  1:1760943/1‑150 (MQ=255)
                                                          tGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGc                                                                                      <  1:1206578/150‑1 (MQ=255)
                                                              ggcCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGcacca                                                                                  <  1:1875833/150‑1 (MQ=255)
                                                               gcCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACtgat                                                                                                            >  1:1714060/1‑123 (MQ=255)
                                                               gcCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACtgat                                                                                                            <  2:1714060/123‑1 (MQ=255)
                                                                      ttGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCa                                                                          >  1:1346853/1‑150 (MQ=255)
                                                                      ttGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCa                                                                          >  1:752764/1‑150 (MQ=255)
                                                                        gCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATc                                                                                                                          <  2:352684/100‑1 (MQ=255)
                                                                        gCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATc                                                                                                                          >  1:352684/1‑100 (MQ=255)
                                                                        gCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTgatgat                                                                                                         <  1:274717/117‑1 (MQ=255)
                                                                        gCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTgatgat                                                                                                         >  2:274717/1‑117 (MQ=255)
                                                                          aGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCtcgt                                                                                                                               <  1:1366018/93‑1 (MQ=255)
                                                                          aGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCtcgt                                                                                                                               >  2:1366018/1‑93 (MQ=255)
                                                                          aGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATgcc                                                                      >  2:1465941/1‑150 (MQ=255)
                                                                           ggCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATgccg                                                                     >  1:1240897/1‑150 (MQ=255)
                                                                           ggCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATgccg                                                                     >  2:2176340/1‑150 (MQ=255)
                                                                           ggCCGACCGCGTCAAACAGTAGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATgccg                                                                     <  1:1155628/150‑1 (MQ=255)
                                                                            gCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATgccgc                                                                    >  2:1030273/1‑150 (MQ=255)
                                                                               gACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGccca                                                                 <  2:1335791/150‑1 (MQ=255)
                                                                                    cgTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCaa                                                            <  1:1028111/150‑1 (MQ=255)
                                                                                     gTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAAc                                                           <  2:752764/150‑1 (MQ=255)
                                                                                     gTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACAGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAAc                                                           >  1:1599040/1‑150 (MQ=255)
                                                                                        aaaCAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTcc                                                        <  1:1465941/150‑1 (MQ=255)
                                                                                         aaCAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCa                                                       <  2:1346853/150‑1 (MQ=255)
                                                                                         aaCAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCa                                                       <  2:1991103/150‑1 (MQ=255)
                                                                                                 cATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCTGTTGCTGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGc                                               <  1:1499565/150‑1 (MQ=255)
                                                                                                 cATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGc                                               >  1:1340373/1‑150 (MQ=255)
                                                                                                 cATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGc                                               <  2:1134512/150‑1 (MQ=255)
                                                                                                   tGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGccc                                             >  2:2090457/1‑150 (MQ=255)
                                                                                                        cgcgGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAgg                                                     <  2:1804027/137‑1 (MQ=255)
                                                                                                        cgcgGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAgg                                                     >  1:1804027/1‑137 (MQ=255)
                                                                                                         gcggcgTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCtt                                       >  2:503615/1‑150 (MQ=255)
                                                                                                               tCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCa                                 >  1:665758/1‑150 (MQ=255)
                                                                                                                    cGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTcgc                            >  2:1599108/1‑150 (MQ=255)
                                                                                                                      cTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCgcgc                          <  2:490591/150‑1 (MQ=255)
                                                                                                                         ccTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGc                       <  1:369703/150‑1 (MQ=255)
                                                                                                                           tACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCa                                                       <  1:1668792/116‑1 (MQ=255)
                                                                                                                           tACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCa                                                       >  2:1668792/1‑116 (MQ=255)
                                                                                                                           tACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCAc                     >  2:1647856/1‑150 (MQ=255)
                                                                                                                               gcggcggcgACTGCAGCGCTGGAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCAGAGCTTGAGGCACTCGCGCTGCACCTTc                 <  2:1240897/150‑1 (MQ=255)
                                                                                                                               gcggcggcgACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTc                 <  1:2176340/150‑1 (MQ=255)
                                                                                                                                    gtcgACTGCAGCGCGGAATCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGgcgc            >  2:1844097/3‑150 (MQ=255)
                                                                                                                                       gACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCaa         >  2:2189655/1‑150 (MQ=255)
                                                                                                                                       gACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCaa         >  2:1921527/1‑150 (MQ=255)
                                                                                                                                           gCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTcg     <  1:1921527/150‑1 (MQ=255)
                                                                                                                                           ccAGCGCTTCAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAg                                    <  1:1176819/118‑1 (MQ=255)
                                                                                                                                           ccAGCGCTTCAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAg                                    >  2:1176819/2‑119 (MQ=255)
                                                                                                                                             aGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTcgcg   >  2:267574/1‑150 (MQ=255)
                                                                                                                                                 cTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGa  <  2:478480/147‑1 (MQ=255)
                                                                                                                                                 cTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGa  >  1:478480/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ACTCATGTTTCGCTCTCCTTGTCGCCAGGCTCCGTAGAGACCGTCGACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGA  >  NC_002947/4795018‑4795309

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: