Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,354,905 G→T 6.5% T189N (ACC→AAC)  sapB ← antimicrobial peptide transport ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,354,9050GT 6.5% 212.9 / 4.8 62T189N (ACC→AAC) sapBantimicrobial peptide transport ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base T (2/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

ACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACC  >  NC_000913/1354759‑1355052
                                                                                                                                                  |                                                                                                                                                   
aCGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTg                                                                                                                                                   >  2:564410/1‑149 (MQ=255)
  gcgTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGAt                                                                                                                                                 >  3:231021/1‑149 (MQ=255)
  gcgTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGAt                                                                                                                                                 >  3:305456/1‑149 (MQ=255)
          cAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCa                                                                                                                                         <  1:300686/149‑1 (MQ=255)
          cAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCa                                                                                                                                         <  3:4075/149‑1 (MQ=255)
           aGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCaa                                                                                                                                        <  2:191033/149‑1 (MQ=255)
                      aCGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGc                                                                                                                             <  4:361748/149‑1 (MQ=255)
                      aCGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGc                                                                                                                             <  4:286603/149‑1 (MQ=255)
                             aTAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATATTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTgc                                                                                                                      >  2:547535/1‑149 (MQ=255)
                              tAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTgcg                                                                                                                     >  1:196560/1‑149 (MQ=255)
                                    aGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATg                                                                                                               <  2:177174/149‑1 (MQ=255)
                                    aGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATg                                                                                                               <  4:10307/149‑1 (MQ=255)
                                         gACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGa                                                                                                                                                  >  2:368082/1‑109 (MQ=255)
                                         gACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGa                                                                                                                                                  <  1:368082/109‑1 (MQ=255)
                                         gACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCAt                                                                                                          >  4:254606/1‑149 (MQ=255)
                                            aaaCCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTc                                                                                                       <  1:139397/149‑1 (MQ=255)
                                              aCCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAAGATCATATGGCGGATTGCGCTCATGATCATTTCAt                                                                                                     >  4:404806/1‑149 (MQ=255)
                                               ccGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATc                                                                                                    <  4:306440/149‑1 (MQ=255)
                                               ccGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATc                                                                                                    >  4:314218/1‑149 (MQ=255)
                                                   gggTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGc                                                                                                >  3:10782/1‑149 (MQ=255)
                                                        gccgccgcTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATgg                                                                                           <  4:311239/149‑1 (MQ=255)
                                                           gccgcTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGa                                                                                        >  3:53811/1‑149 (MQ=255)
                                                           gccgcTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGa                                                                                        <  3:122623/149‑1 (MQ=255)
                                                           gccgcTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGa                                                                                        >  2:26567/1‑149 (MQ=255)
                                                              gcTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTc                                                                                     <  4:455053/149‑1 (MQ=255)
                                                                   tACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGaaa                                                                                <  2:576408/149‑1 (MQ=255)
                                                                      aTAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGcc                                                                             >  3:396212/1‑149 (MQ=255)
                                                                      aTAGTTCGGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGcc                                                                             >  2:317210/1‑149 (MQ=255)
                                                                        aGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAg                                                                           >  4:217604/1‑149 (MQ=255)
                                                                        aGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCAGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAg                                                                           >  4:478875/1‑149 (MQ=255)
                                                                          ttCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGc                                                                         <  3:314218/149‑1 (MQ=255)
                                                                                  aTACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCat                                                                                                                                  <  1:123137/84‑1 (MQ=255)
                                                                                  aTACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCat                                                                                                                                  >  2:123137/1‑84 (MQ=255)
                                                                                    acacTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAAc                                                               >  3:99388/1‑149 (MQ=255)
                                                                                    acacTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAAc                                                               <  2:326149/149‑1 (MQ=255)
                                                                                    acacTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAAc                                                               >  1:76163/1‑149 (MQ=255)
                                                                                    acacTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCACGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAAc                                                               <  4:536430/149‑1 (MQ=255)
                                                                                        ttCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCaa                                                           >  4:11180/1‑149 (MQ=255)
                                                                                           gATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCTTTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAAc                                                        >  1:48885/1‑149 (MQ=255)
                                                                                             tAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCg                                                      <  2:246517/149‑1 (MQ=255)
                                                                                                          cATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCa                                         <  2:562126/149‑1 (MQ=255)
                                                                                                           aTCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCAc                                        >  3:424469/1‑149 (MQ=255)
                                                                                                                        cTTCAGTTGTTGGCGTAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATATCATCCCGCCATGGTGAGTCCGAAATCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGa                           <  2:169984/149‑1 (MQ=255)
                                                                                                                        cTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCaaaa                                                         >  4:562479/1‑119 (MQ=255)
                                                                                                                        cTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCaaaa                                                         <  3:562479/119‑1 (MQ=255)
                                                                                                                         ttCAGTTGTTGGCGCAACCGACAGGTTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATc                                                                                                    >  1:64160/1‑75 (MQ=255)
                                                                                                                         ttCAGTTGTTGGCGCAACCGACAGGTTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATc                                                                                                    >  1:198385/1‑75 (MQ=255)
                                                                                                                         ttCAGTTGTTGGCGCAACCGACAGGTTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATc                                                                                                    <  2:198385/75‑1 (MQ=255)
                                                                                                                         ttCAGTTGTTGGCGCAACCGACAGGTTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATc                                                                                                    <  2:64160/75‑1 (MQ=255)
                                                                                                                         ttCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGAt                          >  4:92807/1‑149 (MQ=255)
                                                                                                                         ttCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGAt                          >  4:92407/1‑149 (MQ=255)
                                                                                                                         ttCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGAt                          >  2:530522/1‑149 (MQ=255)
                                                                                                                         ttCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAAGATCGTATGGCGGATTGCGCTCAGGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAGCCGGTAATCGGTTTCAGTTCGTAGAGCAGAt                          >  3:350243/1‑149 (MQ=255)
                                                                                                                              ttgttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGaaa                     <  1:501057/149‑1 (MQ=255)
                                                                                                                               tgttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAAc                    <  1:401205/149‑1 (MQ=255)
                                                                                                                                gttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACg                   >  3:199763/1‑149 (MQ=255)
                                                                                                                                   ggCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGccc                <  1:264603/149‑1 (MQ=255)
                                                                                                                                    gcgcAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCt               <  3:254606/149‑1 (MQ=255)
                                                                                                                                      gcAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGa             >  3:19348/1‑149 (MQ=255)
                                                                                                                                            gACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACggg       <  3:544995/149‑1 (MQ=255)
                                                                                                                                             aCAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGc      >  1:512459/1‑149 (MQ=255)
                                                                                                                                               aGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa    >  2:532157/1‑149 (MQ=255)
                                                                                                                                                gggTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaac   <  4:314329/149‑1 (MQ=255)
                                                                                                                                                 ggTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaacc  >  4:561727/1‑149 (MQ=255)
                                                                                                                                                 ggTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaacc  >  1:191264/1‑149 (MQ=255)
                                                                                                                                                 ggTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaacc  >  1:272293/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACC  >  NC_000913/1354759‑1355052

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