Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,974,081 T→A 7.9% K523M (AAG→ATG) ‡ flgK ← flagellar hook‑associated protein FlgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,974,0810TA 7.9% 219.8 / 6.7 73K523M (AAG→ATG) ‡flgKflagellar hook‑associated protein FlgK
Reads supporting (aligned to +/- strand):  ref base T (32/35);  new base A (3/3);  total (35/38)
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

CCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTGCG  >  NC_002947/4973934‑4974227
                                                                                                                                                   |                                                                                                                                                  
ncGTCGGCATTGAACGACAGGGTAAAGCNGTCATTGGNCGCAGGGCTNCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATACAGAATCGGGTCGCCGGCCGCATCCAACATAGGCACGGGGATGGACAACtgg                                                                                                                                                  >  1:478836/2‑148 (MQ=255)
  gTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttgtt                                                                                                                                                >  2:894802/1‑150 (MQ=255)
     ggCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTc                                                                                                                                             <  2:1408596/150‑1 (MQ=255)
     ggCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTc                                                                                                                                             <  1:644559/150‑1 (MQ=255)
     ggCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTc                                                                                                                                             <  1:1321611/150‑1 (MQ=255)
     ggCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACAAGTTGTc                                                                                                                                             <  2:1030632/150‑1 (MQ=255)
       cATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACAAGTTGTcca                                                                                                                                           <  1:1353901/150‑2 (MQ=255)
             aCGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAg                                                                                                                                     <  1:911777/150‑1 (MQ=255)
                aCAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCt                                                                                                                                           >  2:271162/1‑141 (MQ=255)
                aCAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCt                                                                                                                                           <  1:271162/141‑1 (MQ=255)
                      tAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttgtt                                                                                                                                                <  1:1032254/130‑1 (MQ=255)
                      tAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACttgtt                                                                                                                                                >  2:1032254/1‑130 (MQ=255)
                      tAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGc                                                                                                                            >  2:640397/1‑150 (MQ=255)
                      tAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGc                                                                                                                            <  1:455540/150‑1 (MQ=255)
                              tCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTa                                                                                                                    >  1:203983/1‑150 (MQ=255)
                              tCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTa                                                                                                                    <  1:640397/150‑1 (MQ=255)
                               cATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTAc                                                                                                                   <  2:1127753/150‑1 (MQ=255)
                               cATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTAc                                                                                                                   <  1:554737/150‑1 (MQ=255)
                                aTTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCCCCACCCTGCCGGTAcc                                                                                                                  >  2:483938/1‑150 (MQ=255)
                                aTTGGCCGCAGGGATGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTAcc                                                                                                                  >  2:1059442/1‑150 (MQ=255)
                                       gCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTgg                                                                                                           <  1:894802/150‑1 (MQ=255)
                                          gggCTGCCACCGATTCTGGTCTCGACGGCAAAGGTACGCGGCTTCCTTGTATTATCCAGTATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGCACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTa                                                                                                        <  2:37302/150‑1 (MQ=255)
                                               gCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCtt                                                                                                   <  2:203983/150‑1 (MQ=255)
                                                ccACCGATTGTGGTCTCGGCGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTg                                                                                                  >  1:321118/1‑150 (MQ=255)
                                                ccACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTg                                                                                                  <  2:1104610/150‑1 (MQ=255)
                                                         gtgGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAAc                                                                                         <  1:778122/150‑1 (MQ=255)
                                                             tctcGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCt                                                                                                             >  1:1412235/1‑126 (MQ=255)
                                                             tctcGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCt                                                                                                             <  2:1412235/126‑1 (MQ=255)
                                                              ctcGACGGAAAAGGTTGGTGTTTTTCCTCTACCATCCAGAAGCGGGTCGCCGTCCGCATCTAGCATCGGCACTTCGATGGACAACTTGCTTTCCTTCTCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCtt                                                                                    <  2:275961/150‑1 (MQ=255)
                                                               tcGACGTCAAAGGGACGCGGGTTACGGCTTGCATCCAGAGTCGGGTCTCCGGCCGCCTCCAGGACCGGCACGACTCTGGACAACTAGTTGTCCTGCCCAGGGACAACGCTGACGGTACCAACCTGGGTACCCTTGGCGTCGAACAGCttt                                                                                   >  2:239994/1‑149 (MQ=255)
                                                               tcGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTg                                                                                   >  2:985372/1‑150 (MQ=255)
                                                                 gACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTa                                                                                 >  1:1008149/1‑150 (MQ=255)
                                                                  aCGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAg                                                                                <  2:1202639/150‑1 (MQ=255)
                                                                     gCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCCACCAGCATGGGCACGGCCATGGGCAACATGTTGTCCTGCCAAGGGACCGCGCTGCCAGGACCAATCTGTGTCCCCTTGACGTCGATCAGCTTGTAGccc                                                                             >  2:499618/1‑150 (MQ=255)
                                                                       aaaGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTg                                                                           >  1:948519/1‑150 (MQ=255)
                                                                          ggTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGa                                                                        >  2:1360130/1‑150 (MQ=255)
                                                                           gTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTa                                                                                 <  1:1321902/140‑1 (MQ=255)
                                                                           gTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTa                                                                                 >  2:1321902/1‑140 (MQ=255)
                                                                                 gggTTACCGCTACCATCCAGAATCTGGTCGCCGGCAGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGGCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGGGCCACCg                                                                 >  1:472869/1‑150 (MQ=255)
                                                                                    ttACCGCTACCATCCAGCATCGGGCCGCCGGCAGCAACCCGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTg                                                              <  2:1008149/150‑1 (MQ=255)
                                                                                           tACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGt                                                       >  2:1360335/1‑150 (MQ=255)
                                                                                               aTCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTa                                                                                                                    >  1:578877/1‑85 (MQ=255)
                                                                                               aTCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTa                                                                                                                    <  2:578877/85‑1 (MQ=255)
                                                                                                tCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGaa                                                                                          <  2:408197/110‑1 (MQ=255)
                                                                                                tCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGaa                                                                                          >  1:408197/1‑110 (MQ=255)
                                                                                                tCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTCTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAac                                                  <  1:1360130/150‑1 (MQ=255)
                                                                                                   aGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATATGGGTCCCCGTGTCGTGGAACAGCTTGTAGCTCTGGGAGCCACCGCTGGCCGCGTCGAAcacc                                               >  2:336992/1‑150 (MQ=255)
                                                                                                     aaTCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAAcaccac                                             <  1:1059442/150‑1 (MQ=255)
                                                                                                        tccatcgccgtgccgatgctggatgcggccggcgaTGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGATCCACCGCTGGCAGCGTCGAACACCAcgcg                                          >  2:254664‑M2/36‑150 (MQ=255)
                                                                                                         ccatcgccgtgccgatgctggatgcggccggcgaTGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTg                                                              >  2:1081614‑M2/35‑129 (MQ=255)
                                                                                                         ccatcgccgtgccgatgctggatgcggccggcgaTGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTg                                                              <  1:1081614‑M2/95‑1 (MQ=255)
                                                                                                           atcgccgtgccgatgctggatgcggccggcgaTGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTg                                                              <  1:1291429‑M2/95‑1 (MQ=255)
                                                                                                           atcgccgtgccgatgctggatgcggccggcgaTGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTg                                                              >  2:1291429‑M2/33‑127 (MQ=255)
                                                                                                             CGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCg                                     >  2:948519‑M1/1‑150 (MQ=255)
                                                                                                                 GGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGACACCTCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGGCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCAt                                 >  1:1155382‑M1/1‑150 (MQ=255)
                                                                                                                  gccgCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATc                                >  1:517912/1‑150 (MQ=255)
                                                                                                                  gccgCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTAGAACACCACGCGTACCGGCATc                                >  2:617654/1‑150 (MQ=255)
                                                                                                                    cgCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTgg                                                                                                           <  1:796256/73‑1 (MQ=255)
                                                                                                                    cgCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTgg                                                                                                           >  2:796256/1‑73 (MQ=255)
                                                                                                                     gCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaa                             <  2:321118/150‑1 (MQ=255)
                                                                                                                      cATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGaa                                                                                          >  1:516377/1‑88 (MQ=255)
                                                                                                                      cATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGaa                                                                                          <  2:516377/88‑1 (MQ=255)
                                                                                                                          cAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaatcgcg                        <  1:1103369/150‑1 (MQ=255)
                                                                                                                           aGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaatcgcga                       >  1:1167469/1‑150 (MQ=255)
                                                                                                                            gCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaatcgcgaa                      >  1:1109276/1‑150 (MQ=255)
                                                                                                                             cATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaatcgcgaat                     <  1:617654/150‑1 (MQ=255)
                                                                                                                              aTCGGAACGGAGATGGACAACTTGTTGTCCTGCCCAGGCACCACGATGCCGGTAGCAATCTGGGTCCCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATcgaatcgcgaatc                    <  2:293188/150‑1 (MQ=255)
                                                                                                                                  gCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGa                                                                        >  1:717481/1‑94 (MQ=255)
                                                                                                                                  gCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGa                                                                        <  2:717481/94‑1 (MQ=255)
                                                                                                                                           aggACAACAAGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCa                                                           >  1:25/2‑98 (MQ=255)
                                                                                                                                           aggACAACAAGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCa                                                           <  2:25/97‑1 (MQ=255)
                                                                                                                                           aggACAACAAGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCa       >  2:832261/2‑150 (MQ=255)
                                                                                                                                            ggACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGActgctg            >  1:531680/1‑144 (MQ=255)
                                                                                                                                            ggACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGActgctg            <  2:531680/144‑1 (MQ=255)
                                                                                                                                              aCAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAAcaccac                                             <  2:411071/109‑1 (MQ=255)
                                                                                                                                              aCAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAAcaccac                                             >  1:411071/1‑109 (MQ=255)
                                                                                                                                                aaCTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTGCg  <  2:1109276/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATGGACAACTTGTTGTCCTGCCCAGGCACCACGCTGCCGGTACCAATCTGGGTACCCTTGGCGTCGAACAGCTTGTAGCCCTGGGAGCCACCGCTGGCAGCGTCGAACACCACGCGTACCGGCATCGAATCGCGAATCGACTGCTGGATCAGTGCG  >  NC_002947/4973934‑4974227

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