Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,758,424 C→G 5.1% V454L (GTG→CTG)  pvdP ← pyoverdine biosynthesis‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,758,4240CG 5.1% 260.9 / 3.2 79V454L (GTG→CTG) pvdPpyoverdine biosynthesis‑like protein
Reads supporting (aligned to +/- strand):  ref base C (36/39);  new base G (2/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

GTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAGCGCGCC  >  NC_002947/4758279‑4758572
                                                                                                                                                 |                                                                                                                                                    
gTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGa                                                                                                                                                    <  4:244935/148‑1 (MQ=255)
    ccGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTgc                                                                                                                                                <  5:14695/148‑1 (MQ=255)
      ggATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTgcgc                                                                                                                                              >  7:334201/1‑148 (MQ=255)
           cATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCAGTCGTGCGCCAgga                                                                                                                                         <  4:125341/148‑2 (MQ=255)
            aTCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTg                                                                                                                                       <  1:301435/149‑1 (MQ=255)
            aTCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTg                                                                                                                                       >  7:220001/1‑149 (MQ=255)
            aTCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTg                                                                                                                                       >  5:84292/1‑149 (MQ=255)
                ccACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAAt                                                                                                                                   >  3:145563/1‑149 (MQ=255)
                 cACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATc                                                                                                                                  <  5:364370/149‑1 (MQ=255)
                 cACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATc                                                                                                                                  <  8:117000/149‑1 (MQ=255)
                 cACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATc                                                                                                                                  >  7:15368/1‑149 (MQ=255)
                 cACGCTACGCCCTGGGCACAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATc                                                                                                                                  >  5:380967/1‑149 (MQ=255)
                  aCGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCg                                                                                                                                 >  6:155441/1‑149 (MQ=255)
                   cGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGt                                                                                                                                <  8:334201/149‑1 (MQ=255)
                           ccTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGgc                                                                                                                        >  2:333002/1‑149 (MQ=255)
                            cTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGgcg                                                                                                                       <  2:197162/149‑1 (MQ=255)
                                 cgcAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGaa                                                                                                                  >  5:124709/1‑149 (MQ=255)
                                   cAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAAcc                                                                                                                <  1:230718/149‑1 (MQ=255)
                                    aGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCa                                                                                                               <  8:318801/149‑1 (MQ=255)
                                    aGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCa                                                                                                               <  5:186877/149‑1 (MQ=255)
                                    aGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCa                                                                                                               <  5:35852/149‑1 (MQ=255)
                                      ccccGCCACGTTACTGCAGCCGAGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATAGATCAAGACCCCGCAACCAGt                                                                                                             >  6:262936/1‑149 (MQ=255)
                                       cccGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTc                                                                                                            >  2:63774/1‑149 (MQ=255)
                                        ccGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCCCCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCATGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGAcc                                                                                                           >  1:230366/1‑149 (MQ=255)
                                        ccGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTc                                                                                                            >  8:311131/1‑148 (MQ=255)
                                         cGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCACGGTGGAATCGTTCATGGGCGCGGAACCAGTCCt                                                                                                          >  6:113532/1‑149 (MQ=255)
                                            cACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGa                                                                                                       <  1:45118/149‑1 (MQ=255)
                                               gTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACTATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCgg                                                                                                                    <  3:264683/133‑1 (MQ=255)
                                               gTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGacac                                                                                                    <  1:47523/149‑1 (MQ=255)
                                               gTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACAATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCgg                                                                                                                    >  4:264683/1‑133 (MQ=255)
                                                ttACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACg                                                                                                   <  5:184889/149‑1 (MQ=255)
                                                    tGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcg                                                                                               >  1:181300/1‑149 (MQ=255)
                                                       aGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgtc                                                                                             >  4:81084/1‑148 (MQ=255)
                                                       aGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTCCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGtcgtcg                                                                                            >  1:247450/1‑149 (MQ=255)
                                                          cgtgcgtgTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATc                                                                                         <  3:81084/149‑1 (MQ=255)
                                                          cgtgcgtgTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATc                                                                                         >  7:34515/1‑149 (MQ=255)
                                                           gtgcgtgTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATcc                                                                                        >  8:142254/1‑149 (MQ=255)
                                                              cgtgTCCGGGCCCAGCCAAGGATCGGCGTCCTCCCCCCAGCGGCCGGGCGCGAACCATGGCACGCCATTGACCTCGTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGc                                                                                     <  8:246264/149‑1 (MQ=255)
                                                              cgtgTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGc                                                                                     <  5:180455/149‑1 (MQ=255)
                                                               gtgtCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGc                                                                                     >  6:164161/1‑148 (MQ=255)
                                                                 gtCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCAt                                                                                  >  4:343903/1‑149 (MQ=255)
                                                                  tCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCAt                                                                                  <  7:311131/148‑1 (MQ=255)
                                                                    cGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATgg                                                                                >  1:191739/1‑148 (MQ=255)
                                                                     ggCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGaa                                                                              <  5:164161/149‑1 (MQ=255)
                                                                     ggCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGaa                                                                              <  2:181300/149‑1 (MQ=255)
                                                                       ccccAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAgc                                                                            <  5:155441/149‑1 (MQ=255)
                                                                       ccccAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAgc                                                                            <  6:31000/149‑1 (MQ=255)
                                                                         ccAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAgcgc                                                                          <  8:220001/149‑1 (MQ=255)
                                                                              cAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGaa                                                                     <  6:380967/149‑1 (MQ=255)
                                                                                       gCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGtt                                                            >  5:143179/1‑149 (MQ=255)
                                                                                       gCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGtt                                                            >  3:193374/1‑149 (MQ=255)
                                                                                                 cccAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGaa                                                  <  7:142254/149‑1 (MQ=255)
                                                                                                  ccAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAAcaca                                                                 <  6:39434/133‑1 (MQ=255)
                                                                                                  ccAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAAcaca                                                                 >  5:39434/1‑133 (MQ=255)
                                                                                                  ccAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTggg                                                    >  6:64750/1‑146 (MQ=255)
                                                                                                   cAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAaga                                                <  8:34515/149‑1 (MQ=255)
                                                                                                       gcCCCGGCGCGAACCATGGCACGGCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGTGCGCGGAACCAGTCCTCTACACGTTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAgagatc                                            <  5:262936/149‑3 (MQ=255)
                                                                                                       gcCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAAc                                            <  7:102251/149‑1 (MQ=255)
                                                                                                       gcCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAAc                                            <  2:191739/149‑1 (MQ=255)
                                                                                                           cGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGt                                        <  1:333002/149‑1 (MQ=255)
                                                                                                           cGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGt                                        >  1:45374/1‑149 (MQ=255)
                                                                                                              cgcgAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTccc                                     <  4:106990/149‑1 (MQ=255)
                                                                                                               gcgAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGTAAGAGAAc                                            <  1:56331/141‑1 (MQ=255)
                                                                                                               gcgAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAAc                                            >  2:56331/1‑141 (MQ=255)
                                                                                                                cgAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCg                                   <  6:124709/149‑1 (MQ=255)
                                                                                                                        ggCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTc                           >  6:88730/1‑149 (MQ=255)
                                                                                                                           aCGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGtt                        >  5:329602/1‑149 (MQ=255)
                                                                                                                           aCGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGt                         >  6:380438/1‑148 (MQ=255)
                                                                                                                                aTTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgcg                   >  6:300769/1‑149 (MQ=255)
                                                                                                                                    aCCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgcgcgc                <  5:143917/148‑1 (MQ=255)
                                                                                                                                     cctcctGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCa              <  6:84292/149‑1 (MQ=255)
                                                                                                                                       ccctGGCGCAGTCGTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATa            >  1:128707/2‑149 (MQ=255)
                                                                                                                                        cctGGCGCAGTCGTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgcgcg                 <  3:87072/143‑1 (MQ=255)
                                                                                                                                        cctGGCGCAGTCGTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgcgcg                 >  4:87072/1‑143 (MQ=255)
                                                                                                                                         ctGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATAc           >  1:260820/1‑148 (MQ=255)
                                                                                                                                         ctGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAACCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATAc           >  5:45072/1‑148 (MQ=255)
                                                                                                                                            gcgcACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCa         <  5:113532/147‑1 (MQ=255)
                                                                                                                                              gcACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAgcgc     <  5:64750/149‑1 (MQ=255)
                                                                                                                                                aCGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACCCAGGGTTGACGTGGGAAGAGAACGGGTACCCGAGAAAGTCGTTCTGCGCGCCATACCAgcgcgc   >  6:196564/1‑149 (MQ=255)
                                                                                                                                                 cGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAGCGcgcc  <  5:380438/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAGCGCGCC  >  NC_002947/4758279‑4758572

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