Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,059,505 A→T 7.3% N443K (AAT→AAA PP_4458 ← opine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,059,5050AT 7.3% 186.2 / 5.3 69N443K (AAT→AAAPP_4458opine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base T (3/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

CCTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCGCGT  >  NC_002947/5059360‑5059647
                                                                                                                                                 |                                                                                                                                               
ccTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTc                                                                                                                                              <  7:171706/149‑1 (MQ=255)
ccTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTc                                                                                                                                              <  5:131352/149‑1 (MQ=255)
        aGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTg                                                                                                                                      <  4:269689/148‑1 (MQ=255)
          aCCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAg                                                                                                                                   >  7:245834/1‑149 (MQ=255)
           ccAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCAACAAATCGAGTTGAGTCTTCGGAACAGCCACGACCCAGGCTATCCGGATGGCATCAACAGCTACTGATCAAAGAACGAAATTCGACATTGCAGc                                                                                                                                  >  8:247406/1‑149 (MQ=255)
            caaacaaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACAtt                                                                                                                                       >  8:342159/1‑143 (MQ=255)
            caaacaaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACAtt                                                                                                                                       <  7:342159/143‑1 (MQ=255)
             aaacaaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCgg                                                                                                                                <  5:310879/149‑1 (MQ=255)
              aacaaaGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCggg                                                                                                                               >  5:45834/1‑149 (MQ=255)
                     ggATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCGAATTTCGGACATTGctc                                                                                                                                   >  7:255711/1‑136 (MQ=255)
                     ggATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCGAATTTCGGACATTGctc                                                                                                                                   <  8:255711/138‑3 (MQ=255)
                          gcgcATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATgg                                                                                                                    >  1:95574/1‑148 (MQ=255)
                            gcATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAg                                                                                                                 <  5:197382/149‑1 (MQ=255)
                                aCGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATAATCAGCTACTTAGCAATGTACGAAATTCTAAATTGCAGCGGGTGCGCCGATGGCTGAAGc                                                                                                             >  8:164555/1‑149 (MQ=255)
                                  gtcgtcgtGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGcc                                                                                                            >  2:189324/1‑148 (MQ=255)
                                      tcgtGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATc                                                                                                       <  7:248865/149‑1 (MQ=255)
                                        gtGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCg                                                                                                      >  4:169265/1‑148 (MQ=255)
                                          gAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGACGAATCGTa                                                                                                    >  1:212364/1‑148 (MQ=255)
                                           aaTCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTaaa                                                                                                  <  8:194324/149‑1 (MQ=255)
                                                 aaCCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAAct                                                                                               <  6:45834/146‑1 (MQ=255)
                                                 aaCCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTg                                                                                            <  8:245834/149‑1 (MQ=255)
                                                  aCCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGt                                                                                           <  1:189324/149‑1 (MQ=255)
                                                  aCCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGt                                                                                           <  7:247406/149‑1 (MQ=255)
                                                    cTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCa                                                                                         <  3:311285/149‑1 (MQ=255)
                                                      gcgcCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGc                                                                                       >  8:32436/1‑149 (MQ=255)
                                                          cTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGaa                                                                                                         >  3:177083/1‑127 (MQ=255)
                                                          cTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGaa                                                                                                         <  4:177083/127‑1 (MQ=255)
                                                          cTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCag                                                                                    >  8:62632/1‑148 (MQ=255)
                                                                 cacaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCgg                                                                             <  2:95574/148‑1 (MQ=255)
                                                                  acaAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGgcg                                                                           >  2:75468/1‑149 (MQ=255)
                                                                       tCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCg                                                                       >  5:191193/1‑148 (MQ=255)
                                                                       tCGAGTAGAGTCTGCGGATCAGCCTCGCCCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGAGATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCCGCCAGTAGTCGGcgcgcg                                                                       >  2:230107/1‑148 (MQ=255)
                                                                            ttGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGCCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATc                                                                 >  3:160441/1‑149 (MQ=255)
                                                                                   tGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGagcag                                                           >  2:119456/1‑148 (MQ=255)
                                                                                            gCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCCATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTaaaa                                                 >  1:114697/1‑149 (MQ=255)
                                                                                             ccTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTaaaa                                                 >  3:193604/1‑148 (MQ=255)
                                                                                                cGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGt                                                                                                     >  3:191787/1‑93 (MQ=255)
                                                                                                cGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGt                                                                                                     <  4:191787/93‑1 (MQ=255)
                                                                                                      tGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATcgcg                                       >  3:120379/1‑149 (MQ=255)
                                                                                                               ggATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCg                                                                                                      >  3:338122/1‑77 (MQ=255)
                                                                                                               ggATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCg                                                                                                      <  4:338122/77‑1 (MQ=255)
                                                                                                                gATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCtt                             >  5:101869/1‑149 (MQ=255)
                                                                                                                gATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCAGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCtt                             <  4:193604/149‑1 (MQ=255)
                                                                                                                     catcatCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGagc                                                             <  6:68221/112‑1 (MQ=255)
                                                                                                                     catcatCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGagc                                                             >  5:68221/1‑112 (MQ=255)
                                                                                                                     catcatCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGt                                >  3:115109/1‑141 (MQ=255)
                                                                                                                     catcatCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGt                                <  4:115109/141‑1 (MQ=255)
                                                                                                                      atcatcAGCTACTGAGCAATGTCAGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTTTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGGCCCGTAAAACACGATCGCGATCTGGTCTTGACTGc                       >  4:220841/1‑149 (MQ=255)
                                                                                                                        catcaGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAACACGATCGCGCTCTGGTCTTGACGGCttt                    >  1:328239/1‑149 (MQ=255)
                                                                                                                          tcaGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCtttt                   >  6:263998/1‑149 (MQ=255)
                                                                                                                          tcaGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCtttt                   <  4:120379/149‑1 (MQ=255)
                                                                                                                           caGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTg                  >  5:267490/1‑149 (MQ=255)
                                                                                                                            aGCTACTGAGCAATGTGCGAAATTCTACATTGCAGCGTGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGGAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTgg                 <  2:212364/149‑1 (MQ=255)
                                                                                                                            aGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGTGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTg                  <  1:75468/148‑1 (MQ=255)
                                                                                                                            aGCTACGGAGCAATGTCCGACATTCGCCATTGCAGCGGGTGCGCCGATGGCAGACGCCGAATCGTAAGACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAACACACGATCGCGCTCTGGTCTTGACGGCTTTTg                  <  4:247276/148‑1 (MQ=255)
                                                                                                                              cTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGcc               <  1:119456/149‑1 (MQ=255)
                                                                                                                                 cTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCt              <  6:191193/147‑1 (MQ=255)
                                                                                                                                 cTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTgg            <  3:169265/149‑1 (MQ=255)
                                                                                                                                 cTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTgg            >  5:123689/1‑149 (MQ=255)
                                                                                                                                  gctgCAATGTCGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCTGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTgg            >  4:106336/4‑148 (MQ=255)
                                                                                                                                  gctgCAATGTCGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGag                                                              >  2:237717/4‑98 (MQ=38)
                                                                                                                                  gctgCAATGTCGAATTTCGGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGag                                                              <  1:237717/95‑1 (MQ=38)
                                                                                                                                     gCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAg        <  7:32436/149‑1 (MQ=255)
                                                                                                                                       aaTGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTc      <  7:164555/149‑1 (MQ=255)
                                                                                                                                         tGTCCGAAATTCGACATTGCAGCTGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTcgc    <  2:120144/149‑1 (MQ=255)
                                                                                                                                           tCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTcgcg   <  2:114697/148‑1 (MQ=255)
                                                                                                                                            ccGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCGCGt  >  1:31327/1‑148 (MQ=255)
                                                                                                                                             cGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGag                                                              >  7:278664/1‑87 (MQ=255)
                                                                                                                                             cGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGag                                                              <  8:278664/87‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CCTGTTTCAGACCAAACAAAGGGATGGCGCATACGTCGTCGTGAATCTGAACCTGCGCCTTTTTCCACAAATCGAGTTGAGTCTGCGGATCAGCCTCGACCCTGGCTTTCCGGATGGCATCATCAGCTACTGAGCAATGTCCGAAATTCGACATTGCAGCGGGTGCGCCGATGGCAGAAGCCGAATCGTAAAACTCTGTCAGCCAGTAGTCGGCGCCCGGATATCGAGCAGCACCGTAAAACACGATCGCGCTCTGGTCTTGACGGCTTTTGGCCTGGTAAGTCGCGT  >  NC_002947/5059360‑5059647

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