Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,224,508 T→C 6.2% intergenic (‑620/+201) PP_3703 ← / ← PP_3704 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,224,5080TC 6.2% 188.6 / 2.1 64intergenic (‑620/+201)PP_3703/PP_3704hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base C (2/2);  total (30/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 = 1.00e+00

TCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGGATATGC  >  NC_002947/4224362‑4224653
                                                                                                                                                    |                                                                                                                                                 
tCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCAt                                                                                                                                                   <  7:269979/147‑1 (MQ=255)
   aTCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTg                                                                                                                                              >  2:349752/1‑149 (MQ=255)
     ccAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGt                                                                                                                                            >  2:330052/1‑149 (MQ=255)
       aaaCGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCAACCAGCAATGCACAGGTGCACTGCTggg                                                                                                                                            <  7:44630/149‑2 (MQ=255)
         aCGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGcc                                                                                                                                        <  5:218880/149‑1 (MQ=255)
              aaTGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACgt                                                                                                                                    <  2:308573/148‑1 (MQ=255)
               aTGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGTGAATACCCCGAGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACgtgt                                                                                                                                  <  5:258926/149‑1 (MQ=255)
                tGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACgtgt                                                                                                                                  >  8:276664/1‑148 (MQ=255)
                    ccTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTaaa                                                                                                                               <  8:49978/147‑1 (MQ=255)
                    ccTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTaaa                                                                                                                               >  7:49978/1‑147 (MQ=255)
                    ccTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGA‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAgg                                                                                                                             >  5:177506/1‑149 (MQ=255)
                     cTAAGCCTGGCGGTGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGc                                                                                                                            <  3:174753/149‑1 (MQ=255)
                     cTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGc                                                                                                                            >  6:350996/1‑149 (MQ=255)
                     cTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGc                                                                                                                            >  5:163279/1‑149 (MQ=255)
                      tAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAATGACGGCCAGATACGTGATAAGCGGAGGGTTGCGGGAAATGCCCGCGGAAAATAATCGTAGGTGGC‑‑CGAGCAGTGCTCCTGTGCATTGCTg                                                                                                                                              <  5:348967/130‑1 (MQ=255)
                      tAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCt                                                                                                                           >  6:38833/1‑149 (MQ=255)
                         gCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTAt                                                                                                                        <  5:38833/149‑1 (MQ=255)
                         gCCTGGCGGCGGGGAACAACCTCACTCCCGTAACAAAGACGGCCAGATACATGATAAGCGGAGCGATGCGGGAAAACCCCGAGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTAt                                                                                                                        <  2:314283/149‑1 (MQ=255)
                                ggcggGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTaa                                                                                                                 <  2:67305/149‑1 (MQ=255)
                                 gcggGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTaaa                                                                                                                <  5:350996/149‑1 (MQ=255)
                                       aacaacCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAACCGGAGCGATGCTGGAAAACCCCGAGGAAAATAAACCTAGGTGGC‑‑ACAGCAGAGCACCTGAGCATTCCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAg                                                                                                           >  3:100878/1‑148 (MQ=255)
                                           acCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTg                                                                                                      <  6:159115/149‑1 (MQ=255)
                                            ccTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTgg                                                                                                     >  8:190469/1‑149 (MQ=255)
                                            ccTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTgg                                                                                                     <  2:338511/149‑1 (MQ=255)
                                            ccTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAATGCACAGGTGCACTGCTgg                                                                                                                                             >  3:293588/1‑109 (MQ=255)
                                            ccTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAATGCACAGGTGCACTGCTgg                                                                                                                                             <  4:293588/109‑1 (MQ=255)
                                             cTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTggg                                                                                                    >  3:271682/1‑149 (MQ=255)
                                             cTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTggg                                                                                                    <  8:140106/149‑1 (MQ=255)
                                               cGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGc                                                                                                                                         <  1:182062/110‑1 (MQ=255)
                                               cGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGc                                                                                                                                         >  2:182062/1‑110 (MQ=255)
                                                  tCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAg                                                                                                                     >  4:276361/1‑127 (MQ=255)
                                                  tCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAg                                                                                                                     <  3:276361/127‑1 (MQ=255)
                                                    ccGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATcgcg                                                                                             <  8:177701/149‑1 (MQ=255)
                                                      gTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGcc                                                                                           <  1:41290/149‑1 (MQ=255)
                                                           aaaGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCt                                                                                      >  1:281247/1‑149 (MQ=255)
                                                              gACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTg                                                                                     >  3:169565/1‑147 (MQ=255)
                                                                       aTACGTGATAAGCTGAGCTATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCGCATCTGTAGCGGAgat                                                                           <  4:100878/148‑1 (MQ=255)
                                                                           gTGATAAGCTGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATAc                                                                      <  1:349752/149‑1 (MQ=255)
                                                                           gTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAgatgat                                                                        <  7:226177/147‑1 (MQ=255)
                                                                           gTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAgatgat                                                                        >  8:226177/1‑147 (MQ=255)
                                                                                 aGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAAcc                                                                >  3:177130/1‑149 (MQ=255)
                                                                                    ggAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCt                                                               >  7:255913/1‑147 (MQ=255)
                                                                                      aGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAgg                                                           <  2:281247/149‑1 (MQ=255)
                                                                                             cGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTc                                                     <  7:190469/148‑1 (MQ=255)
                                                                                              gggAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTGAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCa                                                   >  5:195612/1‑149 (MQ=255)
                                                                                               ggAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCAt                                                  >  7:111286/1‑149 (MQ=255)
                                                                                                  aaaCCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATa                                                                       <  5:112727/125‑1 (MQ=255)
                                                                                                  aaaCCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATa                                                                       >  6:112727/1‑125 (MQ=255)
                                                                                                   aaCCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTaa                                              <  7:276664/149‑1 (MQ=255)
                                                                                                         gcgGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTgg                                        <  8:111286/149‑1 (MQ=255)
                                                                                                               aaTAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAgc                                  <  8:255913/149‑1 (MQ=255)
                                                                                                                aTAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTa                                                             >  2:281018/1‑121 (MQ=255)
                                                                                                                aTAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTa                                                             <  1:281018/121‑1 (MQ=255)
                                                                                                                 tAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGt                                                      >  7:282789/1‑127 (MQ=255)
                                                                                                                 tAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGt                                                      <  8:282789/127‑1 (MQ=255)
                                                                                                                   aaCCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAgcgtc                               <  8:58606/148‑3 (MQ=255)
                                                                                                                       gttctGGC‑‑CCATCATTGCACCTGTGCATTTCTAGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAGCGCCaaggaa                         <  2:333055/145‑1 (MQ=255)
                                                                                                                             gC‑‑CCAGCAGTGCACCTGTGCATTGCTAGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAGCGCCaaggaa                         >  1:333055/1‑143 (MQ=255)
                                                                                                                              c‑‑ccAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTCGGTATCGCGCCCATCTGTAGCGGACATGATACAGAACCTTATGATTGTCCAGTTAAACATGGAAAAGCCCCAAGGAAGGTTCt                   >  8:244278/1‑149 (MQ=255)
                                                                                                                                 ccAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCtg                  >  5:252310/1‑149 (MQ=255)
                                                                                                                                 ccAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAACACAAGACTTGGGATCGCGCCCATCTGTATCGGAGATGATACAGAAACTTATGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAATGTtgtg                  >  1:101735/1‑149 (MQ=255)
                                                                                                                                 ccAGCAATGCACAGGTGCACTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCtg                  >  5:136261/1‑149 (MQ=255)
                                                                                                                                   agcagTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCtgtg                <  1:330052/149‑1 (MQ=255)
                                                                                                                                             cTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGGa       >  2:315947/1‑148 (MQ=255)
                                                                                                                                                  cATTGCTGGTTGCCAAGTGTATAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCGATCTGTATCGGAGTTGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGGATATGc  <  7:244278/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGC‑‑CCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGGATATGC  >  NC_002947/4224362‑4224653

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: