Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 553,495 C→T 6.0% R149H (CGC→CAC)  lpxH ← UDP‑2,3‑diacylglucosamine pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913553,4950CT 6.0% 223.6 / 3.1 67R149H (CGC→CAC) lpxHUDP‑2,3‑diacylglucosamine pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (33/30);  new base T (2/2);  total (35/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 = 9.11e-01

TGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCAACACCCGGCGACCATAAAGTTCGAGC  >  NC_000913/553357‑553635
                                                                                                                                          |                                                                                                                                            
tGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCacaaacaa                                                                                                                                    <  4:315719/149‑1 (MQ=255)
   ccGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACgg                                                                                                                                 >  1:201399/1‑149 (MQ=255)
     gTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACgg                                                                                                                                 >  1:352259/1‑147 (MQ=255)
     gTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACgg                                                                                                                                 <  2:352259/147‑1 (MQ=255)
         aTCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGc                                                                                                                           >  1:125319/1‑149 (MQ=255)
            aGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAg                                                                                                                        <  1:393589/149‑1 (MQ=255)
            aGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAg                                                                                                                        <  1:287734/149‑1 (MQ=255)
                aTTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAAt                                                                                                                    <  1:325584/149‑1 (MQ=255)
                aTTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGATGGCTTCTTTGCTGTTCGCCCGCATTCGCGCGGCAATGCGTTTGCACACAAACAACGGCAGGGCGAGGAAt                                                                                                                    <  4:476542/149‑1 (MQ=255)
                 ttGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATa                                                                                                                   <  1:338515/149‑1 (MQ=255)
                 ttGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATa                                                                                                                   <  4:548506/149‑1 (MQ=255)
                        tGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTCTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAGACAACGGCAGGGCGAGGAATAGCGTCTg                                                                                                            >  4:408265/1‑149 (MQ=255)
                                    aTCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGttt                                                                                                >  3:328998/1‑149 (MQ=255)
                                     tCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtg                                                                                               >  3:444636/1‑149 (MQ=255)
                                     tCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtg                                                                                               <  3:255568/149‑1 (MQ=255)
                                      cGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtgt                                                                                              <  1:554006/149‑1 (MQ=255)
                                       gCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtgtg                                                                                             <  4:328998/149‑1 (MQ=255)
                                          cTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGAc                                                                                          <  3:324607/149‑1 (MQ=255)
                                               caccGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTgg                                                                                     >  4:128296/1‑149 (MQ=255)
                                                accGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGgc                                                                                    >  1:356903/1‑149 (MQ=255)
                                                 ccGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGgcg                                                                                   >  1:123424/1‑149 (MQ=255)
                                                        ttGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAg                                                                            <  1:156114/149‑1 (MQ=255)
                                                         tGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGc                                                                           <  4:21603/149‑1 (MQ=255)
                                                         tGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGc                                                                           <  2:572436/149‑1 (MQ=255)
                                                          ggTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGATGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGcc                                                                          <  3:233425/149‑1 (MQ=255)
                                                           gTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAGACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCt                                                                         <  3:408265/149‑1 (MQ=255)
                                                           gTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCt                                                                         >  1:116815/1‑149 (MQ=255)
                                                            ttAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTTGCGCGCATTCACGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGCCGCGAAAAGCCTg                                                                        >  3:536186/1‑149 (MQ=255)
                                                            ttAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGcccg                                                                        >  3:498465/1‑147 (MQ=255)
                                                            ttAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTg                                                                        >  1:64473/1‑149 (MQ=255)
                                                            ttAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTg                                                                        >  1:335863/1‑149 (MQ=255)
                                                             tAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGa                                                                       <  1:64013/149‑1 (MQ=255)
                                                             tAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGa                                                                       <  3:249648/149‑1 (MQ=255)
                                                                cGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATaa                                                                    <  1:193798/149‑1 (MQ=255)
                                                                 gTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGTGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAAc                                                                   >  1:525125/1‑149 (MQ=255)
                                                                 gTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGTGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAAc                                                                   >  1:86078/1‑149 (MQ=255)
                                                                 gTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAAc                                                                   >  1:354446/1‑149 (MQ=255)
                                                                 gTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAAc                                                                   >  1:143333/1‑149 (MQ=255)
                                                                        aTCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTc                                                            >  4:281567/1‑149 (MQ=255)
                                                                            ccAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCg                                                        >  4:113004/1‑149 (MQ=255)
                                                                              aGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGt                                                      >  4:362171/1‑149 (MQ=255)
                                                                               gCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTg                                                     <  2:125319/149‑1 (MQ=255)
                                                                                   tttACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGcaca                                                 >  1:330374/1‑149 (MQ=255)
                                                                                               ggCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCAt                                     <  2:463138/149‑1 (MQ=255)
                                                                                               ggCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCAt                                     <  3:311125/149‑1 (MQ=255)
                                                                                               ggCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCAt                                     <  1:510778/149‑1 (MQ=255)
                                                                                                  ttcttTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCa                                  >  1:133361/1‑149 (MQ=255)
                                                                                                  ttcttTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCa                                  >  3:529827/1‑149 (MQ=255)
                                                                                                  ttcttTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCa                                  >  4:279942/1‑149 (MQ=255)
                                                                                                  ttcttTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCa                                  >  4:3493/1‑149 (MQ=255)
                                                                                                     tttGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATaa                               >  4:89818/1‑149 (MQ=255)
                                                                                                     tttGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATaa                               >  1:77382/1‑149 (MQ=255)
                                                                                                        gCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCa                            >  4:355410/1‑149 (MQ=255)
                                                                                                          tGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCAac                          <  3:128296/149‑1 (MQ=255)
                                                                                                              cgcgcgCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCAACAccc                      >  2:260436/1‑149 (MQ=255)
                                                                                                                   gcATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCAACACCCGGCGa                 <  2:356903/149‑1 (MQ=255)
                                                                                                                      ttCGCGCGGCAATGCGTTTGTGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCAACACCCGGCGACCa              <  2:86078/149‑1 (MQ=255)
                                                                                                                      ttCGCGCGGCAATGCGTTTGTGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCAACACCCGGCGACCa              <  2:525125/149‑1 (MQ=255)
                                                                                                                         gcgcgGCAATGCGTTTGCGCACAAGCAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCAACACCCGGCGACCATaa           >  1:497236/1‑149 (MQ=255)
                                                                                                                          cgcgGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCAACACCCGGCGACCATaaa          <  2:200687/149‑1 (MQ=255)
                                                                                                                           gcgGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCAACACCCGGCGACCATAAAg         <  1:558695/149‑1 (MQ=255)
                                                                                                                              gCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCAACACCCGGCGACCATAAAGTTc      <  1:295617/149‑1 (MQ=255)
                                                                                                                                aaTGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCAACACCCGGCGACCATAAAGTTCGa    >  2:387376/1‑149 (MQ=255)
                                                                                                                                aaTGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCAACACCCGGCGACCATAAAGTTCGa    >  4:464601/1‑149 (MQ=255)
                                                                                                                                 aTGCGTTTCCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTGTGTGGACCTTTGCGCGAAAAGCCTTTTAACCCGCATCATCGGTGCACAGCGTGTCGCCATGCATAATCAACACCCGGCGACCATAAAGGTCGAg   <  4:505839/149‑1 (MQ=255)
                                                                                                                                  tGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCAACACCCGGCGACCATAAAGTTCGAGc  >  3:470437/1‑149 (MQ=255)
                                                                                                                                  tGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCAACACCCGGCGACCATAAAGTTCGAGc  >  3:21125/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                            
TGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCATGCATAATCAACACCCGGCGACCATAAAGTTCGAGC  >  NC_000913/553357‑553635

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: ‑