Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,464,897 G→T 6.6% E143* (GAA→TAA)  paaY → thioesterase required for phenylacetic acid degradation; trimeric; phenylacetate regulatory and detoxification protein; hexapeptide repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,464,8970GT 6.6% 186.0 / 3.7 61E143* (GAA→TAA) paaYthioesterase required for phenylacetic acid degradation; trimeric; phenylacetate regulatory and detoxification protein; hexapeptide repeat protein
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base T (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

ATGCATTAGTGGGAATGAACGCGGTAGTGATGGACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTGGCAGGAAACGCCTGGTATTTGATGAGAAT  >  NC_000913/1464751‑1465043
                                                                                                                                                  |                                                                                                                                                  
aTGCATTAGTGGGAATGAACGCGGTAGTGATGGACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTgca                                                                                                                                                  >  3:284203/1‑147 (MQ=255)
        gTGGGAATGAACGCGGTAGTGATGGACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGa                                                                                                                                          <  8:369920/149‑1 (MQ=255)
           ggAATGAACGCGGTAGTGATGGACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCa                                                                                                                                       <  8:198735/149‑1 (MQ=255)
           ggAATGAACGCGGTAGTGATGGACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGCGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAg                                                                                                                                         >  2:343101/1‑147 (MQ=255)
              aTGAACGCGGTAGTGATGGACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGa                                                                                                                                    <  8:308158/149‑1 (MQ=255)
                     cgGTAGTGATGGACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCa                                                                                                                             <  6:100527/149‑1 (MQ=255)
                        tAGTGATGGACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATgg                                                                                                                          <  8:334614/149‑1 (MQ=255)
                           tGATGGACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGaa                                                                                                                        >  1:364160/1‑148 (MQ=255)
                           tGATGGACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGaa                                                                                                                        >  7:365506/1‑148 (MQ=255)
                             aTGGACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGaaaaa                                                                                                                     <  2:343597/149‑1 (MQ=255)
                             aTGGACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGaaa                                                                                                                       >  6:355274/1‑147 (MQ=255)
                               ggACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGa                                                                                                                                    <  1:340812/132‑1 (MQ=255)
                               ggACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGa                                                                                                                                    >  2:340812/1‑132 (MQ=255)
                                gACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAACTGCCAGGTAATTACCTGATTGGCGGCAGCCGGGCGAAAGCGATGCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGc                                                                                                                   <  1:74858/148‑1 (MQ=255)
                                gACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCa                                                                                                                  <  8:124116/149‑1 (MQ=255)
                                 aCGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAg                                                                                                                 <  8:151866/149‑1 (MQ=255)
                                       cggtgATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTAc                                                                                                            <  1:364727/148‑1 (MQ=255)
                                           gATTGGCGAGAACAGCATTGTTGGTGAATCCGAATTTGTGAAAGCAAAAGCAGAATTGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGAAGGAGTTTGCATGGAAAAAGCAGGGTACGCAt                                                                                                        >  1:354634/1‑148 (MQ=255)
                                            aTTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGa                                                                                                      <  2:364160/149‑1 (MQ=255)
                                                   ataaCAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAg                                                                                               <  6:303070/147‑1 (MQ=255)
                                                   agaACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAg                                                                                               <  2:88444/149‑1 (MQ=255)
                                                   agaACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAg                                                                                               <  3:34956/149‑1 (MQ=255)
                                                   agaACAGCATTGTTGGGGCATCCGCATTTGTGAAACCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAg                                                                                               <  1:313557/149‑1 (MQ=255)
                                                           aTTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCCGCGAAAGCGATTCTTGAACTCACTGAGCAGGAGTTGGCATGGAAAAAGCAGGGAACGCATGAGTCCCAGGTGCgggg                                                                                       >  5:291236/1‑148 (MQ=255)
                                                                 ggTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGacac                                                                                  <  4:284203/148‑1 (MQ=255)
                                                                   tGCATCCGCATTTGCGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGc                                                                                <  1:343101/148‑1 (MQ=255)
                                                                     cATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTAGGCAGCCCGGCGAAAGCGATTCTTGAACTAAGTGAGCAGTAGTTGGCATGGAAAAAGCAGGGTAGGCATGAGTACCAGGTGATGTGGacaccatg                                                                              >  6:232012/1‑144 (MQ=255)
                                                                     cATCCGCAATTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACATGATTGTCGGCAGCCCGGCGAAAGAGATTCGTGAACTCAGTGAGCAGGAGGTGGCATGGTAAAAGCAGGGTACGAAAGAGTACCAGGTGCTTGTGACACGCTg                                                                              >  6:237267/1‑148 (MQ=255)
                                                                       tCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTTAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGt                                                                             >  1:294254/1‑147 (MQ=255)
                                                                       tCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTTAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGt                                                                             >  1:234760/1‑147 (MQ=255)
                                                                         cGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAgg                                                                                                                <  3:16566/110‑1 (MQ=255)
                                                                         cGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAgg                                                                                                                >  4:16566/1‑110 (MQ=255)
                                                                                   aaaGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACa                                                               <  7:228717/149‑1 (MQ=255)
                                                                                     aGCCATAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGATAGCGATTCGTGAACTCAGTGAGCAGTAGTTGGCATGTAAAATGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATc                                                             <  2:354634/149‑1 (MQ=255)
                                                                                       ccAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTTAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCaa                                                           <  2:294254/149‑1 (MQ=255)
                                                                                       ccAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTTAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGAAGTTACATCaa                                                           <  2:234760/149‑1 (MQ=255)
                                                                                        cAAAGCAGAAATGACAGCTAATTACCAGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTAGGCATGGAAAAAGCACGGTACGCATGAGTACCAGGTGCTGGTGCCACGCTGTAAGCAGACGTTACATCAAg                                                          >  4:194400/1‑149 (MQ=255)
                                                                                          aaGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTc                                                        >  1:341266/1‑149 (MQ=255)
                                                                                                     ccAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGGGATTCGTGAACTCAGTGAGCAGGAGTTTGCATGGAAAAAGCAGGGTAGGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCg                                             <  4:23575/149‑1 (MQ=255)
                                                                                                      cAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAATCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGATGTAAGCAGACGTTACATCAAGTCGAGCAATTGCgt                                            <  8:98695/149‑2 (MQ=255)
                                                                                                       aGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCgg                                            >  4:239022/1‑148 (MQ=255)
                                                                                                         cTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTAGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGa                                          >  7:169419/1‑148 (MQ=255)
                                                                                                                ccTGATTGTCGGCAGCCCGGCGAAAGCTATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAACGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGaa                                   <  1:103941/148‑1 (MQ=255)
                                                                                                                ccTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCGGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAAc                                  >  7:133712/1‑149 (MQ=255)
                                                                                                                 cTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAAcc                                 <  8:365506/149‑1 (MQ=255)
                                                                                                                    aTTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTg                               <  5:355274/148‑1 (MQ=255)
                                                                                                                     ctGTCGGCAGCCCGGGGTAAGCGATTCGTGAAGTCAGTGAGTAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGGTGTAAGCAGACGTTAGATCAAGTCGAGCCATTGCGGGAAATTGAACCTgg                              <  8:83602/147‑1 (MQ=255)
                                                                                                                        tCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTGGCAg                           >  7:83602/1‑148 (MQ=255)
                                                                                                                        tCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTGGCAg                           >  3:53011/1‑148 (MQ=255)
                                                                                                                          ggCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGCCGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTGGCAGGa                         >  2:77315/1‑148 (MQ=255)
                                                                                                                          ggCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTGGCAGGa                         >  5:117920/1‑148 (MQ=255)
                                                                                                                            cAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTaa                                                                           <  6:177781/96‑1 (MQ=255)
                                                                                                                            cAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTaa                                                                           >  5:177781/1‑96 (MQ=255)
                                                                                                                                 cGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTGGCAGGAAACGCCTg                 >  8:366321/1‑149 (MQ=255)
                                                                                                                                     gAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTGGCAGGAAACGCCTGGTAt             <  6:117920/149‑1 (MQ=255)
                                                                                                                                      aaaGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTGGCAGGAAACGCCTGGTAtt            >  8:160986/1‑149 (MQ=255)
                                                                                                                                      aaaGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTGGCAGGAAACGCCTGGTAtt            >  5:30774/1‑149 (MQ=255)
                                                                                                                                       aaGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTGGCAGGAAACGCCTGGTAttt           >  6:268980/1‑149 (MQ=255)
                                                                                                                                           gATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTGGCAGGAAACGCCTGGTAtt            >  2:21602/1‑144 (MQ=255)
                                                                                                                                           gATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTGGCAGGAAACGCCTGGTAtt            <  1:21602/144‑1 (MQ=255)
                                                                                                                                           gATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTGGCAGGAAACGCCTGGTATTtgat        >  4:136107/1‑148 (MQ=255)
                                                                                                                                                gTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTGGCAGGAAACGCCTGGTATTTGATGAGAAt  >  2:366182/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATGCATTAGTGGGAATGAACGCGGTAGTGATGGACGGTGCGGTGATTGGCGAGAACAGCATTGTTGGTGCATCCGCATTTGTGAAAGCCAAAGCAGAAATGCCAGCTAATTACCTGATTGTCGGCAGCCCGGCGAAAGCGATTCGTGAACTCAGTGAGCAGGAGTTGGCATGGAAAAAGCAGGGTACGCATGAGTACCAGGTGCTGGTGACACGCTGTAAGCAGACGTTACATCAAGTCGAGCCATTGCGGGAAATTGAACCTGGCAGGAAACGCCTGGTATTTGATGAGAAT  >  NC_000913/1464751‑1465043

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