Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,758,425 G→T 7.9% V453V (GTC→GTA) ‡ pvdP ← pyoverdine biosynthesis‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,758,4250GT 7.9% 227.8 / 3.5 71V453V (GTC→GTA) ‡pvdPpyoverdine biosynthesis‑like protein
Reads supporting (aligned to +/- strand):  ref base G (37/28);  new base T (3/3);  total (40/31)
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

TTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAGCGCGC  >  NC_002947/4758280‑4758571
                                                                                                                                                  |                                                                                                                                                  
tttCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACt                                                                                                                                                 >  7:192117/1‑149 (MQ=255)
   ccGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTgcg                                                                                                                                              >  3:74498/1‑149 (MQ=255)
    cGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTgcgc                                                                                                                                             >  2:68154/1‑149 (MQ=255)
        tCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCAGTCGTGCGCCAgg                                                                                                                                         <  6:71221/149‑1 (MQ=255)
           aTCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTg                                                                                                                                      <  8:61617/149‑1 (MQ=255)
              tcCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCAGTCGTGCGCCAggaggtc                                                                                                                                    >  6:351445/1‑143 (MQ=255)
                          ccTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCCCCTGGCGCAGTCGTGCGCCAggaggt                                                                                                                                     >  4:124014/1‑131 (MQ=255)
                          ccTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCCCCTGGCGCAGTCGTGCGCCAggaggt                                                                                                                                     <  3:124014/135‑5 (MQ=255)
                           cTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGgcg                                                                                                                      >  4:223624/1‑149 (MQ=255)
                              ggCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGgcgc                                                                                                                     <  3:95900/147‑1 (MQ=255)
                              ggCGCAGACCCGCCACGGTACTGCAGCCGTGCGTGTCCCGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGgcgc                                                                                                                     <  8:332229/147‑1 (MQ=255)
                                 gcAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAAc                                                                                                                <  8:192117/149‑1 (MQ=255)
                                     ccccGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGt                                                                                                            >  5:243923/1‑149 (MQ=255)
                                       ccGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGTCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTcc                                                                                                          >  2:117812/1‑149 (MQ=255)
                                       ccGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTcc                                                                                                          >  1:314391/1‑149 (MQ=255)
                                           cACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGa                                                                                                      >  5:313664/1‑149 (MQ=255)
                                               ttACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACg                                                                                                  >  5:242464/1‑149 (MQ=255)
                                                     cAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGGGCGGAACCAGTCCTCGGCACGgtcgtc                                                                                            >  3:236447/1‑149 (MQ=255)
                                                     cAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgtc                                                                                            <  3:2705/149‑1 (MQ=255)
                                                     cAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCCCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACCGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgtc                                                                                            >  4:144081/1‑149 (MQ=255)
                                                          gtgcgtgTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATc                                                                                        >  5:177977/1‑148 (MQ=255)
                                                           tgcgtgGCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACCGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCa                                                                                      <  3:144081/149‑1 (MQ=255)
                                                              gtgtCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGcc                                                                                   <  6:243923/149‑1 (MQ=255)
                                                                       cccAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgt                                                                                             >  5:269565/1‑130 (MQ=255)
                                                                       cccAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgt                                                                                             <  6:269565/130‑1 (MQ=255)
                                                                       cccAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAgcg                                                                          >  5:80374/1‑149 (MQ=255)
                                                                        ccAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAgcgc                                                                         <  4:27303/149‑1 (MQ=255)
                                                                              aaGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATgg                                                                                                                         >  7:197992/1‑95 (MQ=255)
                                                                              aaGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATgg                                                                                                                         <  8:197992/95‑1 (MQ=255)
                                                                              aaGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAg                                                                            <  4:102921/140‑1 (MQ=255)
                                                                              aaGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAg                                                                            >  3:102921/1‑140 (MQ=255)
                                                                                  aTCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAg                                                               <  2:314391/149‑1 (MQ=255)
                                                                                   tCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAgg                                                              <  4:74498/149‑1 (MQ=255)
                                                                                      tcGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGCTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGtt                                                           <  1:117812/148‑1 (MQ=255)
                                                                                      gCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAggg                                                             >  7:310937/1‑147 (MQ=255)
                                                                                      gCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAggg                                                             <  8:310937/147‑1 (MQ=255)
                                                                                        gACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGTTGAc                                                        <  1:68154/149‑1 (MQ=255)
                                                                                           cTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAg                                                                                                             <  3:97544/94‑1 (MQ=255)
                                                                                           cTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAg                                                                                                             >  4:97544/1‑94 (MQ=255)
                                                                                           cTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGtt                                                           <  6:315762/144‑1 (MQ=255)
                                                                                           cTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGtt                                                           >  5:315762/1‑144 (MQ=255)
                                                                                           cTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGt                                                      >  8:271534/1‑149 (MQ=255)
                                                                                           cTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGt                                                      >  3:349667/1‑149 (MQ=255)
                                                                                            tCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTg                                                     >  1:173487/1‑149 (MQ=255)
                                                                                             ccACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTgg                                                    >  3:337451/1‑149 (MQ=255)
                                                                                              cACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTgg                                                    >  3:181768/1‑148 (MQ=255)
                                                                                              cACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTgg                                                    <  3:223624/148‑1 (MQ=255)
                                                                                               aCCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGAGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAg                                                               >  1:261661/1‑137 (MQ=255)
                                                                                               aCCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAg                                                               <  2:261661/136‑1 (MQ=255)
                                                                                               aCCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTggg                                                   >  4:73696/1‑148 (MQ=255)
                                                                                                cccAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGaa                                                                    >  8:100726/1‑130 (MQ=255)
                                                                                                cccAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGaa                                                                    <  7:100726/130‑1 (MQ=255)
                                                                                                  cAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGaa                                                 >  3:340703/1‑147 (MQ=255)
                                                                                                   aGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGaa                                                 >  2:204903/1‑146 (MQ=255)
                                                                                                   aGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAagag                                              <  6:80374/149‑1 (MQ=255)
                                                                                                      gcCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAAc                                           <  2:173487/149‑1 (MQ=255)
                                                                                                        cccGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACgg                                         <  4:349667/149‑1 (MQ=255)
                                                                                                              gcgAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAAAGAACGGGTcccc                                   >  2:337409/1‑149 (MQ=255)
                                                                                                                    cATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCgaga                               >  7:331363/1‑147 (MQ=255)
                                                                                                                           cGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTc                      >  2:199436/1‑149 (MQ=255)
                                                                                                                            gCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCt                     <  4:340703/149‑1 (MQ=255)
                                                                                                                                   aCCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgcgcgc               <  6:313664/148‑1 (MQ=255)
                                                                                                                                   aCCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgcgcgc               >  2:249284/1‑148 (MQ=255)
                                                                                                                                      tcctGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATa           >  3:241945/1‑149 (MQ=255)
                                                                                                                                      tcctGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCAAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCAt            >  8:254040/1‑148 (MQ=255)
                                                                                                                                      ccctGGCGCAGTCGTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCAt            >  3:219662/2‑148 (MQ=255)
                                                                                                                                      ccctGGCGCAGTCGTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCAt            <  4:219662/147‑1 (MQ=255)
                                                                                                                                            gcgcACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAg       <  4:337451/148‑1 (MQ=255)
                                                                                                                                              gcACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAgcgc    >  7:14679/1‑149 (MQ=255)
                                                                                                                                              gcACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAgcgc    >  6:62868/1‑149 (MQ=255)
                                                                                                                                                aCGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAgcgcgc  <  7:271534/149‑1 (MQ=255)
                                                                                                                                                  gACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTc                                                                                                           >  7:342660/1‑42 (MQ=255)
                                                                                                                                                  gACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTc                                                                                                           <  8:342660/42‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAGCGCGC  >  NC_002947/4758280‑4758571

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: